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A method for measuring the upper limits of soft beta-spectra by use of screen wall tube 
counters filled with hydrogen gas is described. Illustrative results are given for Rb*’, Na®, S*, 


and Au!*, 





INTRODUCTION 


HE necessity of checking ary relation 

between the decay constants of beta- 
radioactive nuclei and the energies of the emitted 
radiations of the type Sargent' has suggested 
makes data on the energies of soft beta-spectra 
valuable. The difficulties of such a measurement, 
arising principally from the high probability of 
scattering for low energy electrons, have pre- 
viously prevented our obtaining of data of 
desirable reliability for the rather numerous 
long-lived soft emitters. 

The following discussion is intended to be a 
brief presentation of a method which apparently 
will serve for beta-radiation harder than about 
ten or twenty kilovolts. Illustrative data are 
given. 


APPARATUS AND METHOD 


The method? employs a screen-wall Geiger- 
Miiller counter placed along the axis of a sur- 
rounding cylinder which has the radioactive 
body mounted on its inner surface. Fig. 1 shows 
the counter. The whole chamber is evacuated 


1B. W. Sargent, Proc. Roy. Soc. A139, 659 (1933). 
*W. F. Libby, Phys. Rev. 46, 196 (1934). 


and filled with hydrogen gas at a few centimeters 
of Hg pressure and is placed in a solenoidal 
magnet with its axis parallel to the field. The 
effect of the sample on the counter (obtained by 
movement of the sample cylinder inside the 
counter chamber) is studied as a function of the 
field and, in general, is found to disappear at a 
certain value of the field. If gamma-radiation is 
present there is a constant residual effect. The 
maximum energy of the beta-spectrum is deter- 
mined by the field at which the effect of the 
sample on the counter is just eliminated. The Hp 
(product of field and radius of curvature of 
electron track) value for the upper limit of the 
spectrum is this limiting field multiplied by half 
of the shortest distance from the sample surface 
to the counter, because any smaller field allows 
some of the curved tracks to penetrate the 
counter. The geometrical arrangement used in 
this method was suggested and used by Boc- 
ciarelli and Occhialini* on the potassium radia- 
tion with a solid-wall counter as detector. 

The amplifying and recording units used with 
the counters are of conventional design for use 
with solid-wall counters. In fact, screen-wall 


3 D. Bocciarelli, Nature 128, 375 (1931). 
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Fic. 1. Screen wall counter assembly. 


counters differ in only one essential respect from 
solid-wall counters insofar as operation is con- 
cerned. It is necessary to provide an electric 
field between the counter and the surrounding 
chamber wall, presumably to prevent the 
accumulation of ions in the intervening space 
some of which would be free to move back into 
the counting volume and trip the counter. 
Under normal operating conditions the chamber 
wall (and sample cylinder) is kept about fifty 
volts positive relative to the counter wall by 
means of an insulated battery. The size of this 
potential drop seems to be rather unimportant 
providing it is sufficiently small to prevent 
sparking across the gap. It is convenient occa- 
sionally to ground the chamber for purposes of 
handling and to allow the full counter voltage 
to serve as a clearing voltage. Fig. 2 shows the 
effects produced by changing the sign and mag- 
nitude of this clearing field for a standard counter 
with diameter of 2.5 cm and sample cylinder 
diameter of 6 cm. The gas was helium at a 
pressure of 6 cm of mercury, with a trace of air 
impurity. The three curves are for the back- 
ground, potassium radiation plus background, 
and aluminum photoelectrons plus background. 
The aluminum photoelectrons were obtained by 
illuminating an aluminum sample cylinder with 
visible light. The rise observed with the chamber 
wall more negative than the counter was due to 
an effective increase in size of the counter. In a 
certain sense, the screen is replaced as the 
counter wall by the chamber wall and sample 
cylinder. Various uses of this arrangement for 
construction of large sensitive counters suggested 
themselves. In particular, a soft radiation in- 


cident at the wall of a counter of large diameter 
might not produce enough ion pairs to initiate 
the counter discharge in the case of the ordinary 
counter. However, in the case of the arrangement 
described above with a sufficiently large negative 
potential on the sample cylinder with respect to 
the screen, the initial ionization will be increased 
more than it would have been in the standard 
counter because of the larger fields near the 
outer wall. The chances of counting are neces- 
sarily increased. Evidence for this appears in the 
figure in the relative effects of changing the sign 
of the clearing field on the potassium beta- 
radiation and the aluminum photoelectrons. The 
effect on the latter is an increase of considerable 
size while that on the potassium-plus-background 
radiation is little more than is observed with the 
background radiation alone. When the clearing 
voltage is removed entirely the counters have a 
tendency to develop high backgrounds which fi- 
nally grow until they prevent the counter’s oper- 
ating at all. This may require only a few seconds 
with small assemblies. 

The magnetic fields used (6000 gauss at most) 
have two effects on the operation of the counter 
which, though small, must be considered. The 
operating voltage is reduced by placing the 


TABLE I. Counter dimensions. 











SAMPLE 
CYL- 
COUNTER INDER 
COUNTER | COUNTER WIRE Diam- | VALUE 
COUNTER DIAMETER | LENGTH | DIAMETER ETER OF p 
NAME (cM) (cm) (mm) (cm) (cM) 
Large 2.5 14 0.22 13.5 | 2.75 
Medium 2.5 14 0.22 10.1 1.90 
Small 2.5 14 0.22 6.1 0.90 
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screen-wall counter in a field; the effect is about 
two percent at 6000 gauss. The explanation of 
this fact is not clear, though it is possible that it 
is due to an increase in probability of recom- 
bination of ion pairs formed in the space between 
the counter wall and the chamber walls before 
the positive ions so formed have time to move 
to the negatively charged counter wall. In the 
absence of the field these ions would tend to 
collect on the screen and reduce the negative 
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Fic. 2. Effects of clearing fields in small counter (6 cm He 
as gas). 


potential. The effect apparently is associated 
with the curving of the paths traversed by the 
ions and the consequent decrease in their chance 
of reaching the walls before recombination occurs. 
The result is that the field allows the screen to 
maintain a higher average negative potential and 
requires a reduction of the applied voltage relative 
to the field-free value. This effect is not observed 
with solid-wall counters. The adjustment is made 
in practice by setting the voltage so the shape of 
the kicks, as shown on an oscillograph, remains 
constant over the whole range of variation of the 
field. 

The other effect of the field is a linear reduction 
of the background as the field is applied which 
has an obvious explanation in the elimination of 
electron secondaries ejected from the chamber 
wall and sample cylinder by background gamma- 
radiation, etc. This effect amounts to a reduction 
of about thirty percent in most cases. The data 
presented later show this effect and are quite 
typical. 


THE MAGNET 


The magnet used was a solenoid with central 
hole six inches in diameter and one foot long. 
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Fic. 3. Relative probabilities of scattering (=5°) and 
production of one ion pair in Hy. 


The coil is cooled by a kerosene circulation 
system. The assembly is of standard type for 
obtaining the fields used in Zeeman effect work. 
The pole pieces were removed and only one of 
the two coils used in any one determination in 
this work. Both coils were calibrated, however 
(by bismuth spiral and flip coil measurements), 
and both used in the course of the work. Maps 
were made of the fields in the cylindrical holes 
which showed them to be homogeneous within 
five percent over a distance of about 20 cm along 
the axis of the solenoid. The variations along the 
perpendicular direction in this region were 
smaller. 
THE COUNTERS 


Counters of the type shown in Fig. 1 with the 
characteristics given in Table I were used. All 
counters were constructed of copper screen with 
fourteen wires (0.010” diameter) to the inch. The 
counter gas was hydrogen at a pressure of about 
4 cm of mercury with a little air, water, or ethy! 
alcohol vapor introduced to improve the counting 
characteristics.‘ The sixth column of Table I 
gives the value of the radius of curvature, p, of 
the electrons which leave the sample surface 
tangentially and hit the counter surface tangen- 
tially, traveling in a plane perpendicular to 
the direction of the magnetic field. This is half 
of the shortest distance between the sample and 
counter surfaces. As given in the table it is 
computed for thin samples. In the cases where 


‘A. Trost, Zeits f. Physik 105, 399 (1937). 
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Fic. 4. Beryllium sample in medium counter. (1.85 
cm H,). Open circle, exposed. Solid circle, unex f 
imax = 5.541 amperes. Hmax = 220440 gauss. (E ?) max 

=410+75. Emax™ 15,000 +7,000 ev. 

Fic. 6. Beryllium sample in medium counter. (5.1 cm H; 
as gas). p= 1.85 cm. imax =5.0+1 amperes. Hmax = 200+40 
gauss. (Hp) max = 370475. Emax™ 12,000 +6,000 ev. 





the samples were of appreciable thickness, the 
proper p value was computed and used. The 
uncertainty of this quantity seems to be about 
0.5 mm. This information was obtained by 
running a given spectrum in counters with 
different p values. The uncertainty of measure- 
ment of the distance from sample to counter is 
of that order. 


~THE STRAGGLING EFFECT AND THEORETICAL 
LIMITATIONS OF THE METHOD 


The problem of fixing the curvature in a mag- 
netic field and recording a low energy electron is 
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Fic. 5. Beryllium sample in small counter. (2.8 cm H; as 


gas). p =0.85 cm. imax = 10+2 amperes. Hmax =400+gauss. 
(Hp) max = 340+70. Emax = 11,000+5,000 ev. 


Fic. 7. Straggling effect. Beryllium sample in medium 
counter, (4.0 cm A as gas). 


difficult because the probability of scattering is 
high. The method used in this research requires 
that the number of collisions along the path from 
the sample to the counter be as smal! as is com- 
patible with the production of one or more ion 
pairs in the counter volume proper. Obvious 
steps in this direction are the reduction of the 
distance from sample to counter and the use of 
hydrogen as counter gas. 

Since the essential quantity involved is the 
ratio of the probability of serious scattering 
(>5°, for example) to the probability of pro- 
duction of an ion pair, this quantity has been 
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Fic. 8. Radioactive phosphorus in large counter. 
(3.8 cm H:+0.2 cm C.H.O ). tmax=65245 amperes. 
Hmax = 26004200. p=2.75 cm. (Hp)max=7150+4550. 
Emax = 1.69+0.08 Mev. 

Fic. 10. RaE in medium counter. (Approximately 4 cm 
H,+trace of C;H,OH). p= 1.90 cm. imax =68+43 amperes. 
Hyax = 27204120. (Hp)max=51504250. Emax =1.16 
+0.06 Mev. 


calculated approximately for hydrogen gas from 
the formula for energy loss given by Bethe,* as 
modified by Bloch’s summation method, and the 
formulae for electronic and nuclear scattering 
given by Mott.® Fig. 3 shows the results as a 
function of the energy of the electrons. This 
curve allows an estimate to be made of the 


5H. Bethe, Handbuch der Physik Vol. 24, No. 1 (Berlin, 
1933), p. 519. 
ONE. 


Mott, Proc. Roy. Soc. 125, 222; 126, 259 (1929). 
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Fic. 9. Mesothorium 2 in medium counter. (Approxi- 
mately 4 cm H;,+trace of C;H;OH). p=1.90 cm. imax = 88 
+5 amperes. Hinax =35204200. (Hp) max = 668024380. 
Emax™ 1.55 +0.07 Mev. 


Fic. 11. RaD in small counter. (Hp) max =750+450. 


Emax = 48,000 +6000 ev. 


energy of the softest beta-radiation for which 
the method will be valid. The ratio of path length 
through the counter to the path length before 
reaching the counter must be at least as large as 
the ratio of the probability of serious scattering 
to the probability of ionization. For the small 
counter used this ratio is about one to twenty; 
the chief uncertainty arises from the estimation 
of the path length through the counter. This 
corresponds to an electron energy of about six 
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Fic. 12. Rubidium nitrate in small counter. (3.6 cm 
H,+0.4 cm A). p=0.9 cm. Open circles, exposed. Solid 


circles, unex -tmax = 36.0+3 amperes. Hmax = 40.0imax 
= 1440+120. (Hp) max = 12964110. Emax = 132,000 
+20,000 ev. 


kilovolts. However, because an error of a factor 
of two may have been made in the estimation of 
the critical ratio and the calculations themselves 
are somewhat approximate, the limit cannot 
safely be set lower than about ten kilovolts. 

To summarize, with the best conditions of 
distance (small) from sample cylinder to counter 
wall, He gas (with a small amount of heavy im- 
purity to improve the counting characteristics) 
and the lowest pressure of gas insuring formation 
of a sufficient number of ion pairs in the counter, 

.it should be possible to measure the upper limits 
of beta-spectra if they are no lower than about 
ten kilovolts. 

These conclusions were tested experimentally 
by means of a source of soft electrons (kindly 
loaned by Professor E. M. McMillan of the 
Physics Department). It was a piece of metallic 
beryllium which had been bombarded for several 
months with deuterons in the cyclotron. The 
radiation consists of negative electrons (deter- 
mined by placing sloping vanes between the 
sample cylinder and the counter and reversing 


TABLE II, Energy of 8-rays. 








ENERGY (OTHER 


ENERGY (THIS RESEARCH) METHODS) 
(Mev) 





ELEMENT (Mev) 
p= 1.69 +0.08 1.69 
MsTh 2 1.55 +0.07 1.60 
RaE 1.16 +0.06 1.16 
RaD 0.048 +0.006 0.046 
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Fic. 13. Potassium nitrate in small counter. (3.6 cm 
H:+0.4 cm A). tmax=110+10 amperes. Hmax =4400 
+400. p=0.85 cm. (Hp) max = 37404340. Emax =725,000 
+100,000 ev. 


the field direction to find the direction of bending 
of the radiation around the vanes.' Runs were 
made in both the small and medium counters. 
Of course, the test for absence of the straggling 
and scattering effect was to find conditions under 
which the product (Hp) for the end point would 
be the same for both counters. Figs. 4, 5 and 6 
show this to be true approximately for He gas 
at pressure of 2 to 5 cm of Hg. Fig. 7 shows that 
the effect of replacing the hydrogen with argon 
is the expected smearing and shifting of the end 
point to higher values. Similar results have been 
obtained with oxygen and air. The (fp) max 
values show the upper limit to be 13,000+5000 
electron volts. 

The mechanism of counting for the electrons 
of maximum energy when the field is set for them 
to hit the edge of the counter may not involve 
the production of ions in the gas by single col- 
lisions so seriously as reflection by the screen into 
the counter volume or possibly the ejection of 
electrons from the screen itself. The data afford 
some evidence for this. It seems difficult to account 
for the observed counting near the end point at 
some of the gas pressures used without invoking 
some mechanism of this sort. Of course such an 
effect makes the effect of scattering less important. 


Cuecks AGAINST KNOWN SPECTRA 


In order to prove that the upper limits deter- 
mined by this method are not seriously in error 
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a number of determinations have been made for 
spectra whose upper limits are obtainable by 
other methods. Figs. 8, 9 and 10 present typical 
determinations for P®, MsTh 2, and RaE. Fig. 
11 for RaD gives the upper limit of 46 kilovolts 
found by the standard semi-circular focusing 
method’: *.* which is applicable in this low 
energy range in this case because the radiation 
is intense. The 46 kilovolt radiation is emitted 
by conversion of a 47.1-kilovolt gamma-ray. The 
semi-circular focusing method which uses a 
photographic plate is inherently a very good 
apparatus for studying low energy spectra 
providing there is sufficient intensity available 
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Fic. 14. Radioactive sulfur (BaSO,) in small counter. 
(3.6 cm H,+0.4 cm C,H,;OH). p=0.9 cm. imax =3243 
amperes. Himax =40.0imax = 12804120. (Hp) max =1152 
+108. Emax = 107,000 + 20,000 ev. 


Tase III. Summary of results. 








ENERGY OF 
Ete- HAtr- Upper Lmirt 
MENT LIFE (Mev) REFERENCE 





Ac 13.4 yr. | 0.220+0.05 | D. E. Hull, W. F. Libb 
and W. M., Latimer, . 

Am. Chem. Soc. 57, 

1649 (1935) 

MsThlI | 6.7 yr. | 0.053+0.005| D. D. Lee and W. F. 

Libby, following article 

s* 88 days} 0.107+0.020) This paper* 

Au'** 2.7 days} 0.77 +0.16 | This paper 

Rb** a yr. | 0.132+0.02 | This paper 

















* The authors are indebted to Professor D. M. Yost of the California 
Institute of Technology for communicating this value of the half-life 
of radioactive sulfur checking that given first by E. B. Andersen, Zeits. 
f. physik. Chemie, 32, 237 (1936). 


7L. Danysz, Le Radium 9, 1 (1911). 
5 L. Meitner, Zeits. f. Physik 11, 35 (1922). 
*D. H. Black, Proc. Roy. Soc. A109, 166 (1925). 
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Fic. 15. Radioactive gold in large counter. (4 cm 
H:+trace of air). Open circles, exposed. Solid circles, 


unexposed. imax =33+4 x toe e=2.70 cm. Hmax 
=41.11+80= 14404165. (Hp)max=38904450. Emax 
= 770,000 + 160,000 ev. 


to give a trace on the plate. The whole space 
between source and plate can be evacuated and 
the scattering by the gas eliminated. There cer- 
tainly is little reason to doubt the correctness of 
its result on the RaD conversion electrons. 

Table II presents the results obtained by the 
present method together with the accepted 
values from other methods. 


RESULTS 


Figure 12 gives data obtained for the natural 
radioactivity of rubidium, due to Rb*’. The 
result is 132,000+ 20,000 ev and is considerably 
lower than the value of 250,000 ev previously 
used.'? 

Figure 13 shows the data for the natural 
activity of potassium, due to K*, which give a 
value of 725,000+ 100,000 ev in agreement with 
the older determinations. No serious attempt was 
made to make this accurate because of the agree- 
ment with the older work and the low intensity 
of the radiation. 

Figure 14 presents data for radioactive sulfur 
(S**) and Fig. 15 for radioactive gold. The result 
on radioactive sulfur is new ; that for radioactive 
gold checks the value given by McMillan, 
Kamen, and Ruben" from absorption. 

Table III is a summary of results obtained by 
this method. 

1° O. Klemperer, Proc. Roy. Soc. 148, 638 (1935). 


" E,. McMillan, M. Kamen and S. Ruben, Phys. Rev. 
52, 375 (1937). 














(1907); B. Boltwood, Physik. Zeits. 8, 556 (1907); H. N. 
McCoy and W. Ross, J. Am. Chem. Soc. 29, 1709 (1907); 
L. Haitinger, K. Peters and St. Meyer, Vienna, Akad. 
Sitz.-Ber. 120 [2a], 1199 (1911); L. Meitner, Physik. Zeits. 
19, 257 (1918). 

*H. O. W. Richardson and A. Leigh-Smith, Proc. Roy. 
Soc. A160, 454 (1937). 








Gray and O'Leary, Nature 123, 568 (1929); L. H. Gray, 
Nature 130, 738 (1932); J. Petrova, Zeits. f. Physik 55, 628 
(1929) ;S. Bramson, Zeits. f. Physik 66, 721 (1930) ; E. Stahel 
and G. J. Sizoo, Zeits. f. Physik 66, 741 (1930); H. O. W. 
Richardson, Proc. Roy. Soc. A133, 367 (1931); G. von 
Droste zu Vischering, Zeits. f. Physik 84, 17 (1933); E. 
Stahel, Helv. Phys. Acta 8, 651 (1935). 
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; FEBRUARY 1, 1939 PHYSICAL REVIEW VOLUME 55 
| . The Beta-Rays of Mesothorium 1 and Radium D , 
D. D. Lee anp W. F. Lipsy . 
Department of Chemistry, University of California, Berkeley, California I 
(Received September 6, 1938) 4 
s 
Mesothorium 1 has been prepared from an old sample tinuous spectrum was found to be 5150+250 gauss cm, c 
of thorium sulfate whose purity was demonstrated by The upper limit apparently to be assigned to the primary ( 
measuring the rate of decay of thorium B prepared there- _beta-ray spectrum of RaD was determined as 546+10 e 
from, The mesothorium 1 was shown to be free from con- gauss cm, or 25,500+1000 ev, in agreement with the de- 
taminations of other thorium decay products and isotopic termination of Richardson and Leigh-Smith. An estimate : 
radioactive elements by measuring the rate of growth of of the observed number of secondary beta-rays arising from c 
mesothorium 2, and proving the absence of alpha-emitters. the conversion of the 47.2-kev gamma-rays of RaD indi- f 
Magnetic deflection and aluminum absorption experi- cates the emission of five to ten such electrons per 100 r 
ments have shown that the particles emitted are electrons _disintegrations, in rough agreement with the results of g 
whose maximum energy is 53,000+4000 electron volts. Kikuchi and the predictions of Fisk for the L conversion 
Samples of RaD and RaE, pure to within 0.3 percent, have _ of dipole radiation, but in complete disagreement with the : 
been prepared. The momentum of the most energetic conclusions of Stahel and von Droste. Reasons for this $ 
conversion beta-ray of RaD has been determined as disagreement are given. l 
765+32 gauss cm, and the upper limit of the RaE con- t 
c 
MONG the naturally occurring radioactive number of electrons having energies below four a 
bodies about which there is very little, or kv, and estimated the upper limit of the spectrum, a 
only controversial, information available at pres- by comparison with the distributions predicted . 
ent are mesothorium 1 and radium D. in the beta-decay theories of Fermi, and of t 
The emission of ionizing radiation from the Konopinski and Uhlenbeck, as either 16, or t 
first of these has never been definitely established, 24, kv. I 
although it was discovered in 1907, and it was The energy of the gamma-radiation emitted a 
for a long time considered to be ‘‘rayless."" How- has been determined by absorption measure- f 
ever, from the chemical properties of meso- ments, and by analysis of the beta-ray line spec- t 
thorium 1! and its product, mesothorium 2 (which trum, the latter leading to the most accurate t 
have been found to be isotopes of radium, and value, 47,200 ev. The more recent of such de- t 
actinium, respectively) and the displacement law, terminations have been done in such a manner t 
it is assumed that a weak beta-particle isemitted as to determine not only the energy of this g 
in the decay of mesothorium 1. gamma-ray, but the energy of the accompanying I 
The decay of radium D to form radium E is x-radiation and to estimate the quantities of s 
known to occur by the emission of primary beta- each emitted. The x-rays have been found to V 
cays, Gare cays, secondary beta-rays from the _ have energies corresponding rather closely to the g 
internal venebteneene of gamma-rays, x-rays, and 7, x-rays of bismuth, with which they have been t 
perhaps tertiary, or Auger, electrons from the identified, although determination of their ener- c 
internal conversion of the x-rays. t 
The first definite information about the pri- “ Danysz, Le Rad. 2, 1, 6 (1912); E. Rutherford and : 
beta- rticles of radium D was furnished . Richardson, il. ag. 26, 324, 937 (1913); C. D. Ellis, 
ee ee ~ was sts. «., Proc. Camb. Phil. Soc. 21, 121 (1921-23); L. Meitner, | 
by a very recent cloud-chamber investigation by Zeit. f. Physik 9, 131, 145 (1922); 11, 35 (1922); I. Curie \ 
Richardson and Leigh-Smith,? who found a large and G. Fournier, Comptes rendus 176, 1301 (1923); D. H. f 
Black, Proc. Roy. Soc. A109, 166 (1925); L. F. Curtiss, 
1. git? ADT Phys. Rev. (227, 672 (1926). C 
10. Hahn, Ber. 40, 1462 (1907), Physik. Zeits. 8, 277 *S. Kikuchi, Jap. J. of Phys. 4, 143 (1925-27); L. H. t 
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gies have given only average values. The number 
of such photons has been found to be 22 to 30 
per 100 disintegrations, and the number of 
47,000-ev gamma-rays, about 3.5 per 100.5 We 
shall assume that the x-rays are excited by the 
conversion of 47,000-ev gamma-rays in L shells, 
(the K shell could not be excited because the 
energy of the gamma-ray is too low). The number 
observed, considered with the fluorescence yield 
of bismuth, has led to the conclusion that each 
primary beta emitted leaves the resulting RaE 
nucleus in an excited state 47 kv above the 
ground state of RaE, and that all such nuclei 
emit gamma-rays of this energy. Of these, 96 or 
97 percent are internally converted in L, M and 
N shells of extranuclear electrons, the remaining 
three or four percent are emitted without 
conversion, 

In a series of experiments, wherein were used 
a point counter covered with a collodion foil and 
a source of RaD, RaE and Po placed in an 
evacuated chamber, Stahel® measured the num- 
ber of electrons from RaD. He distinguished 
these from the beta-rays of RaE, the alphas of 
Po, and the general background radiation by 
means of a magnetic field and various aluminum 
foils. The absorption of the RaD betas in one, 
two and three of his collodion foils led Stahel to 
the conclusion that he was observing, through 
the first permanent foil, only the secondary elec- 
trons due to the conversion of the 47,200-ev 
gamma-ray; but, as it was pointed out by 
Feather and Richardson,’ the absorption which 
Stahel observed did not agree with the absorption 
which would be calculated from the known ener- 
gies and intensities of these electrons, and the 
range-energy relations of Williams.* In fact, the 
observed absorption could be accounted for only 
by the assumption that Stahel was getting a 
relatively large number of low energy (4 to 20 
kv) electrons through the first foil. From the 
work of Richardson and Leigh-Smith,* and the 
fact that Stahel’s first foil could completely stop 
only electrons of less than 4500 ev, it is clear 
that more than half of these might have been 


*H. Lay, Zeits. f. Physik 91, 533 (1934). 

* E. Stahel, Zeits. f. Physik 68, 1 (1931). 

7N. Feather and H. O. W. Richardson, Nature 129, 314 
(1932). 

*E. J. Williams, Proc. Roy. Soc. A130, 310 (1931). 





primary beta-rays. From this fact alone and 
without estimating the number of observed elec- 
trons which might have been caused by Auger 
effect of the x-radiation, one is led to conclude 
that there are not 83 or more gamma-ray con- 
version electrons per 100 disintegrations as 
Stahel believed, but that this number must be 
considerably smaller, and that the number of 
47,200-ev gamma-rays (emitted plus converted) 
is less than one per disintegration. It follows from 
this that since no beta-rays of RaD having ener- 
gies corresponding to a transition from ground 
state to ground state have been observed, there 
may occur the emission of two or more gamma- 
rays of less than 47,200 ev energy (nuclear 
x-rays) in a “‘cascade”’ process involving more 
than the one known excited level of the RaE 
nucleus. It is not at all improbable, that such 
nuclear x-rays may have been observed by 
Kikuchi, Gray and O'Leary, von Droste and 
Stahel, and their corresponding conversion elec- 
trons may have been observed by Kikuchi, 
Petrova, Richardson, and Richardson and Leigh- 
Smith, only to be confused with the L rays and 
Auger electrons arising therefrom. 

From the work of Kikuchi one finds only about 
twenty electrons attributable to the conversion 
of 47,200-ev gamma-rays per 100 disintegrations. 
Since about thirteen of these are from the L 
orbits, the Z conversion coefficient calculated 
from this value is 13/3.5, or about 3.7. A theo- 
retical calculation by Fisk® gives the values 1.8 
and 2.9 for the L conversion coefficients of dipole 
and quadrupole radiation, respectively. Only the 
assumption of multipole, or magnetic multipole 
radiation could have accounted for the coefficient 
17.1 calculated by Stahel. 

The subsequently described experiments with 
RaD were intended to furnish information re- 
garding the upper limit of the primary beta-rays, 
and to permit an estimate of the number of 
secondary electrons produced by the conversion 
of 47,200-ev gamma-rays. 


PREPARATION OF SAMPLES 


All mesothorium preparations were made from 
a thirteen-year-old sample of Kahlbaum thorium 
sulfate, the purity of which was determined by 


* J. B. Fisk, Proc. Roy. Soc. 143, 674 (1934). 
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Fic. 1. Growth and ap | curve for ThB precipitate to test 
purity of thorium sample. 


following the growth and subsequent decay of 
gamma- and strong beta-ray activity in a sample 
of thorium B prepared from about one-half gram 
of thorium sulfate. The activity was measured 
with a solid-wall Geiger-Miiller counter and 
increased in the first four hours after at exactly 
the rate to be expected for the growth of ThC 
from ThB preparation. After six hours, the 
activity decayed with the half-life of thorium B. 
This indicates the absence of more than three to 
five percent of any radium B, actinium B, or 
radium D contamination (see Fig. 1). 

Mesothorium 1 and 2 were removed from a 
solution of thorium sulfate by precipitating 
barium sulfate. To remove from this precipitate 
any traces of occluded thorium or radiothorium, 
the BaSOQ, was converted to BaCl, by heating 
with carbon, digestion of the residue in concen- 
trated sodium carbonate solution, and leaching 
of the resulting solid with dilute hydrochloric 
acid. To the resulting solution a trace of ferric 
chloride was added, and the ferric ion removed as 
ferric hydroxide by the addition of a slight excess 
of sodium hydroxide solution. This procedure 
removed thorium, radiothorium, mesothorium 2, 
and the elements of the thorium active deposit. 
As thorium X is isotopic with mesothorium |, it 
cannot be removed by chemical separations but 
must be allowed to decay. About one month is 
necessary for the removal of 99.6 percent of ThX 
and its decay products. 

From such a month-old sample of meso- 
thorium 1 and 2, either of these two elements 
could be removed by simple chemical precipita- 
tions. The addition of ferric ion as the chloride 
and its removal as ferric hydroxide leaves a solu- 





LIBBY 


tion containing only about ten percent of its 
original content of mesothorium 2. Two further 
treatments of this kind suffice to remove all 
noticeable amounts of mesothorium 2 and the 
precipitation of the barium ion as barium car- 
bonate removes mesothorium 1 for experimenta- 
tion. The growth of activity of such a sample, as 
observed with a screen wall counter, is shown in 
Fig. 2. Comparison of the measured activities 
with the calculated solid curve shows the sample 
to have been pure to within at least two percent. 
The initial activity is not proved to be caused by 
MsTh 1 beta-rays since the same curve would 
have been obtained with a MsTh 1 sample con- 
taining a trace of MsTh 2. Counting of alphas 
possibly emitted by an equivalent sample of 
MsTh 1 and 2 showed certainly less than two 
percent contamination by radium, thorium X, or 
actinium X. 

The ferric hydroxide precipitates containing 
mesothorium 2 contain also appreciable quan- 
tities of mesothorium 1. About eight to ten per- 
cent has been found in all such precipitates ex- 
amined in this work by following the decay of 
activity with an electroscope. However, if such 
a sample be dissolved in dilute hydrochloric acid, 
the hydrogen ion concentration adjusted to about 
10-* M by the addition of sodium acetate, and a 
large excess of sodium dichromate with a trace of 
barium chloride added, the resulting precipitate 
removes enough mesothorium 1 so that the ferric 
hydroxide precipitate which can be prepared 
from the filtrate will contain mesothorium 2 pure 
to within 98 percent (see Fig. 3). 
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Fic. 2. Growth of mesothorium 1 precipitate. 
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Samples of radium D and radium E were pre- 
pared from the contents of an old radon tube, 
which was washed in concentrated nitric acid and 
water, and crushed under about five cc of concen- 
trated HNO;. After heating and allowing the 
glass particles to settle, the solution was decanted 
and diluted by the addition of about three cc of 
distilled water. Individual preparations were 
made from one-cc portions of this final solution 
by the usual methods, i.e., the solution was 
evaporated twice with excess concentrated HCl, 
the residue dissolved in about five cc of 0.1N HCl, 
and heated to 90° or 100°C. Then 99 percent of 
the polonium content was removed by rotating 
in the solution at about 200 r.p.m. a silver foil (of 
approximately one cm? area) for one hour, and 
the remaining polonium removed by the same 
treatment with another silver foil. Radium E was 
removed similarly on a nickel foil. Such a sample 
of RaE can be shown by following the decay of 
activity for 35 days to have been originally free 
from Po to at least 99.9 percent and to have con- 
tained no more than 0.25 percent of RaD. This 
was done by using the screen wall electroscope 
to be described later and various aluminum foils 
in order to distinguish between the beta-rays of 
RaE and the alphas of Po. 

After the above-mentioned removal of RaE, 
0.1 cc of a solution containing 10~* gram/cc each 
of Pb and Bi, as chlorides, were added to prevent 
the loss of RaE or RaD by occlusion on dust 
particles, etc. the solution was allowed to stand 
over night. Then the RaE left by the previous 
treatment and that produced subsequently by 
the decay of RaD was removed by rotating in the 
hot solution a nickel foil of the same size as that 
used in the previous removal of RaE. The result- 
ing residue was evaporated once again with excess 
HNOs, then twice with distilled water, and the 
final residue digested with five cc of a two percent 
solution of NH,NO; to insure the solution of the 
lead content. The BiONO; residue was filtered 
off to remove the last traces of RaE. Radium D 
was removed from the filtrate by precipitating 
the lead as PbS and found to be radioactively 
pure to at least 99.8 percent. 


APPARATUS AND EXPERIMENTAL TECHNIQUE 


The counter tubes, voltage supply, amplifier, 
recorder and magnet used in these magnetic 
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deflection experiments were the same as those 
described in the preceding article. 

In determining the absorption of the meso- 
thorium 1 and 2 beta-rays by aluminum, meas- 
urements were made with a Lauritsen electroscope 
having a portion of its cylindrical case cut out 
and replaced by a piece of copper gauze such as 
is used for the screen of a screen-wall counter. 
In order to determine how far an ionizing particle 
must travel into the space enclosed by this 
modified case to be detected, measurements were 
made of the rate of discharge of the electroscope 
due to the alpha-particles from a Po source 
placed at various distances from the screened 
portion. The results showed that any alpha- 
particle, and presumably any other ionizing 
radiation which passes into the screen-enclosed 
space produces ions leading to the discharge of 
the electroscope. 

Corrections for the growth of mesothorium 2 
in samples of mesothorium 1 which were used to 
determine the aluminum absorption and the upper 
energy limit of the beta-ray spectrum were made 
in the following manner. The activity observed 
on the electroscope for the uncovered sample in 
a standard position was measured first. Then 
absorbers were added in increasing thicknesses, 
as measurements of the penetration were made 
until the absorber thickness reached 7 or 12 
mg/cm?*. Following this procedure the absorber 
thickness was reduced by the same increments 
until the sample was uncovered. The average of 
the two readings taken at each absorber thickness 
automatically gave a value of the activity cor- 
rected for growth of mesothorium. This average 
is nearly the activity which would have been 
obtained had it been possible to make all 
measurements simultaneously at the time of the 
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Fic. 4. Aluminum absorption curve for mesothorium 1 and 
mesothorium 2, 


maximum thickness measurement, since the 
growth of mesothorium 2 is a nearly linear 
function of time within any period of two hours 
or less, and the two readings at each thickness 
are time-symmetrical with respect to the maxi- 
mum thickness reading. To find the activities 
caused by mesothorium 1 alone, this curve was 
compared with the absorption curve for pure 
mesothorium 2 and the portion of the former 
curve which was parallel to the latter was 
extrapolated back. Subtracting the activities on 
this extrapolated curve from those on the meas- 
ured curve gives a third curve (see Fig. 4) which 
is the absorption of the beta-rays of mesothorium 
1 in aluminum. 

The same technique was used in magnetic 
deflection experiments on mesothorium. The 
activities from a pure mesothorium 2 sample were 
measured at various field strengths and com- 
pared with the curve obtained when a sample 
originally containing only mesothorium 1 was 
measured at various fields. The field was in- 
creased to some maximum strength by pre- 
determined increments, then decreased by the 
same steps. The difference between the measured 
curve and the extrapolated mesothorium 2 curve 
shows the effect of field on the beta-rays of 
mesothorium 1 alone (see Fig. 5). 

In the case of the magnetic deflection experi- 
ments on radium D no such corrections were 
necessary provided the measurements were made 
within a few hours after preparation of the 
samples because of the slower growth of ra- 
dium E, 


EXPERIMENTAL RESULTS 


The average of two determinations of the 
mesothorium 1 upper limit by absorption (see 
Fig. 5) is 57,000+5000 ev, while five magnetic 
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deflection experiments gave an average of 771 
+45 Hp, or 51,000+4700 ev (see Fig. 6). 
Because of the greater accuracy inherent in the 
second type of measurement, these values are 
taken as indicating that the upper limit of the 
mesothorium 1 spectrum is 53,000+4000 ey. 
That this is the upper limit of the primary beta- 
ray spectrum does not necessarily follow, for 
although no gamma-rays have been observed in 
the disintegration of mesothorium 1, the present 
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Fic. 5. Magnetic determination of mesothorium 1. 
Upper limit. 


work does not exclude the possibility that soft 
gamma-rays might exist and their conversion 
give rise to secondary electrons having energies 
up to the observed maximum. This might make 
the observed upper limit somewhat higher than 
the true value. However, since the total con- 
version of soft gammas is relatively improbable 
because it must occur in the Z shell in this case, 
it seems that most of the electrons are of nuclear 
origin. This point is illustrated in the work 
on RaD. 

Figure 7 shows the activities observed at 
various magnetic field strengths with a sample 
of RaD. Curve (1) gives activities due to about 
one-half of the sample measured at various 
fields up to 1500 gauss; curve (II) gives the 
corresponding values for a smaller portion of the 
same sample. (Electroscope measurements made 
19 days later to determine the relative activities 
of these portions of the sample due to accumu- 
lated RaE in each indicated a ratio of 21 to 1. 
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As the activity of the small portion at zero field 
is observed to be 1200 counts per minute, that 
of the large portion must have been over 25,000, 
a value too large to be recorded by the equip- 
ment used.) 

Examination of curve (I) reveals an activity 
which falls off very rapidly with increasing field 
strength up to an Hp value of 570 gauss cm, a 
more gradual slope from that value to 770 gauss 
cm, and a very gradual slope from 770 to 1500 
gauss cm. Curve (II), on the other hand, shows 
only the steep portion ending at 530+20 gauss 
cm, which is apparently the upper limit of the 
primary beta-ray spectrum of RaD, since all the 
betas from the conversion of the 47,200-ev 
gamma-rays have larger values of Hp. The 
gradual slope in curve (I) is, therefore, pre- 
sumably due to conversion betas, and gives an 
upper limit of 770+40 gauss cm which is not 
far from the accepted value,’® 738 gauss cm for 
the highest conversion line. The remaining 
activity in the large sample is undoubtedly due 
to RaE betas, and perhaps Po alphas, as well 
as RaD gammas. Similar curves obtained with 
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Fic. 6. Average curve for magnetic determination of 
mesothorium 1, Upper limit. 


RaD samples prepared in the same or a slightly 
different manner, gave the values collected in 
Table I. 

To obtain an independent check on any 
possible systematic errors in the magnetic depth 
determinations of such upper limits, measure- 
ments were made with a sample of mesothorium 
2 and with one of RaE. The observed upper limit 
for mesothorium 2 was 6880+250 Hp, in ex- 


1 L. Meitner, reference 3. 
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Fic. 7. Magnetic analysis of RaD. 


cellent agreement with the magnetic spectro- 
graphic determination of Yovanovitch and 
d’Espine." For RaE, the upper limit was 
5150+250 Hp, a value entirely concordant with 
the latest determinations. 


DISCUSSION OF RESULTS 


From the data of Table I one is led to conclude 
that the end point observed with RaD at 
546+10 gauss cm, or 25,500+1000 ev, is that 
of the continuous primary beta-ray spectrum, 
although it is possible that the conversion of 
nuclear x-rays, if they exist in this disintegration, 
might give rise to electrons of energies up to this 
limit. 

The electrons found at Hp values of 530 to 


TABLE I. Data on RaD. 








Uprer Limit 


(Steer Portion HiGHest CONVER- 


Torta. INTENSITY 

















oF SAMPLE or CURVE SION LINE 
(COUNTS/ MIN) (GAUSS an (GAUSS CM) 
1200 530 Not observed 
1500 540 Not observed 
1600 540 Not observed 
Av. 1433 Av. 537 
15,000 540 Not determined 
20,000 570 Not determined 
25,000 570 770 
35,000 530 760 
Av. 23,750 Av. 553 Av. 765 








uD. K. Yovanovitch and J. d’Espine, Comptes rendus 
179, 1162 (1924); J. phys. radium, [6 8, 276 (1927). 
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770 gauss cm are then to be considered the only 
observed beta-particles from the conversion of 
the 47,200 ev gamma-ray ; and by extrapolating 
this portion of the curves to zero magnetic field 
(under the assumption of the theoretical appear- 
ance of the curve which would be expected with 
this apparatus for a source emitting only the 
gamma-ray conversion-betas of RaD) one finds 


AND H. 


R. CRANE 


about five such conversion electrons per 100 
disintegrations. Since this extrapolation is some- 
what difficult, it seems entirely possible that 
there might be as many as ten such conversion 
betas per 100 disintegrations in which case the 
calculated L conversion coefficient would be 2.1, 
very nearly the value calculated by Fisk for the 
case of dipole radiation. 
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The Absorption Coefficient of 5.8-Mev Gamma-Radiation in Aluminum 


J. HALPERN AND H. R. Crane 
University of Michigan, Ann Arbor, Michigan 
(Received December 7, 1938) 


The absorption coefficient of monochromatic 5.8-Mev radiation in aluminum is found to be 
0.075 +0.009 cm™'. By deducting the theoretical value of the pair absorption coefficient, the 
Klein-Nishina absorption coefficient is found, which is 0.062 +0.009 cm~. This value is in good 
agreement with the theoretical value. The experimental method used is such that all quanta 
which are scattered more than 10 degrees in the aluminum are excluded completely from 


the data. 


EVERAL measurements’ of the absorption 
coeficient of the F+H!' gamma-radiation 
have been made by means of the ionization 
chamber. The experiment to be described gives 
a result by a different method, which does not 
involve the assumption that the radiation con- 
sists of a single line and which reduces the effect 
of scattering and geometrical errors to a very 
small amount. The method was first used by 
Delsasso, Fowler and Lauritsen? for the measure- 
ment of the absorption coefficient of the Li+H’ 
and the F+H! radiation in lead. 
Our experimental arrangement is shown in Fig. 
1. On alternate expansions of the cloud chamber 
a block of aluminum 10.2 cm thick was placed 
in the path of the gamma-rays, between the 
target and the lead collimator. A plot of the 
number versus energy of the recoil electrons 
ejected from the carbon sheet was made for each 
of the two cases: with and without the aluminum 
block in the gamma-ray path. These are shown 


1E,. McMillan, Phys. Rev. 46, 325 (1934); H. R. Crane, 
L. A. Delsasso, W. A. Fowler and C. C. Lauritsen, Phys. 
Rev. 46, 1109 (1934). 

*L. A. Delsasso, W. A. Fowler and C. C. Lauritsen, 
Phys. Rev. 51, 391 (1937); Phys. Rev. 51, 527 (1937). 


in Fig. 2. If, in both curves, we consider only 
tracks in the shaded area (above 4.6 Mev) we 
can be sure that the gamma-ray quanta re- 
sponsible for these have not been scattered more 
than 10 degrees before reaching the carbon sheet. 
It is calculated from the angle-energy relations 
in the Compton effect that if a quantum of 
initial energy 5.8 Mev were scattered more than 
that amount in the aluminum, it could not 
produce a recoil electron in the carbon sheet with 
an energy greater than 4.6 Mev, and hence could 
not fall within the shaded area of the curve. 
Similarly, gamma-ray quanta of initial energy 
less than 4.85 Mev are automatically excluded, 
because they cannot give rise to recoil electrons 
of more than 4.6 Mev. 





NS 


Zao / 


ML khhd 





Yj 


X WY 





YY” 1/14) 


yy ALUMINUM ° 
L 102 x10x5 4 


Vdd 


SO 


CHANNEL 





Yer Y 
YY Lean 


Zi YOM 


- 


TARGET 


























CLOUD CHAMBER 


Fic. 1. Experimental arrangement. 
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Fic. 2. Energy distribution of recoil electrons ejected 
from carbon. 


By cutting off the distribution in Fig. 2 at a 
higher energy than that shown, the maximum 
allowed angle of scattering can be made as small 
as desired, but on the other hand the uncertainty 
in the result becomes greater due to the smaller 
number of tracks included. By cutting off at 4.6 


Mey, as is done in this computation, we obtain 
results which are equivalent, as far as geometry 
is concerned, to the usual ionization experiment 
in which the maximum angle subtended at the 
ionization chamber by the absorber is only 10 
degrees. 

The numbers of tracks in the two shaded areas 
are 207 and 98. These give a value (98+6.6)/ 
(207 +9.6) =0.47+0.04 for the attenuation of 
the radiation in 10.2 cm of aluminum. The total 
absorption coefficient is given by writing 
0.47+0.04=e—"-** which gives 4=0.075+0.009 
cm". The pair absorption coefficient is, according 
to Bethe and Heitler, 0.013 cm for aluminum. 
The photoelectric absorption is negligibly small. 
The coefficient due to Compton scattering alone 
(Klein-Nishina absorption coefficient) therefore 
becomes 0.062+0.009 cm~'. The Klein-Nishina 
theory gives 0.060 cm™. 

It is interesting to note that there are rela- 
tively more low energy electrons in the lower 
curve than in the upper curve in Fig. 2. This is 
undoubtedly due to radiation softened by scat- 
tering in the aluminum, and it would have 
vitiated the results, had we been measuring the 
intensities with an ionization chamber. 

The authors are indebted to the Rackham 
Endowment Fund for the support of this work. 
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The Pair Internal Conversion Coefficient in the F+-H’ Reaction and Measurements 
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on the Gamma-Ray Spectrum 


J. HALPERN AND H. R. Crane 
University of Michigan, Ann Arbor, Michigan 
(Received December 6, 1938) 


The gamma-radiation from F+H’ was studied by the cloud-chamber method and was 
found to consist of a single line at 5.8 Mev. The group of electrons at about 4 Mev found by 
Gaerttner and Crane was confirmed, and it is shown that this group probably consists of the 
negative members of internal conversion pairs. The coefficient for pair internal conversion was 
measured and found to be about 1 pair per 100 gamma-ray quanta, which is greater than the 
theoretical value. A search was made for gamma-radiation down to 0.5 Mev, and for delayed 
emission of electrons by the CaF; after bombardment, with negative results in each case. 





T is a well established fact that fluorine bom- 
barded with protons emits a strong gamma- 
ray line at 5.8 Mev. There has also been some 
indication of a lower line, at about 3.8 Mev, but 
the evidence on this point is conflicting. As a 
result of the experiments to be described here, 
it has become clear that the coefficient for the 
internal conversion of the 5.8 Mev gamma-ray 
into pairs, in the field of the emitting nucleus, is 
unexpectedly high, and that it was this process 
which gave rise to electrons which were inter- 
preted as indicating a second gamma-ray line. 
We shall give the evidence for this and also 
arrive at a value for the interval conversion 
coefficient. 


CONFIRMATION OF PREVIOUS RESULTS 


The cloud-chamber experiment of Delsasso, 
Fowler and Lauritsen' was repeated with the 
arrangement shown in Fig. 1, which is essentially 
the same as theirs. We used a 6” diameter cloud 
chamber* with a magnetic field of 1850 gauss. 


CHANNEL 


1CMx5SCM 


So } 
TARGET ee 


10. cw —a 


Fic. 1. Experimental arrangement. 








CLOUD CHAMBER 


1L. A. Delsasso, W. A. Fowler and C. C. Lauritsen, 
Phys. Rev. 51, 527 (1937). 
2H. R. Crane, Rev. Sci. Inst. 8, 440 (1937). 


The sheet of scattering material placed across 
the center of the chamber was carbon, 1.5 mm 
thick. The target (CaF:) was bombarded with 
protons of approximately 0.6 Mev energy from 
the high voltage accelerating tube.* The tracks 
of the electrons ejected from the carbon sheet by 
the gamma-rays were measured, and the results 





° 
o 
N 
' 
: F+H 
= | 1.5 MM CARBON 
780 TRACKS 


NUMBER OF TRACKS 
100 


50 














ELCTRON ENERGY IN MEV 


Fic. 2. Energy spectrum of recoil electrons from 1.5 mm 
carbon, obtained with the arrangement shown in Fig. 1. 
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e- TRACKS 


A 


Fic. 3. A Diagram showing the target inside the cloud 
chamber. The cylindrical piece of carbon, 1.5 mm thick, 
yields recoil electrons. B The carbon is supported on a shaft 
which can be rotated from outside the chamber. It is 
placed in the two positions shown in B on alternate expan- 
sions of the chamber. 





*H. R. Crane, Phys. Rev. 52, 11 (1937). 
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are shown in Fig. 2. The distribution clearly 
indicates a single gamma-ray line, and is in 
complete agreement with the results of Delsasso, 
Fowler and Lauritsen. 

The experiment of Gaerttner and Crane,‘ 
which indicated the presence of a lower line, was 
next repeated. The experimental set-up is indi- 
cated as A in Fig. 3; the essential feature is that 
the target is inside the cloud chamber, and is 
immediately surrounded by the scattering ma- 
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Fic. 4. Energy spectrum of the negative electrons obtained 
with the arrangement shown in Fig. 3 A. 


terial, 1.5 mm of carbon. The distribution of 
negative electrons found is shown in Fig. 4, and 
confirms that found by Gaerttner and Crane. 

A small number of tracks obtained by 
Gaerttner and Crane with the arrangement of 
Fig. 1 seemed also to indicate a second line, but 
this seems to have been incorrect in view of the 
distribution shown in Fig. 2, which contains a 
much greater number of tracks. 


New EXPERIMENTS 

To determine whether the lower energy peak 
in Fig. 4 is composed of electrons originating in 
the carbon scatterer or in the CaF; target itself, 
the following experiment was performed. In an 
arrangement shown as B in Fig. 3, the carbon 
scatterer was attached to a moveable arm so 
that it could be placed close to the target on 
alternate expansions of the cloud chamber, and 
could be removed on the other expansions. When 
the carbon was removed, there remained only a 
thin aluminum wall between the target and the 
chamber. Both positive and negative electron 


*E. R. Gaerttner and H. R. Crane, Phys. Rev. 52, 582 


(1937). 
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Fic. 5. Positive and negative electrons obtained with the 
arrangement shown in Fig. 3 B. Equal numbers of photo- 
graphs were taken with the carbon in the two positions 
shown. 
tracks are recorded in Fig. 5. The higher energy 
peak in the negative electron distribution is 
decreased by nearly a factor of two by removal 
of the carbon. The lower energy peak does not 
seem to be decreased at all. Our conclusion is 
that the higher energy peak is composed of 
Compton recoil electrons and is therefore de- 
pendent upon the amount of material surround- 
ing the target; that the lower energy peak is 
composed of electrons which have originated in 
the CaF, target and have passed through the 
carbon scatterer. 

A glance at the plots of positrons will indicate 
two things: that the number of them is insen- 
sitive to the amount of material surrounding the 
target, and that the plot formed by them is not 
significantly different from that of the negatives 
comprising the lower energy peak, C or D. This 
approximate equality of positives and negatives 
makes it unnecessary to postulate any other 
process than pair formation to account for the 
effect in question. Gaerttner and Crane have 
already demonstrated that pair formation by 
internal conversion takes place in the F+H’ 
reaction. The only item which does not seem to 
fit this explanation is the somewhat excessive 
sharpness of the lower energy peak (especially 
in the data of Gaerttner and Crane) as compared 
with that expected from the theoretical distri- 
bution of energy between members of pairs. 











262 J. 
8 ' 
® pai ol 

; g 174 TRACKS 

80 

fe a 








° as 2 3 4 6 , 6 
ELECTRON ENERGY IN MEV 


Fic. 6. Energy spectrum of negative electrons obtained 
with the arrangement shown in Fig. 1, except that a 0.05- 
mm lead scatterer was used instead of the carbon scatterer. 


This however, is too uncertain to cause serious 
concern. 

As a further check upon the notion that the 
lower energy peak is due to negative members 
of pairs, a series of photographs were taken under 
the arrangement of Fig. 1, but using a 0.05-mm 
lead sheet instead of carbon. The negative elec- 
trons appearing in this experiment were a 
mixture of Compton electron and pair members, 
but the pairs were formed in the lead instead of 
in the target as before. The plot of negative 
electrons, given by Fig. 6, shows that the 
negative members of pairs give rise to a second 
peak in the distribution in this case also. 


THE COEFFICIENT FOR PAIR 
INTERNAL CONVERSION 


By comparing the curves in Fig. 5 the ratio 
of the number of pairs to the number of gamma- 
rays can be determined. The area E+A in Fig. 5 
represents approximately those electrons due to 
Klein-Nishina absorption in the 1.5-mm carbon 
absorber plus the other material surrounding the 
gamma-ray source. The difference in areas, 
(E+A) minus (F+B) gives the absorption in 
1.5 mm carbon alone. By using the Klein- 
Nishina absorption coefficient, the number of 
gamma-rays emitted by the source (in the interval 
of time and solid angle determined by the experi- 
mental arrangement) is found. The number of 
pairs emitted under the same conditions should 
be given by either of the positron plots. It should 
also be given by either of the areas C or D, since 
these are thought to be composed of negative 
members of pairs. 

It is desirable to introduce a refinement in the 
above computation, for the following reason. 


HALPERN AND H. R. CRANE 


There is no doubt that many of the electrons and 
positrons of very low energy are lost in the 
experiment, due to stopping and scattering in the 
material surrounding the source. The fact that 
the positron plots fall off rapidly at low energy 
is an indication that this is true. Therefore better 
estimates of the true areas under the curves are 
obtained by using the areas A, B, and the upper 
halves (2.25 to 4.5 Mev) of the pair plots, and 
calculating the total areas from these with the 
aid of the theoretical forms of the Klein-Nishina 
and pair curves, respectively. This treatment 


_ gives a ratio of 11 pairs per 1000 gamma-ray 


quanta, which is the pair internal conversion 
coefficient. 

Several theoretical calculations of the coef- 
ficient have been made,® and indicate that for 
the fluorine radiation the value should lie between 
1.2 and 2 pairs per 1000 gamma-rays, depending 
only slightly upon whether dipole or quadrupole 
radiation is assumed. Our value of 11 per 1000 
does not seem to be compatible with the theo- 
retical value, on the assumption that the angular 
momentum change is not greater than 2 units 
and that either Ne” or O'* emits the radiation. 


SEARCH FOR Low ENERGY GAMMA-RADIATION 
AND FOR DELAYED ACTIVITY 


The reaction giving rise to the gamma- 
radiation remains a mystery. A number of 
possibilities have been discussed,* but none 
seems to give a satisfactory energy balance. We 
have searched for gamma-radiation of low energy 
which might have escaped attention previously, 
and have found nothing down to 0.5 Mev. We 
have also tested for delayed emission of electrons 
or gamma-radiation by placing the CaF; inside 
the cloud chamber about 5 minutes after bom- 
bardment. In this case electrons of energy as low 
as 5 or 10 kev could have been detected. No 
activity was observed. 

The authors are indebted to the Rackham 
Endowment Fund for the support of this work. 


’ L. Nedelsky and J. R. Oppenheimer, Phys. Rev. 44, 948 


(1933); M. E. Rose and G. E. Uhlenbeck, Phys. Rev. 48, 
211 (1935); J. C. Jager and H. R. Hulme, Proc. Roy. Soc. 
148, 708 (1935). 

*A complete discussion has recently been given by 
Professor C. Breit, in a paper published by E. J. Bernet, 
R. G. Herb and D. B. Parkinson, Phys. Rev. 54, 398 
(1938). 
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An Independent Determination of the Binding Energy of the Deuteron 
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This paper is an account of an experiment designed to determine the binding energy of the 
deuteron by a method which is relatively insensitive to uncertainties in the energy-range 
relation for protons of low energy. The protons, produced by the disintegration of Th C” 
y-radiation, were observed in a low pressure cloud chamber in a strong magnetic field. The 
curvatures of the tracks allowed the calculation of the corresponding kinetic energies. The final 
value of the binding energy as given by this experiment is 


Vs =(2.17+0.05) X 10* electron volts, 


which is in reasonable agreement with previous determinations by others. The uncertainty 
0.05 Mev is the maximum uncertainty to be expected in this value of the binding energy. 





I. INTRODUCTION 


N 1935 Chadwick and Goldhaber' determined 

the binding energy Vz of the deuteron from 
the disintegration of deuterium by ThC” 
y-radiation of high energy. The particles taking 
part in this disintegration react according to the 
equation 


1H?-+-y—,H!'+ on 


Chadwick and Goldhaber detected the product 
protons with a proportional counter and linear 
amplifier, and estimated the kinetic energies of 
these protons by the amounts of ionization they 
produced in the counter. In this way they got 


Vg=(2.14+0.16) X 10° electron volts. 


Subsequently Chadwick, Feather, and Bret- 
scher*® performed with a cloud chamber a more 
precise experiment designed to measure the 
ranges of the protons which are liberated when H? 
disintegrates. After converting the observed pro- 
ton ranges to the corresponding proton kinetic 
energies by means of the best range-energy data 
available at that time, they got 


Ve=(2.25+0.05) X 10° ev. 


Their data, recalculated by Bethe’ in the light of 
more recent data on the energy-range relation 
for protons, yield, however, 


V»=(2.17+0.04) X 10° ev. 


1 J. Chadwick and M. Goldhaber, Proc. Roy. Soc. A151, 
479 (1935). 

?J. Chadwick, N. Feather and E. Bretscher, Proc. 
Roy. Soc. A163, 366 (1937). 

*H. A. Bethe, Phys. Rev. 53, 313 (1938). 


The discrepancy between these last two values 
of Vz is caused by previous uncertainties in the 
energy-range data for protons of low energy. 

At about the same time, J. R. Richardson and 
L. Emo‘ obtained 


Ve=(2.18+0.07) X 10° ev 


from the disintegration of H* by the high energy 
y-radiation from radioactive Na™. This value of 
Vz also involves the use of the energy-range rela- 
tion for protons, for Richardson and Emo de- 
termined the range of the product protons with a 
cloud chamber and then got the corresponding 
kinetic energy by the energy-range relation. They 
used almost the same relation as that which 
Bethe* used. 

Since these previous values of Vg depend upon 
the somewhat uncertain energy-range relation 
for protons of low energy (except for the first 
value quoted above, which is not of high pre- 
cision), it was thought desirable to determine Vz 
by a method which yields results that are rela- 
tively insensitive to this dependence. An attempt 
to do this with an ordinary magnetic spectro- 
graph with photographic recording was not 
successful, probably because the available Th C”’ 
y-ray source was not sufficiently strong. 

It was decided for the present experiment to 
use a deuterium-filled cloud chamber in which 
the total pressure was about 8 cm Hg and to 
place this cloud chamber in a strong magnetic 
field. Then the observable curvatures of the pro- 


‘J. R. Richardson and L. Emo, Phys. Rev. 53, 234 


(1938). 
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264 F. T. ROGERS, JR., 
ton tracks should enable the kinetic energies of 
the protons to be determined. 


II. BinpiInG ENERGY IN TERMS OF 
PROTON ENERGY 


In the following derivations we shall not con- 
sider either the effects of the relativistic varia- 
tions of mass or the small difference between the 
mass of the proton and the neutron. 

Let V be the energy in Mev of the Th C” 
y-ray which produces the disintegration of a 
deuteron ; Vg the binding energy of the deuteron 
in Mev; W the kinetic energy in Mev of the 
proton liberated by the disintegration of the 
deuteron ; p the ratio of the momentum of the in- 
cident y-ray to the momentum of the product 
proton (obviously p= 10*K'V/c(2W)'); K=e/m 
= specific charge of the proton, in e.m.u. per g; 
6=angle between the direction of the incident 
y-ray and the direction of emission of the product 
proton. 

Then if the equations for the conservation of 
energy and for the conservation of two per- 
pendicular components of momentum during 
the collision-disintegration process are properly 
solved, we have the simple expression 


Ve= V—2W(1—p cos 6+ ?’). (1) 


In this equation we treat V as a known quantity. 
We determine several pairs of values of @ and W 
(hence of p also) experimentally by making 
suitable observations. 


III. THEORETICAL CONSIDERATIONS FOR 
PROTON ENERGY 


If a product proton were moving in a vacuum 
in a uniform magnetic field,' its trajectory would 
be a right-circular helix, the axis of which would 
be parallel to the direction of the magnetic field. 
Actually, in the present experiment, the protons 
moved in an almost uniform magnetic field (to 
within 0.7 percent); but, moving in the cloud 
chamber, they did not traverse a vacuum. The 
loss of energy by a proton to the chamber gases 
thus produces a trajectory which is not at all 
like the right-circular helix mentioned above. 
Allowance must be made for this. 


+See, e.g., J. J. Thomson's Conduction of Electricity 
Through Gases (Cambridge Univ. Press, 1906), p. 106. 
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Let H be the strength in gauss of the (uniform) 
magnetic field in which the protons move; z a 
plane which is perpendicular to the direction of 
H. Let y be the angle between the direction of 
motion of a proton in H and a z-plane; y is of 
course measured in a plane perpendicular to a 
x-plane and will be constant (except for the 
effects of scattering by the chamber gases) over 
the length of a proton path. Let p be the radius 
of curvature of the z-plane projection of a proton 
path; p will of course decrease as the proton 
traverses its path. Let p’ be the value of p at the 
beginning of a proton path; and v the kinetic 
energy in Mev of a proton when its projected 
path has a radius of curvature p; v=W when 
p=p’. 

Then if for any point on a proton’s path we 
equate the centrifugal force acting on the proton 
to the force on it due to its motion in H, and if 
we solve this equation for v, we get 


v=5X10-"K (Hp sec yp)’. (2) 


Now if we multiply this by dv and integrate from 
v=0 to v= W, we get 


W =10-"K(He sec y)?, (3) 


where 
Ww 
=(1 w)| p*dv. (4) 
“0 


Eq. (3) gives W in terms of ¢, which is a function 
of the shape of the entire proton path, rather 
than in terms of a radius of curvature (e.g., as p’) 
at some one point on the path. The outstanding 
difficulty in the calculation of W from the ‘“‘one 
point’’ equation, 


W=5X10-"K(Hp’ sec y)?, 


is the determination of p’; in this experiment the 
accurate determination of p’ was not possible 
because the observed proton tracks were: (1) 
only slightly curved by the magnetic field; (2) 
often appreciably broadened by diffusion (charac- 
teristic of low pressure cloud chamber containing 
mostly deuterium gas) ; (3) occasionally scattered 
through appreciable angles. Eq. (3) allows these 
three obstacles to radius of curvature measure- 
ments to “average themselves out’’ over the 
length of a proton track. 
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Fic. 1. P;'P’P; is an actual proton path. PPP; is the 
t» projection of this path. a is measured from the Ps 
direction as shown. J;P;/; is the direction perpendicular to 
the & plane. The lengths of the several arcs are: P’P, =r’; 
P,'P’=L'—r'; PP,=r; P\P=L—r. 


IV. APPLICATION OF ENERGY-RANGE RELATION 


We now seek a more useful expression for ¢. 
To this end, we shall discuss the equation of a 
proton’s path in the cloud chamber, assuming K 
to be a constant over the whole path. Suppose we 
have a proton moving along its path. Since we 
wish to investigate ¢, which contains p, consider 
the projection of this path on a x-plane which 
passes through the end-point P: of the path 
(see Fig. 1). In this plane establish mutually 
perpendicular ¢ and » axes with origin at Ps, 
with the & axis tangent to the path’s projection 
at P». Let P; be the &y projection of the beginning 
point P,’ of the path, and let P be the £&y pro- 
jection of any point P’ of the path. To be associ- 
ated with P are the kinetic energy v at P’ and 
the radius of curvature p of the projected path 
at P. Let the range of the proton at P’ be 1’, 
and let the total length of the path be L’; 
then let 


r=r'cosy and L=L’'cosy (5) 
be the & projections of r’ and L’. Since the 
proton loses energy as it traverses the path 
P,'P:, the value of p decreases, a fact which we 
indicate by the functional relation 


p*= g(r). (6) 


The differential equation of the path in terms of 
the — and 7 coordinates is the trio of equations: 


dt=drcosa, dn=drsina, 
da=dr/[g(r) }, (7) 
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where @ is an angle (in the &y plane) defined 
adequately by Fig. 1. 

Let the chord P;MP,, M being the mid-point, 
be of length 2d. The perpendicular to this chord 
erected at M will intersect the projected path 
P:PP, at some point N;:let MN=h. Let the 
radius of the circle through P:, N, and P, be R. 
Then 

R= (h?+d*)/2h. (8) 


R, through the measurable 4 and d of a proton 
track, will serve as our measure of the radius of 
curvature of the proton track. We now seek a 
relation between R and @. 

For the relatively simple case in which p* is 
proportional to r’, we may take 


= Cy sec ¥, (9) 


where C, is a suitable constant. Substituting this 
into Eq. (7) and integrating, we get for the path 
of the proton 


§=}C, sec ¥(cos a+a sin hey (10) 
n=4C, sec ¥(sin a—a cos a), 


which are the equations of the evolute of a circle.* 
If from these equations R? is calculated by Eq. 
(8) for several paths of various lengths (i.e., for 
various values of a» (see Fig. 1)), and if o? is 
calculated for the same values of ao, and if then 
the ratios o?/R?® are formed, we get the values 
tabulated in Table I. Obviously, for proton 
tracks formed under the condition that pp’ is 
proportional to r’, 


o? = (1,230.03) R?, 
0.2=2(L’/C, sec ¥)'=1, 


independently of C,. (The relation between a» and 
L’ comes from Eqs. (7) and (9).) 

Now Eqs. (11) cannot be used unless Eq. (9) 
is known to be valid. Actually, the energy-range 
relation for protons of low energy is, within five 
or ten percent of energy, that given in Table II. 
These energy-range data are those of Livingston 
and Bethe’ as revised by them’ to include the 
recent experimental data of Parkinson, Herb, 
Bellamy, and Hudson.* It can be seen from this 

*See, eg., W. F. ood's Differential and Integral 
Calculus (Macmillan, 1927), p. 143. 

7M. S. Livingston and H. A. Bethe, Rev. Mod. Phys. 9, 
268 (1937). 


* D. B. Parkinson, R. G. Herb, J. C. Bellamy, and C. M. 
Hudson, Phys. Rev. 52, 75 (1937). 


(11) 








266 F. T. ROGERS, JR., 
table that the strictly linear energy-range rela- 
tion assumed in Eq. (9) is entirely consistent, for 
v£0.3, with the presently accepted relation. 
Therefore, we should not be in great error in 
taking Eq. (9) to be true. We must admit, how- 
ever, that p? may deviate somewhat from the 
linear relation of Eq. (9). 

We can see that even an appreciable deviation 
from Eq. (9) produces only an insignificant 
change in the value of o?/R*. To see this, we may 
work out the case in which (¥ = 0 for convenience) 


g(r) = Cyr/(aor®*+a,+a2r)’, (12) 
in which we shall take 


C= 120, 
a,=0.717, 


ao = 0.236, 
a2 = 0.0368. (13) 
This form for g(r) is quite similar to that sug- 
gested for protons of very low energy by the 
energy-range graph of Livingston and Bethe.’ 
Its average deviation from the function Cyr is, 
for proton energies up to 0.3 Mev, of the order 
of 10 percent. The constants in it are chosen to 
fit the following case: the stopping power of the 
gas is 0.075 that of air, and H=3520; the other 
condition on the constants is that Eqs. (12) and 
(13) represent (for this case) the data given in 
Table IT. 

If we substitute Eqs. (12) and (13) into Eqs. 
(7) and integrate, and if we then calculate several 


TABLE I. o*/ R? vs. ao. 











ao o?/R? ao @/R 
0.150 1.29 0.600 1.21 
0.200 1.26 0.750 1.22 
0.300 1.22 0.875 1.23 
0.450 1,21 1.000 1.24 








values of o*/R? as a function of L (instead of ao) 
as was done to obtain Table I, we find that 


o? = (1.27+0.05)R’, 

2.5£L<4. (14) 
Obviously the difference between Eqs. (11) on one 
hand and Eqs. (14) on the other hand, is essen- 
tially a small one relative to the difference 
between Eq. (9) and Eggs. (12)—(13). 
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V. EFFECT OF THE VARIATION OF PROTON’s 
CHARGE 


All our previous considerations, now, have been 
based on the assumption that the charge on a 
proton is a constant e=mK over the whole of 
the proton’s path. It is known, however, that the 
charge of a proton traveling through a gas may 
be neutralized for short intervals of time if the 
kinetic energy of the proton is small. As Living- 
ston and Bethe® point out, only if a proton has a 


TABLE II. Energy-range relation for low energy protons. 








RANGE OF PROTON IN CM Kinetic ENERGY 
In Dry Arr at 15°C OF THE PROTON 
AND 760 mu oF HG In Mev 





0.000 
0.202 
0.355 
0.463 
0.550 
0.628 


=Sesosss 
SCHARNS 











kinetic energy of 0.1 Mev or more, can one be 
reasonably safe in taking the charge to be the 
usual value, e= 4.8 X 10~"° e.s.u. Since the protons 
which we observed in the present experiment had 
energies less than 0.1 Mev over almost half their 
path lengths, we must make an allowance for 
their reduced charges at low energies if we intend 
to utilize their entire tracks for calculations. We 
shall derive a relation analogous to Eqs. (3) and 
(14) to allow for this variation of charge. 

Let E be the kinetic energy in ergs of the pro- 
ton when liberated by the disintegration of a 
deuteron ; 6 the ratio of the true charge on the 
proton to e (this, of course, varies as the proton 
traverses its path); ¢ the kinetic energy in ergs 
of a proton which has a range r’; i=E when 
r'=L’. Now just as Eq. (2) was obtained, we 
may also find 


i= (LTdep sec y)*/2m; 
therefore, 
v=5X10-“K (//dp sec y)?*. (15) 


Note that 
p= Cw/8?, (16) 


where C; is a suitable constant. Hence, just as 


*M. S. Livingston and H. A. Bethe, reference 7, p. 262. 
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Eq. (3) was obtained, 
W=10-"K(Hr sec y)’, (17) 


where 
Ww 
a= (1/W) f 5? p*dv. (18) 
0 


Note that for v2 0.18, 5=1 very nearly and this 
last reduces to 


r= hp’, 


To evaluate this, we need to know essentially 
how 6 varies with v and how p varies with v. 

We can get an idea of how 6 might be expected 
to vary with v from the data of Bartels’® on the 
variation of charge of positive rays of hydrogen 
passing through hydrogen. Bartels’ data, as 
given for gas pressures of 1 mm Hg, are sum- 
marized in Table III, in which m,/mo is the 
ratio of the number of protons to the number of 
neutral hydrogen atoms in a beam of such posi- 
tive rays. The third column in this table is the 
ratio n,/(n,+m9). These data are the circles in 
Fig. 2; the (dotted) curve for 4 in Fig. 2 is taken, 
by analogy with the variations of a-particle 
charges, to pass near the point (0.1, 0.9); its 
equation will be of the form é6=4(v). This par- 
ticular variation of 6 with v, while of course indic- 
ative of the nature of the true variation, cannot 
be expected to give much more than the order 
of magnitude of the true variation under the very 
different conditions which existed in the cloud 
chamber. Here the pressure was 8 cm (80 times 
as great as the pressure for which the data in 
Table III are quoted) ; the gas mixture was only 
? hydrogen (the remainder H?,O vapor), and the 
whole was in the presence of moderately intense 
y-radiation. We shall see shortly how this varia- 
tion of 6 with v must be replaced by another law 
for the present experiment. 

Continuing with the case (for simplicity) for 


v2 0.18. (19) 


TABLE III. Bartels’ data for n,/no as a function of v at 1 mm 








Hg pressure. 
v 0.010 0.035 0.060 
n/N 0.16 0.61 1.30 
6 0.14 0.38 0.57 








1° H. Bartels, Ann. d. Physik 13, 373 (1932). 
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Fic. 2. Values of 6 for various values of the kinetic energy v. 


which ¥=0 and for which v is proportional to r: 


where C,;=0.0713 (this is appropriate to the 
example for H=3520 and a stopping power of 
chamber gases of 0.075 that of air, which we have 
used in Section IV), we have 


p=([g(r) P= 4104/5(») ; (21) 


C:= 1687 from Eq. (16). Substituting this into 
Eqs. (7) and integrating, we get the path which a 
proton follows under these conditions. If from 
the path computed using the dotted-curve varia- 
tion of 6 with » (Fig. 2) we compute for several 
values of (r=)Z the values of R*, we find that 
for 25£L¢ 4.0, R=16.5+0.5. But the values of 
R which were found in the present experiment 
were in the neighborhood of 12 cm. Since the 
values of R calculated from Eq. (21) are largely 
independent of C3, this means that the true varia- 
tion of 6 with v is not the one given by Table III 
for 1 mm Hg pressure. 

The true variation of 6 with v must obviously 
be such that at very small values of », 6 is ap- 
preciably greater than is indicated in Table III. 
Furthermore, the true variation of 6 must be such 
that the values of R calculated from Eq. (21) by 
means of it, are somewhat smaller than 16.5 cm. 
An approximation to what must be the true 
variation of 6 with v is given by the solid curve 
of Fig. 2. With this law for 6(v) in Eq. (21), the 
calculated values of R in the neighborhood of 13 
cm, are more nearly equal to the observed values 
of about 12 cm. In view of the necessary experi- 
mental uncertainties in the observed values of R, 
this variation of 5(v) is probably sufficiently near 
to the true variation of 6(v) for the purposes of 
the present experiment. The values of r*/R?® got 
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by this 6(v) are given in Table IV. Obviously we 
can take r?=(1.11+0.09)R*® provided 2.5<L 
¢£4.0. Allowing for deviations from Eq. (20) as 
was done in Section IV, we may take 


7? = (1.15+0.09) R?, 
2.52 L¢ 4.0. 


Converting the bounds on LZ to bounds on W 
and substituting in Eq. (17), we thus get as our 
final expression for W, 


W =(1.15+0.09) X10-“K (HR sec y)?, 
0.175 < W<0.275. 


*(22) 


| as 


This expression for W is largely independent of 
uncertainties in the energy-range relation for 
protons of low energy, as are Eqs. (11) and (14). 


VI. Error THEORY 


Let ¢;, j=1 through 7, be the absolute errors 
which might be present in single determinations 
of R, ¥, W, V, 0, Ve, and in the mean value Vz 
of m determinations of Vz, respectively. Then it 
can be shown from Eq. (23), that 


€3 22W(Ak/2k+6:/R+e tan y), (24) 


where & is the coefficient 1.15. In a similar 
manner, 

€6 e,+5 K 10%; + es. (25) 
Hence, substituting Eq. (24) into Eq. (25) and 
noting that in the mean, some of the errors will be 
reduced by a factor of n! (according to the usual 
theory of errors), we have 


€7 Se, +2W(Ak/k) +[4Wei/R+ee tan v) 
+5X10'esJn; (26) 


here ¥ is the mean of the m values of y. 


TABLE IV. r?/R? vs. L. 








2.5 3 3.5 
1.02 1 1.14 


oo Om 
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VII. OBSERVATIONS AND CALCULATED 
RESULTS 


A low pressure (8 cm Hg) cloud chamber was 
used to observe the protons in this experiment. 
It was of a depth of about 2.5 cm and an inside 
diameter of 11 cm. A mixture of H*,O vapor and 
deuterium was used in the chamber, and a y-ray 


AND M. M. ROGERS 

source (of Ra Th) was placed as close as possible 
to the glass face of the chamber. The y-ray ac- 
tivity of the Ra Th was equivalent to about 5 
mg of Ra; this source was in a 4-mm thick lead 
container. In spite of the extreme propinguity of 
the y-ray source to the cloud chamber, no great 
trouble was experienced with the fogging of the 
chamber by the y-radiation. This was probably 
due to the low average atomic weight of the 
gas in the chamber. The light source was a carbon 
arc with photographic-flame carbons at 100 
amperes d.c. The photographs were taken with a 
stereoscopic camera having f:1.9 lenses. The 
magnetic field at the region of the cloud chamber 
was measured by the inductive method, with in- 
ductances, resistances, and standard cells which 
had been calibrated at the National Bureau of 
Standards. Under the conditions which prevailed 
when the photographs were made, H was found 
to be 3762 gauss to well within one percent. 

In all, about 2300 pictures were made, among 
which were found fifteen proton tracks which 
were suited to accurate measurement. The 
measurement of them was effected as follows: 

(1) the tracks were projected on a screen by 
running the (well illuminated) film back through 
the camera and moving the screen until the 
stereoscopic images coincided ; 

(2) with a Brown and Sharp protractor, ¥ and 
6 for each track was measured (see Table V)—the 
6 of a track was measured between the line from 
the center of the y-ray source to the beginning of 
the track and the tangent to the track at its 
beginning point; 

(3) the two ends and the mid-point N of each 
track were carefully recorded on small strips of 
paper with a very sharp pencil; this was done 
several times for each track; 

(4) the chord (of length 2w) P2P;’ (see Fig. 1) 
was measured and its half-length recorded (see 
Table V); 

(5) the mid-chord perpendicular (length /) 
from each track to its chord was measured and 
recorded (see Table V); 

(6) from w, d was calculated and recorded: 


d=w cos y. 


This method, far from being inaccurate, af- 
fords a very reliable method of measuring such 
curved proton tracks. In fact, it is believed that 
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TaBLe V. Data on proton tracks. 



































Ve ||\V¥e—Val 
Track | ¥( #® | w (cm) | & (cm) | d (cm) | (Mev) | (Mev) 
1 37 137 | 1.52 | 0.086} 1.215/| 2.220} 0.05 
2 30 173 | 1.25 | 0.062} 1.082] 2.197} 0.03 
3 17 63 | 1.59 | 0.096/ 1.520) 2.153| 0.02 
4 18 134 | 1.76 | 0.125 | 1.681 | 2.145] 0.02 
5 35 97 | 1.61 | 0.089) 1.319} 2.166} 9.00 
6 2 47 | 1.325 | 0.069) 1.350} 2.117} 0.05 
7 3 38 | 1.54 | 0.097| 1.540/ 2.198} 0.03 
8 42 159 | 1.33 | 0.064 | 0.975 | 2.261} 0.09 
9 0 76 | 1.45 | 0.086 | 1.450 | 2.167} 0.00 
10 27 131 1.30 | 0.062 | 1.160} 2.201} 0.03 
11 12 97 | 2.00 | 0.174/ 1.955|2.198| 0.03 
12 23 86 | 2.03 | 0.151 | 1.869] 2.201) 0.03 
13 9 O | 1.79 | 0.118] 1.769} 2.118) 0.05 
14 7 94 | 1.20 | 0.057} 1.191} 2.123} 0.05 
15 21 93 | 1.62 | 0,112} 1.512) 2.144] 0.03 
Mean 2.174} 0.03 














the determinations of w are accurate to within 
one percent and that those of 4 are accurate to 
within four percent. 

In the seventh column of Table V are the 
values of Vg calculated by Eqs. (23) and (1) from 
the data in the earlier columns of the table." In 
the eighth column are the deviations of the 
various values of Vs from the mean value Vz. 
The mean value of Vz, as seen on the table, is 
2.17 Mev. The average deviations of the several 
Vz's from this mean is 0.03 Mev, corresponding 
to a “probable error” of 0.02 Mev. 

The maximum error ¢; to be expected in 
this value of Vg we get from Eq. (26) taking 
€:,/R=0.05, 2 =0.025, «,=0.005, «6 =0.02, P= 19°, 
Ak=0.09, and m=15. Substituting these values 
and taking W =0.2 Mev, Eq. (26) gives e;=0.05 
Mev. Accordingly, we take as the final value 
given by this experiment 


V »=2.17+0.05 Mev. (27) 


It is to be noted that the uncertainty 0.05 is a 
conservative value indicating the maximum error 
to be expected in this value of Vg; as such, it may 
be expected to be of the order of two or three 
times as large as the “probable error” of Vz. The 
precision of this experiment thus compares well 


"With V=2.623+0.005 Mev. 





with the precisions of the previous determina- 
tions of the binding energy of the deuteron. 


VIII. Conclusions 


It will be noticed that the value V,=2.17 
Mev determined by this experiment is in excel- 
lent agreement with the value 2.17 found by 
Bethe’ on recalculation of the data of Chadwick, 
Feather, and Bretscher,* and with the value 2.18 
got.by Richardson and Emo‘ from the disintegra- 
tion of H* by the high energy y-radiation from 
radioactive sodium. This concurrence of values 
tends to indicate that the presently accepted 
energy-range relation (used by Bethe and by 
Richardson and Emo) is correct for protons of 
energies in the neighborhood of 0.2 Mev. It also 
tends to fix the mass of the neutron at the pres- 
ently accepted value of 1.0089 atomic mass units 
rather than at some higher value.” It also tends 
to indicate very definitely that the variation of 
protonic charge with energy as the proton 
traverses a gas such as the presently used 
H?,0 — H*, mixture, has 5(v) about twice as great 
(for »<0.1) as the experimental data for pres- 
sures of 1 mm Hg would predict—note that this 
conclusion depends chiefly upon the agreement of 
calculated values of R with the observed values 
of R, and is only confirmed by the agreement of 
V», with the other values of 2.17 and 2.18 Mev.* 


Finally, we wish to acknowledge our indebted- 
ness to Dr. H. A. Wilson for his having suggested 
to one of us that the binding energy of the 
deuteron be redetermined, and for his having 
provided us with every desired facility for carry- 
ing out the necessary experiment. 

We are also grateful to Miss Eby Nell Mc- 
Elrath for assisting in the taking of the necessary 
photographs, and to Mr. Stanley Heaps for help- 
ing with the development of the negatives. 


12 See, e.g., E. Hudspeth and T. W. Bonner, Phys. Rev. 
54, 308 (1938). 

* Note added in proof: In addition to the values of Vz 

uoted in the Introduction, G. Stetter and W. Jentschke 
Ceits. {. Physik 110, 214 (1938)) have more recently given 
the concordant value Vg =2.189+0.022 Mev. 
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New measurements of the structure of the green lines of Mn I combined with the data 
already available make possible a determination of 24 hyperfine interval factors. Interval factor 
formulas derived by the method of energy sums are in reasonable accord with much of the data. 
This permits an approximate determination of the nuclear magnetic moment of Mn®, which 


is found to be 3.0 nuclear magnetons. 





INTRODUCTION 


EASUREMENTS on the hyperfine struc- 

ture of the visible lines in the arc spectrum 
of manganese were first reported in 1909 by 
Janicki' who measured twenty-one structures 
using a Lummer Plate. These lines and others 
have since been remeasured by Wali-Moham- 
mad,? by White and Ritchl* and by Wali- 
Mohammad and Sharma.‘ White and Ritchl, 
whose measurements are the most extensive 
gave an analysis of the hyperfine structure 
showing the nuclear spin to be 2} units. They 
reported values for the total separations of 
several of the hypermultiplets and pointed out 
that the prominent splittings are associated with 
configurations in which an unbalanced s electron 
occurs. Using the vector model of the atom they 
developed formulas for the separations expressing 
the intervals in terms of coupling factors related 
to the individual electrons of the configuration. 
Following White and Ritchl’s work Wali- 
Mohammad and Sharma measured many of the 
lines used by White and Ritchl in their analysis 
and reported disagreements in some of the 
measured structures. We have remeasured a 
number of the lines in the green region concerning 
whose structure the existing data seemed least 
conclusive. An analysis of these data along with 
those of the previous observers yields values for 
the hyperfine intervals differing in some details 
from those obtained by White and Ritchl and 
permits calculations leading to an approximate 
value for the nuclear magnetic moment. 


1L. Janicki, Ann. d. Physik 29, 833 (1909). 

* Wali-Mohammad, Astrophys. J. 39, 185 (1914). 

+H. E. White and R. Ritchl, Phys. Rev. 35, 1146 (1930). 

*Wali-Mohammad and P. N. Sharma, Phil. Mag. 18, 
1144 (1934). 


EXPERIMENTAL DATA 


Table I gives the results of our measurements 
upon thirteen lines in the green region together 
with the data of the previous observers. The 
interval in wave numbers separating each 
measured component from the strongest line in 
its pattern is always given. The data of the other 
observers have been converted to the same units. 
Since many of the structures are only partially 
resolved, the measured positions may refer only 
to points of maximum intensity in the pattern. 
Components and maxima are designated by 
numbers specifying their relative intensities in 
the pattern. When a structure terminates in an 
unresolved “‘tail’’ the terminus is designated by 
an é in the intensity column. White and Ritchl 
also have indicated the extent of unresolved 
patterns. 

For our observations the manganese was 
excited in helium in a water-cooled metal hollow 
cathode discharge tube. Resolution was by 
means of a ruled Fabry-Perot interferometer 
and small Hilger glass prism spectrograph. 
Interferometer spacers of 2.37, 5.28 and 10.03 
mm were employed. Three or more spectrograms 
were taken with each spacer so that each pattern 
reported in Table I was verified by its appearance 
under different resolutions and in different 
intensities, with the exception of the lines 
4A5399.5 and 5413.7 which were measurable on 
only two plates because of their low intensities. 

It may be seen (Table I) that in most cases in 
which the new measurements differ materially 
from the older the new observations reveal 
additional components not before resolved. This 


*R. C. Machler and R. A. Fisher, J. Opt. Soc. Am. 25, 
315 (1935). 
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is illustrated by the first line of the table, nection with \5341.4 for which White and Ritch] 
\5537.8, in which the pair of components at reported the same pattern as we, but upon a 
0.129 and 0.221 cm-, respectively, were observed scale some 10 percent larger while Wali-Moham- 
as a single component in an intermediate position mad and Sharma found a similar pattern on a 
by Wali-Mohammad and by Wali-Mohammad -—————— = nN ain Janicki seen 
° . interferograms. It is evident from Janicki’s pu 
and Sharma. The same is true of the pair at 0.320 spectrograms as well as from his measured separations 
and 0.377 cm~'. White and Ritchl did not resolve oa oe wr yey re lines —_, by or as 
: a * ses 7, 5399.5 a A, respectively, are in reality 
structure in this line and so gave only the initial \A5399.5, $407.4 and $4204. respectively. Wali-Moham- 
and terminal points of the pattern. The most mad apparently followed Janicki’s example in this identi- 
: : : . : _ fication. We have taken the liberty of shifting these 
— disagreement in data* appears in con measurements by both Janicki and Wali-Mohammad in 
‘The close grouping of lines in the vicinity of 5400A the table, thus bringing them into close agreement with 
has led to some confusion in their identification in the those of the later observers. 
TABLE I. Hyperfine structure of green lines of Mn | according to different observers.* 
ents . 
ther 4r(cm™~) 4r(cm™~') 
The MA) Ist New M&S WER ™M ) || »a) |inr New M&S WER M cz 
>ach | 6 0.000 0.000 0.000 0.000 1 —0.10 
oa | 3 .129 6 000 0.000 
.154 .150 2 .144 
ther 5537.8 | 3 .221 5413.7 171 
. |} 2 .320 2 .237 
uts. | 346 340 e 275 
ally | 1 376 366 
—— —__—_—_——_— 5 .000 =—0.,000 .000 §=60.000 30.000 
only | 6 .000 000 .000 000 0,000 5 ‘182.190 190 «185,189 
ern. | 1 .061 5407.4 3 345 341 .357 .350 
b 1 .150 : 2 486 ~— 506 505 490 
y 5516.8 | 164 556 
: in 7 3.220 e 623 
247 246 = 230.240 ||— — 
an 1 .264 1 —.114 
by 2 372 410395388 || s399.5 | 4 000 000.000 
chl 5 000 000 = .000S——=.000 3 .230 .234 = .223 
ved 3.149 AS20184 154 -— ~ — 
c 3 .270 x j d 
a. 295 290 293 4 — 14 = 112 
vas _ 349 4 —.21 —.210 
al | ¢ 408 392 ST 1 3 = 308 —|289 
Ow —_ > _______—_ ————$__—_—__—_— 2 — 383 — .357 
by 4 000 000 000 §©.000 ~§6.000 1 —.410 — 395 
5 3 .207 .206 221 .206 217 || - 
er 5481.4 2 373 373 391 383 407 | 1 — .098 
oh. |} 1.463 6 000 000 ©=—-.000 
03 le 516 512 || 5377.6 | 3 101 
. LE ee st — 121 
ms | 4 000 000 = 000.600 || 2 130 130 
5470.7 | 3-191 182 184 187 
rn , 3 350 327 347 351 6 .000 .000 .000 .000 .000 
—_——_——_|—_—— ee — || 4 374 .309 406 376 .379 
nt | 6 000 .000 5341.1 | 3 514 «6.428 «=—.560 Ss S23—s«523 
val | 5 —-106 —.105 3 620 .530 680 
$432.6 | 4 —.195 —.186 667 642 
on | 3 —.266 — .254 2 .680 .762 
2s, e —.310 — .298 | 
in | 5 000 000 ©=.000 000 | 
ly 4 190 197 185 195 
3 341 .356 .356 .359 
al $420.4 2 440 
: 475 I} 
a 1 510 503.492 || 
e .59 .627 | 
5, 
* M&S indicates data of Wali-Mohammad and Sharma;,W&R_indicates data of White and Ritchl; M indicates data of Wali- Mohammad : 
J indicates data of Janicki 
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TABLE II. Hyperfine structure of violet lines of Mn I from 
measurements of White and Ritchl.* 












































(A) Ar(cm~") (A) 4v(cm~) 
0.0000 0.0000 
2147 4048.8 1683 
4018.2 
2946 2251 
"3565 
4029e 0000 
0742 
0000 1605 
—|0930 4055.5 2207 
~"1724 2633 
4030.7 ~'2352 2925¢ 
—'2826 
— "3138 0000 
0619 
0000 4063.5 1426 
~'0873 2209 
4033.1 — "1654 iin 
~'2207 0000 
~'2576 0834 
4079.2 1494 
0000 1986 
— ‘0891 26166 
4034.5 —'1610 
—'2083 0000 
0000 _ 
0829 0000 
4035.7 1805 0672 
2462 4082.9 1218 
3113¢ 1638 
0000 0000 
1010 4083.6 0654 
1867 19560 
4041.4 roo 
"3146 
"3513 








* White and Ritchl gave intervals between successive components 
in cm™~!, These intervals have been added to give positions with respect 
to the st st component of the pattern. ¢ following a position indi- 
cates the end of the pattern rather than a component. 





scale 17 percent smaller. The new measurements 
on this line are, however, in close agreement with 
those of Janicki and of Wali-Mohammad. 
Other data of use in this analysis are the 
excellent measurements of White and Ritchl 
upon the group of lines in the vicinity of 4000A 
listed in Table II. The measurements of other 
observers in this region are too fragmentary for 
inclusion in the table. White and Ritch! have 
also given microphotometer curves depicting the 
complex unresolved patterns found for \\4754.1, 
4783.4 and 4823.5 with a wave-length scale. 


ANALYSIS OF HYPERFINE STRUCTURE 


The term classification of the Mn I lines’ here 
of interest is shown by Table III. The wave- 


™See R. Bacher and S. Goudsmit, Atomic Energy States 
(McGraw-Hill Book Company), for references. 


lengths of the lines appear within the rectangles 
and the terms from which they originate along 
the margins ; even terms are above and odd terms 
to the left. 

All of the information now available is in 
complete accord with White and Ritchl’s assigned 
value for the nuclear spin, J=2}. Accepting the 
value of J as established and assuming that the 
interval rule for hypermultiplets is obeyed 
within the limits of observation, we may proceed 
to apply the graphical method of hyperfine 
structure analysis described by Fisher and 
Goudsmit.* The scheme employed throughout 
the analysis of the twenty-rfine hyperfine 
structures has been to construct for each pattern 
what may be called a “pattern graph” repre- 
senting the observed positions of components on 
a varying scale by means of properly spaced 
diverging pencil lines. This pattern graph which 
is On semi-transparent paper is then laid over the 
appropriate analysis graph and shifted about 
until the position is found at which the corre- 
sponding lines of the two graphs coincide. This 
determines the ratio between the interval factors 
of the two hypermultiplets. The relation of the 
two graphs for \5537.8 is illustrated in Fig. 1. 


TABLE III. Classification of Mn I lines. 
































3d* 4s *D 
J 44 34 24 1} j 
1} 5481.4 5516.8 5537.8 
3d5 4s(5S)4p *P® 24 5420.4 5470.7 5505.9 
34 | 5341.1 5407.1 
44 | 4041.4 4079.2 
34 | 4018.2 4055.5 4083.6 
3d* 4p *D® 24 4035.7 4063.5 4082.9 
1} 4048.8 4079.4 
J 3d5 452 *§Sxy 3d* 4s 5s *Sy 
24 5432.6 4754.1 
3d 4s4p*P® = 33} 5394.7 4783.4 
4} 4823.5 
1} 4034.5 
3d* 4s(7S)4p *P® 24 4033.1 
3} 4030.7 
J 3d* 4s Ss *Sy 
24 5377.6 ; 
3d* 4s 4p *P® 1} 5399.5 
j $413.7 











*R. A. Fisher and S. Goudsmit, Phys. Rev. 37, 1057 


(1931). 
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The indicated ratio of interval factors is 
A(®D,)/A(®P;,;°) = —0.87. From the measured 
spacings we then get A(*D,)=0.028 cm™ and 
A(®P,;°) = —0.032 cm“. 

The pattern of \5537.8 is the most completely 
resolved of the observed structures, since five 
of the expected six components are measurable. 
Using this line as key we have carried through 
the analysis and obtained values for the interval 
factors of all states involved. These interval 
factors appear in Table IV. Since most of the 
patterns are less completely resolved than 45537, 
an unambiguous determination of the interval 
factors cannot always be made from a single 
structure. However, the requirement of con- 
sistency in the interpretation of related patterns 
provides in every case more conditions than are 
necessary for a unique assignment of values to 
the interval factors. It is to be recognized that 
the interferograms contain significant informa- 
tion having to do with the contours of the 
patterns which cannot conveniently be tabulated, 
but which nevertheless can be used to advantage 
in the graphical method of analysis when the 


TABLE IV. Hyperfine interval factors, Mn 1. 








INTERVAL Factor (cw) 








3d 432 6S Se 0.000 k 








} 0.0278 
} 0165 | 
3d* 4s *D } 0145 } +0.0005, A 
} 0135 | 

} 0145 

} .0222 

} 0205 >+ .001, A 
} .018 





Sd*4s4p"P* 





} — .056 
5 —.023 »>+ .001, A 
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3d* 4s 4p *P® 
—.020 } 





.0230 
0185 > + .001, A 


1} 
3d* 45(7S)4p *P® 2) 
34 0160 | 





— .0315 | 
—.0215>+ .001, A 


14 
3d54s(5S)4p *P® | 24 
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*S. Goudsmit, Phys. Rev. 37, 663 (1931). 
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Fic. 1. Graphical analysis of \5537.8 of Mn I. Above: 
Analysis graph showing theoretical line distribution for 
transition J = }->J=1} with varying interval factor ratio, 
A(4)/A(14). Theoretical relative intensity of each compo- 
nent is indicated. Below: Pattern graph showing observed 
pattern on varying scale. Estimated relative intensity of 
each component is indicated. Coincidence between the 
two graphs occurs at A(})/A(1}) = —0.87. 


plates are available for study. Aside from the 
number of significant figures with which they are 
expressed the interval factor data cannot all be 
given equal weights because of the varying 
element of judgment which is represented in the 
interpretation of the patterns. We have indicated 
the degree of our confidence in the data of 
Table IV by grading them A, B and C, A indi- 
cating that we regard the figure as comparatively 
reliable, C that we regard it as somewhat dubious. 


NUCLEAR MAGNETIC MOMENT 


In order that the nuclear magnetic moment 
may be found from hyperfine structure data it is 
necessary that formulas expressing the observed 
interval factors in terms of the coupling between 
the nucleus and the individual electrons be 
available. Goudsmit* has shown how such for- 
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mulas may be obtained for configurations through 
application of the method of energy sums. Breit 
and Wills’ in extending the theory have shown 
that the electron ‘‘coupling constants” used by 
Goudsmit are in general not constant for all 
states of a configuration but depend upon the 
inner quantum numbers of the states. They 
derived improved interval factor formulas for a 
few two- and three-electron configurations and 
Crawford" and Crawford and Wills” have 
applied the same method to the derivation of 
formulas for the d*s and p*s configurations. Com- 
putations by the method of Breit and Wills have 
so far not proved feasible for the important 
configurations of Mn I, namely d'‘s, d‘s p, d's-s 
and d*p. We then have no choice but to rely 
upon the admittedly approximate method of 
Goudsmit. There is some justification for con- 
fidence in this method as applied to Mn I in view 
of the fact that the hyperfine structure may 
largely be attributed to the 4s electron and that 
Mn I isa celebrated example of Russell-Saunders 
coupling. 

The method of energy sums depends upon the 
assumption that an externally applied magnetic 
field may remove some couplings completely and 
leave others entirely undisturbed. Thus explicit 
expressions may be written for the various inter- 
actions in each stage of decoupling. It is assumed 
that the interactions between the individual 
electrons and the nucleus may be expressed by 
coupling factors associated with individual 
electrons. These factors are assumed to be 
constant within a given configuration but may 
vary from configuration to configuration. With- 
out giving in detail the rather laborious algebraic 
steps leading to the interval factor formulas we 
will mention briefly the stages of decoupling 
which it has been necessary to consider in order 
to obtain each set of formulas appearing in 
Table V. 

3d*(*D)*s *D. Here it is necessary to consider 
two stages of decoupling. The first is that in 
which merely the coupling between the nuclear 
magnetism and the total electron configuration 
is broken, the configuration itself being coupled 


1G, Breit and L. A. Wills, Phys. Rev. 44, 470 (1933). 
" M. F. Crawford, Phys. Rev. 47, 768 (1935). 
985); F. Crawford and L. A. Wills, Phys. Rev. 48, 69 
1 ’ 
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as before. The second stage assumed is that in 
which the 4s electron is broken away from the 
configuration leaving the strongly coupled 
3d°(*D) group acting as an entity having the 
quantum numbers characteristic of a °D state. 

3d*(*S)4s 4p *P®. Here also two stages must be 
considered, i.e., (1) the nuclear spin decoupled 
from the configuration which retains the proper- 
ties of an *P, and (2) both the 4s and 4 electrons 
decoupled from the 3d5(*S) group which retains 
the properties of a °S. 

3d*(*S)4s 4p *P*(a) and *P°(8). Two *P* states 
arise in this configuration. The symbols a and 8 
will be used to designate the *P states of lower 
and of higher absolute energy respectively. One, 


TABLE V. Interval factor formulas derived by method of 
energy sums. 








d*(*D)s *D 


A(4) =(7/15)a(s) — (4/3)a(d*) — (4/3) b(d*) 
A(14) =(13/75)a(s) +(2/15)a(d*) + (74/105)b(d*) 
A(24) =(23/175)a(s) +(12/35)a(d*) + (164/245) d(d*) 
A(34) =(37/315)a(s) +(26/63)a(d*) + (122/441)d(d*) 
A(4}) =(1/9)a(s) +(4/9)a(d*) — (20/63)b(d*) 
rA =a(s) 
d5(*S)sp spe 
A(24) =(9/49)a(s) — (88/245)a(p) 
A(3}) =(59/441)a(s) + (80/441 )a(p) 


A(44) =(1/9)a(s)+(8/45)a(p) 
LA =(3/7)a(s) 


d(8S)sp «Pe 
A(4) = — }a(s) — (8/15) a(p) 
A(14) = —(11/75)a(s) +(32/375)a(p) 
A(24) = —(3/25)a(s)+(56/125)a(p) 
LA = — ja(s) 


d*(*S)sp a *P® and B *P® 


aA(14)+8A(1$) =(2/25)a(s) — (104/125) a(p) 
aA (24)+8A (24) = (62/1225)a(s) + (1696/6125)a(p) 
aA(34)+ 8A (34) =(2/49)a(s) + (136/245) a(p) 
(aA +BA) =(6/35)a(s) 
a tpe 
A(14) =(4/15)a(s) —(8/25)a(p) 
A(24) =(124/735)a(s) —(8/1225)a(p) 
A(34) = (20/147)a(s) + (16/49)a(p) 
LA =(4/7)a(s) 
pepe 
A(14) = —(14/75)a(s) — (44/75)a(p) 
A(24) = —(62/525)a(s) + (248/875)a(p) 


A(34) = —(2/21)a(s)+(8/35)a(p) 
LA = — fa(s) 


d*(*D)p spe 


A(4) = —(8/225)a(p) — (10/9)a(d*) — (8/15) b(d*) 
A(14) = —(86/1125)a(p) +(1/9)a(d*) +(148/525)d(d*) 
A(24) = —(32/875)a(p) + (2/7)a(d*) + (327/1225)d(d*) 
A(34) = (142/4725)a(p) +(65/189)a(d*) + (244/2205)b(d*) 
A(44) =(16/135)a(p) + (10/27) a(d*) — (8/63) d(d*) 
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the lower in energy, may be described as d5(*S)s 
(7S)p *P°, the other as d*(*S)s (*S)p *P®. While 
these designations are useful as a means of 
distinguishing the states they cannot be regarded 
as an accurate statement of the manner in which 
the states are built up since each *P® must be 
thought of as built in part upon the “7S and in 
part upon the °S of the ion by addition of the p 
electron. If we are to obtain interval factor 
formulas for these states and avoid this am- 
biguity we may consider only the sums of the 
interval factors of states having the same J value 
in the two multiplets. These formulas for the 
sums of interval factors may be obtained by 
considering the same stages of decoupling as in 
treating the *P°. 

In Mn I the two states *P°(a) and *P°(8) are 
found to be some 11,000 cm™ apart. This sug- 
gests that the designation mentioned above may 
be fairly accurate. Proceeding on this assumption 
we may obtain interval factor formulas for the 
two multiplets separately. For this purpose it is 
necessary to consider an intermediate stage in 
which the 4p electron is decoupled from the 
3d§ 4s (7S) or (5S) group as the case may be. 
While this is admittedly a questionable procedure 
it leads to relations which seem to have some 
significance when compared with the experi- 
mental data. A variation calculation too involved 
for inclusion here indicates that values for the 
coupling factors obtained by introduction of the 
experimental interval factors into formulas for 
a®P® and 8*P* will represent lower and upper 
limits, respectively, for the ‘‘correct’’ value. 

4d‘(°S)4s 4p *P®. The same considerations as 
applied to the *P® of this configuration lead to 
the formulas. 

3d*(°D)4p *D®. Since several states arise from 
this configuration the formulas can be obtained 
only by making a further assumption not used 
in deriving the other formulas, i.e., that the sums 
of energies of all states characterized by the 
same values of M, and Ms are independent of 
the degree of decoupling. 

In the formulas of Table V thie interval factors 
of the hypermultiplets are designated by A 
followed by the inner quantum number of the 
particular state in parentheses. Each electron 
coupling factor is designated by a followed by 
parentheses indicating the particular electron or 
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electron group to which it refers. The interaction 
of the d* group is expressed by two factors, one 
relating to the spin and the other to the orbital 
interaction. The orbital factor is represented by 
a(d‘) and the spin factor by b(d*) in the formulas. 

Introducing the experimental interval factors 
into these formulas we obtain numerical values 
for the coupling factors. In most cases more 
relations exist than there. are constants to be 
determined so that consistency provides a test 
of the accuracy of the formulas. A convenient 
means of determining the coupling constant for 
the s electron is provided by the common 
property of all formulas of Table V, that the 
sum of the interval factors of all states within a 
given multiplet depends only upon the contri- 
butions of unpaired s electrons. The values given 
below are those found to give best consistency 
within a given set of formulas. These results are 
graded A, B and C in accord with their relative 
reliability. The coupling constants determined 
for the different states aré as follows: 


3d* 45 *D 
a(4s) = 0.0863 cm= 
a(3d*)=0.0094 “ }A 


b(3d*) <0.0001 “ 
3d® 4s 4p *P* 

a(4s) =0.143 cm=", A 

a(4p)=0.009 “ ,B 
3d* 4s 4p *P? 

a(4s) =0.162 cm“, A 

a(4p)=0.006 “ ,C 


3d° 4s 4p a *P® 
a(4s)=0.095 cm, B 
a(4p), small, values inconsistent. 
6 *P° 
a(4s)=0.175 cm=", B 
a(4p), small, values inconsistent. 


As an example of the relative consistency 
within a set of formulas we may consider the 
3d*°4s*D terms, for which the experimental 
interval factors are rather accurately known. 
Substituting the values, a(s)=0.0863 cm™, 
a(d*)=0.0094 cm, b(d*)=0, into the five 
formulas, we obtain for the interval factors 
A(4) to A(44), respectively, 0.0278, 0.0162, 
0.0145, 0.0140 and 0.0138 cm~'. These are to 
be compared with the experimental values of 
Table IV. 
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The only coupling factor that adapts itself to 
a calculation of the nuclear magnetic moment is 
that of the 4s electron. In addition to the fact 
that our formulas are probably less reliable in 
reference to a p electron, the smallness of the 
indicated values for a(4p) makes them extremely 
sensitive to errors in measurement. The value of 
a(4s) obtained from the 3d* 4s *D state is 
probably the most reliable. The values for the 
same factor determined from the 3d° 4s 4p *P° 
and ‘P® are considerably larger and in only 
moderate agreement with one another. That they 
should be larger is expected, since the screening 
of a 4p electron must be rather less than that of 
a 3d electron. The interval factor formulas for 
the two *P* states are the least satisfactory when 
applied to the observed interval factors. How- 
ever, if the values of a(4s) obtained from these 
two multiplets are regarded merely as an upper 
and a lower limit for the correct value as men- 
tioned above, then these values are in accord 
with those obtained from the *P® and *P® of the 
same configuration. The observed interval factors 
from the 3d° 4s 5s 8S and ‘S terms are not useful 
because formulas derived by the method of 
energy sums do not distinguish between the con- 
tributions of the two s electrons. 

Goudsmit’s relation” expressing the nuclear 
g factor in terms of the coupling between an s 
electron and the nucleus may be written in the 


form, 
R'Z,1838 


a?Z,WIK (4 Zi) 





g(I) = ga(s) 


Here W represents the energy in cm™ necessary 
18S. Goudsmit, Phys. Rev. 43, 636 (1933). 
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to remove the s electron from the atom in that 
stage of ionization in which it is the only outer 
electron. The other symbols have their conven- 
tional meanings. The computation must be made 
for each configuration since a(s), Zo and W are 
different for different configurations. 

Considering the state 3d*° 4s *D we have, 
a(4s)=0.0863 cm“, W=53,500 cm, Z9=1, 
Z;=25, K=1.06. W is the approximate energy 
necessary to remove the 4s electron from an 
atom in the above state leaving the ion in the 
state 3d* *D, the centers of gravity of both states 
being taken into account.’ Putting these numer- 
ical values into the formula we obtain 


g(J) = 1.13. 


Considering the configuration 3d° 4s 4p we 
have two indicated values of a(4s), i.e., 0.143 
and 0.162 cm~. Since the larger one is to be 
favored in such cases" we suppose that the value 
of 0.160 cm is satisfactory for use in computa- 
tion. Here W is the energy needed to remove the 
4s electron from the configuration 3d5 4s of Mn IT 
leaving the configuration 3d° of Mn III. If 
centers of gravity are taken into account, 
W =123,000 cm, approximately. Also Z»=2. 
We now obtain from the formula, 


g(Z) = 1.20. 


In view of the comparative agreement between 
these two independent calculations we believe 
that the value 1.2 may be assigned as the nuclear 
g factor of Mn® until more precise means for its 
determination are available. Multiplying by the 
value of J, we obtain a nuclear magnetic moment 
of 3.0 nuclear magnetons. 
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In an attempt to explain the discrepancy between Mott's 
theory and the observational results on the polarization of 
electrons in the double scattering experiment, we have in- 
vestigated the depolarizing effect of (1) multiple elastic 
scattering, (2) inelastic scattering with spin change of the 
incident electron and (3) exchange scattering in which the 
exchanged electrons have opposite spins. Although mul- 
tiple scattering actually takes place in the Au foils used in 
the experiments, the depolarization caused by such scatter- 
ing is negligible. For a screened Coulomb field and electrons 
of 100 kv the depolarization in a 2.5X10~* cm Au foil is 
only 2.2 percent whereas a 90-100 percent depolarization 
is needed to reconcile theory and experiment. The de- 
polarization caused by inelastic scattering is about 100 
times smaller than that caused by elastic scattering. The 
small depolarization in these two cases is due to the fact 
that while there are many elastic and inelastic collisions, 
most of these collisions take place with small scattering 
angle 3 and for small angles the depolarization per collision 


is small (of order 3). The main contribution to the de- 
polarization by inelastic collisions comes then from large 
scattering angles at which inelastic scattering is relatively 
infrequent. The exchange scattering contributes practically 
no depolarization (~10~"). While the depolarization per 
collision in each exchange process is relatively complete, 
the number of exchange processes is extremely small due to 
interference (~10~™ times the number of elastic collisions). 
Therefore the absence of polarization in electron scattering 
cannot be explained on the basis of the processed con- 
sidered here. 

A discussion of the validity of the first Born approxima- 
tion is given (Appendix §2). It is shown that in the case 
that the Born expansion parameter ¢*Z/Av is of order 
unity, the first Born approximation gives a result for the 
scattered intensity which is accurate to within 20 percent 
or better. The depolarization is given with an accuracy of 
18 percent or better. 





I. INTRODUCTION 


NE of the important consequences of the 
Dirac theory of the electron is the pre- 
diction that an initially unpolarized beam of 
electrons will become partially polarized as the 
result of an elastic scattering process. If this 
partially polarized beam of electrons is made to 
undergo a second scattering, the intensity of the 
twice scattered beam, according to the theory, 
should be asymmetrical in the azimuth about the 
direction of incidence on the second target. 
In fact, the calculated scattering consists of two 
parts: The first, independent of the azimuth, is 
the product of the relativistic Rutherford intensi- 
ties corresponding to the scattering angles at the 
two targets. The second part arises from spin- 
orbit interaction in the scattering field; it de- 
pends on the orientation of the magnetic moment 
with respect to the second scattering plane, that 
is, it is azimuth dependent. Obviously this term 
yields a nonvanishing contribution to the scatter- 
ing only after the spin axes have been given a 
preferential direction. 
The quantitative theory of the polarization of 
electrons by double scattering has been given by 


Mott.' According to his exact treatment of the 
scattering in a Coulomb field the expected 
azimuthal asymmetry is large enough to be 
observed experimentally if the following condi- 
tions are fulfilled: (1) only single scattering 
should take place in the two targets,* (2) the 
atomic number of the scattering atoms of both 
targets should be large (Z/137 comparable with 
unity), (3) both scattering angles should be 
large (of order +/2), and (4) the energy of the 
incident electrons should not be too small.’ 

The intensity of the doubly scattered beam is 
proportional to J where 


J=1+46(8;, 32) cos ¢. (1) 


Here J, and #2 are the two scattering angles, 
¢ is the azimuthal angle and 4 is a function of 
the two scattering angles which has been calcu- 
lated by Mott. The asymmetry, that is, the 
quantity given in the experiments is 


26=(Jo—Jz)/JIm, (2) 


1N. F. Mott, Proc. Roy. Soc. 135, 429 (1932). 

2As will be seen in the following, this condition is 
perhaps too stringent. 

*For Au the theoretical asymmetry is negligible below 
30 kv, attains a,(rather flat) maximum at 130 kv and then 
very slowly approaches zero for infinite energy. 
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where the subscripts refer to the azimuthal angle. 
For 38;=0,=7/2, the scattering in Au of elec- 
trons of 90 to 160 kv should be characterized by 
an asymmetry of 12 to 16 percent. 

Very careful experiments for the purpose of 
observing the asymmetry have been carried out 
by Dymond.‘ With electrons of 90, 125 and 
160 kv energy scattered through angles 7/2 at 
two Au foil targets the asymmetry observed was 
of the same magnitude as the experimental 
error of one percent. That is, if there is any 
asymmetry ai all, it is small enough to be com- 
pletely masked by the one percent experimental 
error due to the small deviations from ideal 
geometry, etc. 

In the attempt to reconcile the theoretical 
results of Mott and Dymond’s experiments we 
may first inquire whether all the conditions 
mentioned above were fulfilled. The only con- 
dition the fulfillment of which may be open to 
question is that of single scattering. It can be 
shown (Appendix §3) that an appreciable amount 
of multiple scattering (about 50 collisions) takes 
place even in the thin foils used by Dymond 
(2.5 10-* to 5 10~* cm thickness). It is evident 
that multiple scattering will result in some 
depolarization of the beam but just how much 
remains to be seen. 

In addition to multiple scattering we have 
investigated the depolarizing effect of other 
processes. These are : (1) inelastic scattering with 
change of spin of the incident electron and 
(2) exchange scattering in which the exchanged 
electrons are of opposite spin. 

The result of these considerations may be 
stated very briefly. Unfortunately, none of the 
effects considered produces any appreciable de- 
polarization of the electrons and the discrepancy 
between theory and experiment remains—per- 
haps more glaring than before. 


II. THe THEORY OF POLARIZATION 


Before discussing the depolarizing effect of 
multiple scattering it will be most convenient 
to recapitulate some of the general theory of 
double scattering.' We consider an electron with 


*E. G. Dymond, Proc. Roy. Soc. 136, 638 (1932), 145, 
657 (1934). For further references and a critical survey of 
the experimental work see H. Richter, Ann. d. Physik 28, 
533 (1937). 
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wave number & incident along the z axis and 
scattered elastically in the direction 3, ¢ by 
some field of force. The wave function describing 
the electron is asymptotically 


v~se"+u,(d, gle*"/r, v=1---4, (3) 


where s, are the components of the amplitude of 
the Dirac plane wave. In particular we have 


$3=A, ss=B, (4) 


where A and B are constants related to the 
direction of the spin axis (xo, wo) in the incident 


beam. 
—B/A=cot (xo0/2)e‘*®*. (5) 


We may further impose the normalizing con- 
dition 
|A|?+|B\/?=1. (6) 


The differential cross section for scattering into 
the solid angle dQ is then given by 


Hod 2 = (| w3|*+ | ug|*)dQ. (7) 


The general solution of the scattering problem 
may be obtained by a superposition of the two 
waves corresponding to A=/, B=0 and A=0, 
B=1. Mott has shown that in these two cases 
the asymptotic forms of the wave functions yp; 
and y, are given by: 


A=1, B=0: 
vsrwe™+f(dje"'/r 
7 (8a) 
vir~wg(d)e've*’/r, 
A=0, B=1!: 
y ~— —g(d)e—*ve'*r ry 
; (8b) 


ysrwe*+ fide”, Tr, 


where f and g are independent of ¢ and are the 
same functions in (8a) as in (80). By superposi- 
tion we find for the scattered beam 


u3=Af—Bge'*, us=Bf+Age'*. (9) 


For the scattered electron the direction of the 
spin axis x, w is given by (cf. Eq. (5)) 


—u,y/us=cot (x/2)e"*. (10) 
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The scattering of an initially unpolarized beam 
may be described by considering the beam to be 
composed of two equally intense completely 
polarized beams with opposing directions of 
polarization. This direction is arbitrary and it is 

‘convenient to choose it so that the direction of 
polarization in each beam is the same after one 
scattering as before. That is, from (5) and (10) 


u,/us=B/A. (11) 
From (9) we find 
B/A=xie'*. (11a) 


Choosing the scattering plane as the plane ¢=0, 
we have 


B/A= +14, (11b) 


which means that the direction of polarization is 
normal to the scattering plane. The two polarized 
beams may now be considered separately. 

The amplitude of the once scattered beam is 
then 


u3/A=u,/B=f¥ig (12) 


and the intensities for the two scattered beams 
are given by 


2= | f\?+\g|?Fa(f*g—fe*). (12a) 


Thus after one scattering the beam is polarized 
with a polarization given by 


64 =1( fg*— f*g)/(|f|?+\21*)- 


The existence of a polarization therefore depends 
not only on the existence of g but also on a 
difference in phase between f and g. 

The scattered intensity averaged over the spin 
directions is from (12a) 


(12b) 


(9) =| f|2?+/g/%. (12c) 


Since the two beams are treated independently 
the intensity for the double scattering is obvi- 
ously equal to (compare (1)) 


Do( 3) Ho( F2)[ 1+ (5( 01) 5(d2))! cos ¢ |. (12d) 


III. DEPOLARIZATION BY MULTIPLE ELASTIC 
SCATTERING 


We consider that the scattering in each target 
consists of many small angle scatterings (mul- 





tiple scattering) plus one large angle scattering 
so that the total deflection is +/2. The small 
angle scatterings taking place before the first 
large deflection do not polarize the initial beam 
and, of course, we are not interested in the de- 
polarizing effect of the small angle scattering 
after the second large deflection. Therefore our 
problem is to determine the depolarization of the 
beam between the two large deflections after it 
has been polarized by the first of these large 
deflections. 

We define the polarization p in terms of 
the probability amplitudes for the two spin 
directions : 

p= |A|?—|Bl*. (13) 


Then the change of polarization per collision is 
bp=(|us|*— | u4|*)/(|us|*+|us|*)—p. (14) 


From (13), (9), (7) and (12c) we find for the 
change in polarization averaged over the azimuth® 


bp= —2p|g|?/(\f|*+\8!*)- (15) 


Multiplying by the number of collisions with 
scattering angle 3 which take place in a distance 
x, in the first target and a distance x, in the 
second target, which from (12c) is 2rN(x 
+x2)(|f|?+|g!*) sin ddd, we find for the net 
polarization 


P(x, X2) 
=exp| ~4evini+sy) |g? sin oao| (16) 
0 


if the polarization is complete initially. Here N 
is the number of scattering atoms per unit vol- 
ume. Averaging over the effective thickness of 
each target for 45° incidence (and emergence) 
we have 

1 


vid v2d 
=— dxf dx2P(x, X2), (17) 
2d*J/ 4 0 


where d is the actual thickness of the target. Thus 
we find for the degree of polarization 


as 18 
-( yd ) (18) 


‘ There is no correlation between the direction of the 
spin axis and the azimuth of the scattering plane. 
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where 


y=4v2en f |g|* sin ddd. (18a) 
0 


From the above treatment it is apparent that 
P represents the ratio of the polarization of the 
beam just before the second large deflection to 
that just after the first large deflection. Therefore 
the asymmetry is now given by 


26 = P2dbyott- (19) 


In the following we shall refer to 1—P as the 
degree of depolarization. 


Calculation of degree of depolarization 


The calculation of g may be carried out by 
formulating the scattering problem in terms of 
the phases of the scattered (partial) waves at 
infinity. We denote the radial wave functions of 
the third and fourth components of the Dirac 
wave function by Gx. For the two types of solu- 
tions of the wave equation, that is, for the total 
angular momentum quantum number j=/+} 
and /—}, x is given by 


j=l+} 
j=l-}. 

At large distances r from the scattering atom 
G. has the asymptotic form 


x=l for 


rere (20) 


Ge~1/r cos (kr—(1+1/2)r+ nx). (21) 


Defined in this manner, the phase «x vanishes for 
the free particle.® 

Mott! has shown that superpositions of the 
partial waves, Gx multiplied by a suitable 
spherical harmonic, may be found which have the 
asymptotic form (8) with g(#) given by 


g(8) =i/2kEP:(cos 8)[e%—e2*-1-1], (22) 


where 

P}(cos 8) =sin 8(d/d cos 3)P;(cos 8) (22a) 
and P; is the ordinary Legendre polynomial. 
With this expression for g the integral occurring 


* For the free particle G,~r™+J|,,4;(kr), where J is the 
Bessel function 
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in (18a) becomes 


f igi? sin odd 


=2/k? >> (U(l+1)/(2/+1)) sin? 6m, (23)° 
1 


where 
5:= m1 — 9-1-1. (24) 


To determine the phase difference 5n; we con- 
sider the differential equation obtained for Gx by 
eliminating the radial function of the first and 
second Dirac components from the well-known 
pair of simultaneous first-order equations.’ The 
equation obtained may be written in the form 


@G' /dr?+Q(r)G’ =0, (25) 
G=(2me+E— V)-'G’/r. (25a) 


E is the kinetic energy of the electrons and V is 
the scattering potential. 


x(x+1) 2V(E+mc?) 
gf re 
x+1 dV/dr 
r 2mc?+E—V 
1 (dV /dr)*+2(2mc?+ E— V)d?V/dr’ 
4 (2mc?+E—V)? 





O(r) = 








For the scattering potential V we must use a 
screened Coulomb field since at the small angles 
in which we are interested, the scattering by the 
atomic electrons is just as important as the 
nuclear scattering. We shall take 


V=—(e?Z/rje—"”*, (27) 


where 6 is the screening length (see Eq. (39) 
below). 

We may solve (25) by the W. K. B. method. 
If rv; is the (largest) root of Q(r) we have for 
r>?ri, 


const. 


G, =———(2mc?+ E— V)-! 


rQ} 
X cos lf Oldr—r/4| (28) 


m1 


7 Cf. e.g., H. A. Bethe, Handbuch der Physik, Vol. 24/1, 
p. 311. 
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In order to evaluate the integral in (28) we note first of 
all that the last two terms in (26) are small compared to 
the third term. As is evident from (26) the relative magni- 
tude of the small terms will be largest for small r. In this 
case we have 


dV/dr~V/r, #V/dFr~2V/r. (29) 
Then in the last term the important part is that involving 


' @V/dr*, viz.: 


(dV /dr)?/4m2d V /dr?~ V/8mce 
S(Z/137)(kh/mc)1/8IKK1, (29a) 

since for 100 kv kh/mc~1. In (29a) we have evaluated V 
at the minimum ,, i.e., r:~1/k. (See Eqs. (33a) and (34a) 
below.) The ratio of the fifth to the third term is 

4(A/mer)* h(kh/mc)*1/P <1/P. (29b) 
Similarly, the ratio of the fourth term to the third is less 
than 1// which is also small compared to one. 

The third term is therefore the leading term of those 
involving the potential V. This term will be comparable 
with the free particle terms ~# only for the smallest 
values of /, viz.: 


third term __ 


<2Ze m a(R ‘ 
first term ri) eke 


eo kb =) eke (29c) 


where Ex is the ionization potential for the K shell; 
Ex ~40 kv for Au. 


We may now expand Q in powers of the small 
terms. We set 


Q=¢+9q’, (30) 
where 
g=k*?—x(u+1)/r?+2kee*'*/r, (30a) 
ko=e®Z(E+mce?) /hec?. (30b) 
Then 
én: = f (Q,3 aes Qt_.)dr 
=} ["r\@i-a-0er (31) 


"1 


+4 = . . 
=— f q-'(dV/dr)dr/r, 
E+2mcv ,, 


since g is the same for both values of x. 
Denoting the integrand in (31) by J we use the 
identity 


f tar=f tair+(f lar ler), (31a) 


where Jo=J(ko=0) and ro=17;(ko=0) is the root 
of g for ko=0. We denote the contribution of the 








first integral to the phase difference by én)“ 
and that of the last two integrals by 6n;“); the 
latter will be important only for small / (cf. 
above). When the last two integrals are 
considered together, the difference between the 
integrands is appreciable only when the term 
2koe~*'*/r is important compared to the free 
particle terms in (30a). This is the case for 
r~r,~ro. Then, since we consider / small, in 
I—Iy we may set the exponential e~"/*=1 and 
neglect the term 1/6 in dV/dr since these two 
neglections partially compensate each other and 
the error thus induced is ~(ro/b)*~(l/kb)* 
(1/20)? for energies of the order 100 kv, (cf. (33a) 
below). 
With 
bn. = bn, +46n,, (32) 
we have 


1+} 


"dr 
E+2mc*? k =f ~(-+ ene 


In the following we shall replace x(x+1) by 
(1+ 4)? for both values of x, this procedure being 
appropriate for the W. K. B. solutions. Then 


ro= (1+4)/k. (33a) 


The contribution to the phase difference 4n;“ 


is given pt 
= dr 
Af 5 ~¢-\(bs) 
Hi ‘(ko=0 
0 aay 


re 


6,6 = ———— 





bm, PS 
E+2me* 


It may be noted that r; differs very little from ro: 
r, = ro( 1 —ko/k*ro) = rol 1 —(Z/137])) kh/me }. (34a) 


The integration in (33) is carried out in 
the following manner. With the substitution 
r=r, cosh x the integral in (33) becomes 


1 @ 
—{ sech? x(1+ cosh x)e %*" “dx, 
0 


where \=1ro/b=(1+4)/kb. Integrating the first 
term by parts we find 
Ak eZ je 
bn, =— ——_f cosh xe °° *dx. (33b) 
1+4 E+2mc 
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From Heine's integral representation of the 
Hankel function,’ the integral in (33b) is equal 
to —2/2H,™(ir). Then we obtain 


OZ Hy (id) 
8g)? 2 — —— —_— (35) 
2 6 E+2me? 


The integral (34) may be evaluated in terms 
of elementary functions. We find 


(= “) hina (36) 
(1+4)? E+lme° 


r/2<O<r. (36a) 





bn, = 


tan O= —(/+})k/ko, 


From (35) and (36) we note that the phase 
difference 64n; is small for all 7. Since both én,“ 
and 4n:“? decrease monotonically with increasing 
l, the largest 4n; will be 5n). For energies of the 
order 100 kv and for any reasonable screening 
length 6 the argument of the Hankel function in 
(35) is small for small /. 


noe a 


—. (small /) (37) 
E+2me?] 1+} 


For /=1 and E=100 kv (37) and (36) give 
én, =0.116 and bn; =0.176 so that 6n,;=0.29. 
Therefore we may replace sin dy; in (23) by 6m 
for all 1. For very large 1, >kb the phase differ- 
ence behaves like e~'/**. For small /, kb the 
summand in (23) is comparatively small because 
of the factor (1+ 1)/(2/+1). For the important 
values of /, 1~kb the summand varies very slowly 
and we may replace summation by integration. 

For the quantity y occurring in the polariza- 
tion factor P (compare (18)) we find from (18a), 
(23), (35) and (36): 


2 
= ~1ten(——) 
X (log 2kb— C—$+2¢S1+¢°S2), (38) 
where C is Euler’s constant = 0.577 and 
eZ E+mc* 


=— . (38a) 
he EX E+2mc?)! 





* Jahnke-Emde, Tables of Functions, p. 218. 
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mw @ df, ")(4X) 
-> —, (38b) 
2 (44)? 
S2:= > 07(1+4)-'. (38c) 
1 


The sums S,; and S:; may be evaluated numeri- 
cally. For 6 we use the Thomas-Fermi atomic 
radius 


b=h?/meZ' =1.23X10-* cm for Au. (39) 


Then for E=100 kv we have kb=21, ¢=1.05, 
S,=1.98 and S:=1.76. Then 


7 = 890 cm~". (40) 


For the thickest foil used by Dymond, d=2.5 
X 10-§ cm and 
yd =2.2X10™. (41) 


Since the main contribution to the sum 5S, (and 
to S: as well) arises from small* ] where HH,“ (a) 
~2/xx, it follows that S,; and therefore y is 
practically independent of the screening length 6. 
In addition, for the energies of interest, y is 
insensitive with energy and in general yd will be 
of the order 10-*. For such small values of yd the 
degree of depolarization 1—P (Eq. (18)) is just 
yd. Hence multiple elastic scattering produces a 
negligible depolarization of the electron beam. 
To resolve the discrepancy between theory and 
experiment a degree of depolarization of about 
95 percent or more is necessary, i.e., yd should 
be > 4. It is very unlikely that the approxima- 
tions made in the above derivation can introduce 
an error of more than a factor two or three, 
certainly not a factor 100. 

The explanation of the small depolarization we 
have found lies in the fact that while several 
small angle collisions occur the depolarization in 
each collision is too small to give an appreciable 
total depolarization. (Compare further Ap- 


pendix §1). 


IV. INELASTIC SCATTERING 


In this section we shall investigate the effect 
of inelastic scattering on the azimuthal asym- 


*In the sum 5S, the first term contributes about half 
and the first five terms about 90 percent of the total value 
of the sum. In S; the contribution of the first term is 
75 percent and of the first five terms about 99 percent 
of the total value of the sum. 
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metry of the doubly scattered electron beam. In- 
elastic scattering may reduce the asymmetry in 
two ways: (1) a simple reduction of the energy 
before the large deflections take place and (2) 
depolarization caused by change of spin of the 
incident electron in the inelastic collision. 

The first effect would reduce the asymmetry 
since the large deflections would then occur at 
lower energies at which the polarization is indeed 
smaller (cf. reference 3). However, the rate of 
energy loss'® in Au is two Mev per g/cm? and 
with Dymond’s foils the total energy loss is only 
one kv. Therefore this reduction in the asym- 
metry is entirely negligible. 

We consider now the inelastic scattering in 
which the incident electron changes its spin. If 
we denote the differential cross section for this 
process by ®,’ sin ddd, the depolarization per 
collision is 2%,’/#, where ®, is the total scattering 
cross section for excitation of the mth atomic 
state. In a manner similar to that used in the 
derivation of (18) we have for the degree of 
depolarization : 


1—e-74\? 
1~P=1-( ) . (42) 
yd 





where now 


y=42rNE [ ®,'(3) sin ddd. (42a) 
"eo 


The summation in (42a) is to be taken over all 
excitations. 

To calculate the differential cross section we 
may use the Born approximation. That this pro- 
cedure will yield a sufficiently good result may 
be concluded from the fact that the values of the 
angular momentum of the incident electron for 
which the inelastic scattering is important are 
larger than the corresponding momenta for 
elastic scattering. Since larger angular mo- 
mentum implies smaller phases of the scattered 
partial waves, it follows that the Born approxi- 
mation is better for the inelastic scattering than 
for the elastic scattering—and for the latter the 
Born approximation is quite accurate even for 
the values of atomic number and energy with 
which we are concerned (see Appendix §2). 





” See e.g., H. A. Bethe, reference 7, p. 522. 


The differential cross section (per solid angle) 
in Born approximation is given by" 


4e* 
*  C(e—k’)?+(E—E'/he)*}? 


E+mc?*\? 
x( ) asl (43) 


h*c? 








Ma= f drut exp (iq°1T;)(@o+@-@j)tm, (43a) 
| 


ao=s(k)s*(k’), a= —s*(k’)as(k) (43b) 


and q=k—k’, q=2k sin 8/2. (43c) 


Here E, k and E’, k’ are the energies and wave 
vectors before and after the scattering; s(k) and 
s(k’) are the amplitudes of the Dirac plane waves 
for positive and negative spin, respectively; u» 
and u, are atomic wave functions for the initial 
and final state, respectively, the components 
being Schrédinger-like wave functions since very 
little energy is given to the atom; the Dirac 
operator @; acts on the spin coordinate of the jth 
electron. In (43) and (43c) we have used the fact 
that the energy loss is small: k~k’, E=E’. 

In §4 of the appendix it will be shown that the 
second term of (44a) is much smaller than the 
first term over practically the entire range of 
integration over 8. Then the retardation term in 
the denominator of (43) may be neglected and 
the matrix element (43b) becomes 


M,=4@0F,, (44) 


where F,, is the “generalized form factor” 





Fem drut exp (1q-Tj)Um. (44a) 
ae 

With ao=-— sin de‘? (45) 
2 E+mc* 


the differential cross section (43) becomes 


1 e 2 
o,'=-(_—_) cot? (8/2)|F.|*. (46) 
4\ E+2mc? 


 H, A. Bethe, reference 7, p. 495. 
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Summing over all excitations we have, using the 
Hartree model” 


LX | F.l?=L(1—| Fal*), (47) 


ngim m 


where F,, is the atomic form factor for the mth 
electronic state 


F.= f dr exp (iq:t)|Ym(t)|?, (47a) 


wherein y, is the wave function of a single 
atomic electron in the mth state. The integration 
in (47a) is now to be taken over the coordinates 
of a single electron. The summation on the right- 
hand side of (47) is to be taken over all occupied 
states. 

For large g the atomic form factor F,, is neg- 
ligible due to interference. We may set 


F,,=0 for q>1/bn, (48) 


where b,, is the radius of the mth electron shell. 
For smaller g we may expand the exponential in 
(47a) and obtain 


Fa=i—ig' fdr ra|yml?, G<1/bn (48a) 


where 7, is the component of r in the direction 
of q. The term linear in r, vanishes from sym- 
metry. The matrix element of 7,2 may be set 
equal to b,,?. Thus we have 
2p 2 / 
1—[Falt=[?" aiadad (49) 
1 g>1/b» 


Actually 1—| F,,|* is less than the smaller of the 
two values on the right-hand side of (49). 

For the evaluation of the total cross section for 
excitation of the atomic electron in initial state 
m we have 


® e 2 
fo! sin ddd = (=) Lot), (50) 
9 E+2mc? m 
where 


Im 
Iy=bat f cot? (8/2)g¢? sin ddd =(50a) 
0 


% The transitions which are forbidden by the exclusion 
principle (state initially occupied) have not been sub- 
tracted. Therefore the cross section is overestimated. 
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is the contribution to the scattering from small 
angles and 


sin &,,/2=1/(2kb,,). (50b) 


The contribution from the large angle scatter- 
ing is 


n={ cot? 3/2 sin ddd. (50c) 
Im 


Carrying out the integrations in (50a) and (50c) 
we find 


In carrying out the summation over m occurring 
in (50) we obtain the geometric mean of the shell 
radii 6,, which we may replace for order of mag- 
nitude purposes by ), the screening length intro- 
duced in Section III (cf. Eq. (39)). That is, 


+ (o+11) =4Z log 2kd. (50e) 


For the quantity y we obtain from (42a) 
vy =4v24NZ(e?/E+2mc?*)* log 2kb. (51) 


For E=100 kv and the foil thickness d=2.5 
X< 10-5 cm we find 


y=5.0cm™", yd=1.3X10™. (52) 


We see then that the depolarization by inelastic 
collisions is negligible. We observe that the de- 
polarization caused by inelastic scattering is 
smaller than that caused by elastic scattering by 
about a factor 1/Z (cf. Eqs. (38) and (51)). This 
can be understood from the fact that the small 
angle scattering gives a small contribution (cf. 
Eqs. (50a), (50c)) to the depolarization. At large 
scattering angles, it is well known that the total 
elastic scattering is about Z times larger than the 
inelastic, and the depolarization in one scattering 
process is the same for elastic and inelastic 
collisions at a given scattering angle, (see Ap- 
pendix, Eqs. (107), (107a)). 


V. EXCHANGE SCATTERING 


A change in the polarization of an electron 
beam can take place as a result of exchange 
scattering only if the exchanged electrons have 
opposite spins. Since the spin forces are weak for 
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the atomic electrons the spin of each electron will 
not change in the capture and emission processes. 
This also has the important consequence, that 
such exchange processes as we are considering 
can take place only with atomic electrons in in- 
complete shells. Exchange between the incident 
electron and an atomic electron of opposite spin 
in a closed shell will result in two atomic elec- 
trons having the same set of quantum numbers 
and this is forbidden by the Pauli exclusion 
principle. Thus in Au only the 6s valence electron 
can participate in the type of exchange scattering 
which leads to depolarization. 

In contrast to the case of depolarization by 
elastic and inelastic scattering, the depolarization 
per collision due to exchange scattering will not 
be small but the met depolarization will be ex- 
tremely small because very few exchange col- 
lisions of the type considered will take place. This 
is due essentially to the strong interference be- 
tween the wave functions for the valence state 
and the continuum. 

To obtain a quantitative estimate of the im- 
portance of exchange scattering we make the 
following (necessarily) rough calculation of the 
exchange cross section. The differential cross 
section is given by™ 


$..(8) sin ddd = (8r%a9")—'| M.x|? sin ddd, (53) 


where d)9=f*?/me® is the Bohr radius. A factor } 
has been inserted in (53) because the spin of the 
incident electron must be opposed to that of the 
valence electron and this is true in half the cases. 

The matrix element for exchange is essentially 
given by 


Mea f drdrovz"(t)xw"(t2) 
X |r—tz|~'W2(tz)xx(r),- (54) 


where r, rz denote the coordinates of the elec- 
trons initially in the continuum and valence 
state, respectively. ~z is the wave function for a 
single electron in the valence state and moving 
in the field of the rest of the atom. x, and x, 
are the wave functions of the incident and out- 
going electrons, respectively. We take 


xu(r)=exp (¢K-r), xx-(tz)=exp (tk’- rz), (55) 


4 See e.g., H. A. Bethe, Ann. d. Physik 5, 325 (1929). 


k and k’ representing the initial and final wave 
vectors. 

Because of the rapid oscillation of x the 
contribution to the matrix element will be small 
in the region where yz is slowly varying, but 
will be appreciable near the nucleus where ¥z 
also varies rapidly. We therefore represent the 
valence state wave function by 


¥z(r) = vz(O)e-*, B = Z/ ao. (56) 


With the wave functions (55) and (56) the 
matrix element (54) may be evaluated to give 
approximately 


32x? 


M x = ——|2(0) |* (57) 
Rp? 


(B*+q")* 
where q=k—k’ (cf. Eq. (43c)). 

The value of the 6s wave function ~z at the 
nucleus may be obtained from considerations of 


the type used in the theory of hyperfine struc- 
ture.'* For s electrons we find 


Z(1+2)*«(4, Z;) B*« 


’ 
rajn*? xZ2n*3 


(58) 





¥2z*(0)= 


where z is the degree of ionization (here =0), 
Z; the effective nuclear charge (here =Z), n* is 
the effective quantum number and x«(}, Z) is a 
relativistic correction factor given in Eq. (246) 
of reference 14. For the 6s state of Au, n*=1.21 
and «= 2.0. 

From (53), (57) and (58) we find for the 
differential cross section 


x? R,\* ae? 
..(8) =128 (—) , (59) 
Zn *6 E (1 + £7) 








where §=(2kao/Z) sin 8/2 and R,=e'm/2h* is 
the Rydberg energy. 
From (59) we obtain for the total cross section 


32% x? sRy\*/a0\? 
wn AC) 
3 n**Y\E Z 
x [1—(1+4k*/p*)-*]. (60) 


4 See e.g., H. A. Bethe and R. F. Bacher, Rev. Mod. 
Phys. 8, 82 (1936), Chapter VIII. 
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For energies of order 100 kv, k/8~1 and the 
square bracket in (60) is practically equal to one. 

For the degree of depolarization we have 
(compare (42a)) 


yd = 2v2 Nd®,,*** 


_ Aves x? 
~(- ‘\ (- *) wa. (61) 
n*®\ E 


For E=100 kv, d=2.5X10~ cm we find yd=2 
x<10-", an entirely negligible depolarization. 
This result is due to the rarity of exchange 
collisions. Since most of the collisions are elastic 
we see by comparison of (60) and the total 
cross section for elastic collisions (Appendix, 
Eq. (112)) that the fraction of collisions which 
are of the exchange type is 





?,,*°* 


420) Qe 
,,*°* 3 n*é 


Thus for E = 100 kv only one in 4X 10'* collisions 
is of the exchange type. 





APPENDIX 


§1. Depolarization per collision at small angles 


In order to see that in each elastic collision 
only a smail depolarization occurs, it is ad- 
vantageous to consider the angular dependence of 
the depolarization. The depolarization in a de- 
flection at angle 3 was given in Eq. (15) as 


—dbp/p=2|g\*/(|f|*+/¢!*). 


The function g is given by (22). The Faxen- 
Holtzmark expression for f is 


(101) 


f() =i/2k¥P,(cos 8) 
0 


 [(L+1) (e%!— 1) +-U(e%-"-1-1) J. (102) 


From (21) and (28) the phase 7x is given by 


m=Lim( f Qtdr— [ euar), (103) 


where the subscript 0 refers to Z=0. We can 
expand the integrand in terms of the most 
important potential term and find for both 
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values of x (compare (30))" 





© dr eZ ri 
eo ah, f —e-""Qul=—— —Ho'(id). (103) 


» v 
For the sake of simplicity we first eliminate 
n-1-1.. We can see very easily that the phase 
difference 9;— ~:-,1= 6; is small compared to 9; 
for all 1. The ratio 6n:/m is largest for ]/=1, 
From (103a) and ((37) et seg) we find for the 
usual values of the constants, 6n;/n,;=0.10. 
Therefore in (102) and (22) we expand e***" and 
f and g become 


f=(i/2k) XO (21+1)Pi(cos #)(e*"—1), (104) 
0 


g=—(1/k) DP :\(cos 8) bne2. (104a) 
0 


Now since most of the collisions take place 
for small scattering angles we are interested in 
the values of f and g for 3 small. Making the 
appropriate expansions we have from (22a) 


f=(t/2k)> (21+1)(e?*"—1), (105) 
0 


g=(0/2k) S1(1+1)dn,e?*". (105a) 
0 


Since m, and én, behave like e~'’** for large /, the 
series (105a) converge and we find from (101) 
that at small scattering angles the depolarization 
per collision is small of order J*. This result is, 
therefore, true for any scattering field for which 
the scattered intensity is finite for small angles. 
For scattering in a pure Coulomb field, the 
series (105a) will, of course, diverge. In this 
case, for small #,'* 


| f | =(2e?Z/hv)(1/8°R), (106) 
|g| =(e*Z/hv)(E/E+mc*)(1/dk), (106a) 


18 The error incurred by this expansion is small for the 
important values of /(~kb). The error is greatest for /=1; 
with E=100 kv the phase from (103c) is 2.62 whereas 
a numerical integration of (103) gives the exact value 
m =2.87. The iene between the exact value of n, 
and the value given by (103a) decreases rapidly with 
increasing /. 

16 Reference 1, Eq. (6.8). 
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so that the depolarization per collision for small 
deflections is of order # for all scattering fields. 

For the numerical value of the depolarization 
in one elastic collision we find for 100 kv (see 
Eqs. (109) and (109a) below) 


2\2/f|?=0.01382. (107) 


For inelastic collisions we have, of course, the 
same result. The depolarization per inelastic 
collision is from Section IV (cf. (45) and remarks 
preceding Eq. (42)), 


26,’ /@, =2|ao|?=}(E/E+mc?)? sin? 3 
~0.0138?. (107a) 


The small depolarization found in the fore- 
going may now be understood as follows. The 
depolarization depends on the product of the 
scattering cross section and the depolarization 
per scattering (101). The former is large only for 
small # and just in this region the latter is 
small. 


§2. Validity of the Born approximation 


In the foregoing the statement has been made 
that the Born approximation is accurate for the 
elastic scattering even though the Born expansion 
parameter e’?Z/hv~1 for the energies and atomic 
number in which we are interested. We shall now 
show that this is so not only for the scattered 
intensity |f|?+|g/* but also for the depolariza- 
tion~ |g|*. For large scattering angles the Born 
approximation is known to give a result for the 
scattered intensity which is in agreement with 
the exact formula of Mott (nuclear scattering). 

In the Born approximation the phases 7; and 
én; are assumed small and for f and g we have 
from (105) and (105a) 


f= — (2/B)E0+4)m, (108) 


gp =(0/2k) > (1+4)*5m. (108a) 
0 


Replacing sums by integrals and using (103a) 
and (35) we find!’ 


7 It is obviously consequent in the first Born approxi- 
mation to use for the phase and phase difference the 
limiting forms for small scattering potential, hence 


bn: = bn. 


ATT 9‘? (dd) dA 


1/kb 
= —x(e*Z/hv)bH, (i/kb) 
= 2(e*Z/hv)kb*, (109) 


fa= —rikb*(e?Z /hv) 


since kb>1; 
r eZ E i ~ 
£a= ——kb?*— aaa) of M277," (td) dd 
4 he \E+2mc*? L/ kb 
rie’Z E ig 
= —— — ———) —H,(i/kb) 
4 hc \E+2mc?7 k 


= (e*Z/hc)(E/E+2mc*)'kb?d. (109a) 


Of course gp and fs depend on the screening 
length 6, but their ratio is essentially inde- 
pendent of the screening and in agreement with 
(106) and (106a) as it should be, cf. Sauter."* 

In order to test the Born approximation for 
our problem, the exact value of the sums 
occurring in (104a) were found numerically.'* 
For kb=100 and eZ/hv=1 corresponding to E 
about 100 kv we find 


|f| =1.81X104/k, |g) =0.167XK10*8/k. (110) 
The Born approximation yields 
fe=2.0X10*/k, ge=0.182K10'8/k. (110a) 


The use of the Born approximation therefore 
overestimates f by only ten percent and g by 
9 percent. For the scattering the error is 21 
percent and for the depolarization 18 percent. 
The screening radius }b has been taken rather 
large in the above example. For smaller screening 
lengths the phases for a given / are smaller and 
the Born approximation better. 

This result is perhaps not surprising insofar as 
it is well known that the condition for the 
validity of the Born approximation eZ/hoe<1 is 


18 F, Sauter, Ann. d. Physik 18, 61 (1933). 

1*For the calculation of these sums the part 4m; of 
the phase difference may be neglected since this part of 
én: contributes to g the amount 

1/2Z\* E+me oo, 

a~3(< 25 — feelin 
The absolute value of the sum is of order or less than rkb 
and for the numerical case considered in (110) g,:300/k 
as compared to g~17000/k from 4n;. 
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a sufficient but not a necessary condition. In the 
problem considered here it is possible to see 
qualitatively the reason for the fair agreement 
between (110) and (110a). In (104) and (104a) 
the main contribution to each of the sums arises 
from those values of / for which the phase n, is +/2 
or smaller and in this region the assumption of 
small phases is sufficiently well fulfilled. 

It may be pointed out that the Born approxi- 
mation is not accurate for the calculation of the 
polarization of the electrons, but only for the 
depolarization (or scattered intensity). The 
polarization as was mentioned in the foregoing, 
depends on the relative phase of f and g whereas 
the depolarization depends on the value of g 
itself. The relative phase of f and g is in fact 
zero in the first Born approximation®® (in con- 
trast to the value of g itself) and at least the 
second approximation is necessary for the calcu- 
lation of the asymmetry. 


§3. Total number of collisions 


a. Elastic collisions.—The average number of 
elastic collisions taking place in a foil of thick- 
ness d is 


ne=2V2eNd f .:(8) sin ddd, = (111) 
0 


where a factor v2 has been inserted because of 
the obliquity of the path. With @=/f|? and 
f given by (104) we have in the Born approxi- 
mation 


buiot— de f .i(8) sin ddd 
0 
=4r(e*Z/hv)*b*. (112) 


With d=2.5 X10-* cm and b=h?/me*Z! we find 
that on the average 40 elastic collisions take 
place in each foil. 

b. Inelastic collisions.—The total number of 
inelastic collisions in each foil is* 


2vix Nd : > or ~— 113) 
1, = T aan "te ’ ( 
mv? " A, Aa 





where A, is an energy of the order of magnitude 


*°N. F. Mott, Proc. Roy. Soc. 124, 425 (1929). Also 
F. Sauter, reference 18. 


E. ROSE AND H. A. 


BETHE 


of the ionization potential and », the number of 
electrons in the mth shell. For Au we may take 


A,= 6 volts for the 6s shell 
san SS ear 
=100 “ “ “ 5s, 5p and 4f shells. 


The contributions to the sum in (113) from 
the 6s, 5d and the 5s, 5p and 4f shells are 1.85, 
4.95 and 1.7 reciprocal volts, respectively. The 
contribution from the remaining shells will be 
about 0.5. The total number of inelastic collisions 
is then n;~16. This number may be in error by 
a factor 3 or so in either direction but in any 
case it is indicative of the fact that several 
collisions can take place in the foils under 
consideration. 


§4. Relative magnitude of transition “current” 
and “density” terms 
In Section IV it was stated that the current 
term in (43a), viz: 


a-J=a: far exp (1q-Tj)aju, (114) 
7 


makes a contribution to the inelastic scattering 
which is much smaller than that due to the 
density term aoF, (cf. (44), (44a)). Since the 
relative magnitude of the current term is largest 
for small scattering angles we may expand the 
exponential and obtain 


F,,=1Q° Tn, (115) 


where 


ram f Lristnds (115a) 
| 


is the dipole moment corresponding to the 
transition u,,—«,. For the current term we have, 


with a;= —f;/c 
J=i(AE/ho)t., (116) 


where AE is the energy loss. Therefore the ratio 
of the current to density term is 

a-J AEa-r, AE (a 
= = —, (117) 
doF, hcaoq-tn hcay g 
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since the direction of r, is not significant. From 
(43b) we find 


|a| =Aclq?—(AE/hc)* }*/2(E+mc*). (118) 
With ao given by (45) we find 
a-J AE 


<p=———. (117a) 
aoF, E sin 3 





The right-hand side of (117a) is appreciable only 


for 


dS AE/E~#,, (118) 


where #,, is the angle introduced in (50b). In 
the average #,,~1/kb so that p~1 for 


d~1/400=0.14°. 


For scattering angles 8>9d,", a-J<aoF,. 
Therefore over practically the entire range of 
integration over # the current term need not be 
considered at all. 
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The unsymmetrical broadening of the Q branch toward higher frequencies in the parallel type 
bands near 1100 cm™ in the absorption spectrum of deutero-formaldehyde is taken as a mani- 
festation of convergence toward lower frequencies rather than, as is usual, in the case of infra-red 
bands, toward higher frequencies. This phenomenon is accounted for by taking into account an 
interaction between rotation and oscillation of the Coriolis type arising from the accidental 
degeneracy of the three vibration frequencies, belonging to three different symmetry classes, 


falling near 900 cm™, 1000 cm™ and 1100 cm", 


I. INTRODUCTION 


HE absorption spectrum near 10y in the 

deutero-formaldehyde spectrum consists of 
three bands which badly overlap. Two of these 
are characteristic of oscillations of the electric 
moment normal to the axis of symmetry of the 
molecule, while the third is characteristic of a 
vibration parallel to the axis of symmetry. The 
first two of these consist essentially of a group of 
prominent lines protruding above a background 
which at present cannot be resolved spectro- 
scopically, and the third contains P, Q and R 
branches. It has already been pointed out else- 
where! that the spacings between the prominent 
lines in the perpendicular bands in this region are 
noticeably greater than the spacings between the 
lines in a similar band near 4.6z and it has been 


suggested that the explanation for this is un- 


1E. S. Ebers and H. H. Nielsen, J. Chem. Phys. 6, 311 
(1938). 


doubtedly the same as that for the similar effect 
observed in two overlapping perpendicular bands 
in the spectrum of ordinary formaldehyde near 
8.0u.2 The effect is here quite satisfactorily ex- 
plained by taking account in the energy of a 
Coriolis interaction between the two frequencies. 
In the third of the above bands referred to in the 
spectrum of deutero-formaldehyde the con- 
vergence is in the direction of lower frequencies 


*H. H. Nielsen, J. Chem. Phys. 5, 818 (1937). The 
type of interaction between rotation and oscillation 
discussed recently by W. H. J. Childs and H. A. Jahn 
(Nature 141, 916 (1938)) in connection with the methane 
spectrum is entirely similar to the one discussed in this 
work relative to HCO. The formaldehyde case is especially 
interesting in that both perturbing frequencies are optically 
active so that the convergence in opposite directions in 
the two bands, described by C. and J., can here actually 
be observed. In the case of the methane molecule the 
perturbing frequencies are two- and three-fold degenerate, 
respectively, while in this case they are nondegenerate. 
The importance of such terms in the interpretation of 
spectra was already emphasized by E. B. Wilson, Jr., in 
one of his early papers on this subject. 
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rather than in the direction of higher frequencies 
as is usually the case for bands in the infra-red. 
Although the P and R branches are too incom- 
pletely resolved to satisfactorily study the con- 
vergence it may be established that this is true 
by observing that the Q branch is unsymmetri- 
cally broadened toward higher frequencies. For 
this to be so, what corresponds to the classical 
moment of inertia of the molecule must be smaller 
in the upper state than it is in the ground state. 
An explanation suggests itself when one reviews 
the case of ordinary formaldehyde. Here, due to 
Coriolis forces, the two nearly coincident fre- 
quencies which, however, belong to different 
symmetry classes, suffer perturbations to their 
accompanying rotational levels in such a manner 
that the band of lesser frequency is strongly con- 
vergent toward higher frequencies while the 
other band is strongly convergent toward lower 
frequencies. In this case we have, not two, but 
three different oscillations of nearly coincident 
frequencies, and belonging, here also, to three 
different symmetry classes. We shall investigate 
in this note what the nature of the interactions 
may be here between rotation and oscillation to 
see whether or not they may give rise to effects 
in the spectrum such that the Q branch of the 
parallel type band will be unsymmetrically 
broadened toward higher frequencies and at the 
same time the type of convergence retained 
which is characteristic of the overlapping perpen- 
dicular type bands in the spectrum of ordinary 
formaldehyde. 


II. SOLUTION TO THE QUANTUM MECHANICAL 
PROBLEM 


The three oscillations which give rise to the 
three bands in question are the following ones 
given in the order of ascending frequencies: the 
deformation frequency of the molecule in and out 
of its own plane; the frequency caused by the 
deformation in the plane of the molecule of the 
O—C—D bond angle and the frequency caused 
by the deformation in the plane of the molecule 
of the CD, bond angle. The first two of these will 
give rise to changes in the electric moment per- 
pendicular to the axis of symmetry, while the 
third will produce an alteration of the electric 
moment parallel to this axis. As a simplified 





V4 Vs Vs 


Fic. 1. The three modes of oscillation which the molecular 
model may execute. 


model of the molecule which will approximate to 
the actual case when it is oscillating in these 
three modes we shall adopt one where the C and 
O atoms and the C and D atoms, respectively, 
are rigidly joined together, an assumption which 
is justified when further we adopt force fields of 
the valence force type, bearing in mind that the 
frequencies in which we are interested are all 
deformation frequencies where no changes in the 
interatomic distances take place. This simplifica- 
tion is in reality tantamount to suppressing the 
other three modes of oscillation of the molecule. 
The xz plane has been taken as the plane of the 
molecule, the z axis having been made to lie 
along the axis of symmetry and the origin of the 
system to coincide with the center of gravity of 
the molecule. In this system the coordinates of 
the C, O and D particles when they are at rest, 
will be respectively, (0,0, 23°), (0,0, 2,.°), and 
(x,°, 0, 2:°) and (x2°, 0, 22°) where x,°= —x,° and 
z;°=22°. The three modes of oscillation which this 
model may execute are described in Fig. 1 and 
are respectively identified with the frequencies 
of the molecule, v4, v5 and vg where v4> v5 > v¢. 

The quantum mechanical equation for this 
molecular model, oscillating harmonically in 
these three modes and rotating, has been de- 
veloped with the neglect, however, of all terms 
in the moments of inertia which arise from de- 
formation of the molecule. This can be justified 
on the basis that these will almost certainly con- 
tribute amounts to the energy small in compari- 
son with the interaction terms in which we are 
interested, namely those arising from the Coriolis 
forces. The derivation of the quantum mechan- 
ical equation will be omitted here since the 
method is a perfectly straightforward one. The 
equation which is obtained may be divided into 
zeroth and higher order terms of which only the 
two first sets will here be set down: 
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(Hy— Eo)Wo= (hv/2) {2Eo/hy+ (vs/v)(d?/ dn? — 9?) 
+ (v4/v)(82/ag? — ¢°) + (v6 /v) (82 /dE* — £*) 
+u?(2e/h)*(P2+P,*)+(29/h)*(Ao/Co)P.*) 
— 2d3u?(A o/Co)[(ve/vs)'nd/dE 

— (vs/v—)*€0/An \(2i/h)P.}Yo=0, (1a) 

Hy = —whv {dal (va/v6)'€0/E 

— (ve/v4) 0 /0E(2ri/h) |P.+del (5/4) 'f0/On 
—(v4/v5)'nd/0E ](2et/h)P,}. (1b) 


In the above u?=(h/47*Aoyv), 


h,=[2m(M.+ Mp) tan a 

+4m(M.o— Mop) tan a ll (M.+Mo)uteus'}", 
he=[—4m(M.+ Mo)e tan a 

+4m(M.e— Mod) tan a \[(M.+Mo)piius' }', 
As=(2m+4M.€6+4M dp) (uitus'), 


P,, P, and P, are the usual angular momentum 
operators’ (functions of the eulerean angles @, ¢ 
and y) conjugate to the x, y and s axes, Ao, Bo 
and Cy are the principal moments of inertia. In 
the above m, M. and M, are respectively the 
masses of the deuterium, carbon and oxygen 
atoms, o, p, « and 6 are constants defined as 
follows: 


(mz — m(2,;° —23°))(M_-(24 — 23°) + 2m(2,° —23°))— 
(m/M>,)(M.23°— (M.+2m)(2,°—2;°)) 
XK (M(2¢ — 23°) +2m(2,°—23°)) 
m(24° — (2,9 —23°+x;° tan a))(M.(24° — 23°) 
+2m(z,°—z2;°+x;,° tan a)) 
(m/M,)(M.23°—(M.+2m)(2;°—23°+x,° tan a)) 


xX (M.(2.8 —23°) + 2m(2;° —2;°+x;)° tan a)). 


The variables », £ and ¢ are related to the dis- 
placements of the D, C and O atoms from equi- 
librium in the x, y and z coordinates in the 
following manner : 


n= (4% yivs/36h)'(x1+-x%2— (x3/€) +(x4/8)), 
£= (492 ugv6/36h)'(yitye— (¥3/0)+(¥4/p)), 
f= (492 yov,/4h)*(2; +22) cot a, 


7E. C. Kemble, Fundamental Principles of Quantum 
Mechanics (McGraw-Hill Co., 1937), p. 232. 


where: 
ui = ((2m+4M 62 +4M 95") 

+2m(1+-4e+4e*) tan*® a), 
pe = (2m(M,.+ M,)(1+ tan?’ a) 

+4m? tan* a)(M.+Mp»)-', 
us = (2m+4M.07+4M op") 


and may be regarded as reduced masses respec- 
tively to be associated with one of the three 
modes. 

The last term in (1a) is an interaction term 
between rotation and oscillation. Normally it 
would be retained as a part of (1b) with the other 
similar terms. Because of its coefficient (Ao/Co), 
which for D,CO is of the order six, this term will, 
however, be considerably larger than the other 
interaction terms so that for values of the rota- 
tional quantum numbers up to about ten it will 
be actually of the same order of magnitude as the 
rotational energy. It is therefore included as a 
part of H». The second-order terms of H, which 
are not given here, consist of terms which are 
essentially the squares of the oscillational angular 
momentum components along the x, y and z axes. 
Their contributions to the energy of the molecule 
will be entirely independent of the rotational 
quantum numbers. Since they are therefore not 
an interaction energy we shall not compute their 
value. 

Following the procedure suggested in the 
earlier work referred to herein we take as a solu- 
tion to the zeroth-order Hamiltonian for the 
ground and first excited states v4, v5 and v., the 
functions : 


1(0) = | (642? ui pousrareve/h*)* 

Xexp— (n?+£+¢") /2}0(6, J, K, M) 

Xexp (iKe+iMy), (2a) 

¥(vs) = N(vs) | —th[(ve—v5)/2 

+ ((ve—vs)?/4+(KAsh/4x*Co)*)*) 

X (2§ exp— (n?+£?+$7)/2) 

+ (Kygh?/42? Co) (v5/v6)* 

X (2n exp — (n? ++ ¢7)/2)} 

x 0(6, J, K, M) exp (tKg+iMy), (2b) 
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¥( v6) = N( v6) | (iK Ash? /44?Co) (v6 /v5)! 

X (2E exp — (n?+8?-+n*)/2) +h[(ve—vs)/2 

+((ve—vs)?/4+(Krsh/4?Co)*)! ] 

X (29 exp — (0? +?+ 5") /2)} 

x 0(0, J, K, M) exp (iKe+iMe), (2c) 
where 
Niu? = | (KAgh?/42?Co)?(ve/r;) +h?L(v6—v5)?/4 

+ (Krgh/42?Co)*)* }?} (1628 nwo svarsve/h?)! 


j and k taking the values 6 and 5 and 5 and 6, 
respectively. 
¥(¥4) = { (160? wiper arsre/h*)? 
X (2 exp— (n?+£?+¢*)/2)} 

x 0(0, J, K, M) exp (tKe+iMy). (2d) 
0(0, J, K, M) exp (iKe+iMy) in the above 
functions is the normalized eigenfunction for the 
symmetric rotator problem ; K and M are integers 
less than, or equal in the limiting case, to J. The 
functions (2) give the energies of Ho as diagonal 
matrices except for a small contribution arising 
from (Av/v) which may be regarded as a left over 
from the term representing the Coriolis inter- 
action. This term may actually be included as a 
part of H; and its contribution to the energy 
taken into account at a later point if this is 
desirable. 


The zeroth-order energies for the molecules in 
the corresponding states are readily found to be: 


E(0) = (h/2)(vs+va+v6) +(J(J+1) —K?) 
X (h?/82*?A 9) +K*h?/8x7Co 
E(v4) = (h/2)(vs+3va+v6) +(J(J+1) —K?) 
X (h?/82°A o( 1 —a—B))+K*h?/82°Cy 
E(v5) =h(vs+v4/2+ 46) +h((ve—vs)?/4 
+ (KXsh/4x?Co)*)'+ (J(J+1) —K*) 
X (h?/8*A o(1+8))+K*h?/8x°Cy 
E (ve) =h(vs+4/2 +6) +h((ve—vs)*/4 
+ (Kysh/4x*Co)*)'+ (J(J+1) —K*) 
X (h?/82°A o(1+a))+K*h?/82°C, 


(3) 


where a and 8 may be regarded as correction 
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terms to Ao which to this approximation are, 
however, both equal to zero. 

To evaluate the corrections to the energy con- 
tributed by H; we may without hesitation make 
use of the nondegenerate perturbation theory and 
proceed in the usual straightforward manner to 
evaluate E, and Ez, the first of which vanishes. 
We shall include in E, only the second-order 
contributions from 7/7, as already suggested, and 
of these we may actually neglect all the terms 
except those where the frequency differences 
(vs—vs5) and (v,—¥v¢) occur in the denominators, 
since all others will be small when compared to 
these. The term (Av/v)nt of which we have al- 
ready made mention will contribute an amount 
entirely too small to justify its evaluation and it 
also will be omitted. The remaining terms of E; 
when computed and added to Ey change the 
expression (3) only in this respect that now a 
and 8 are no longer zero but take the values: 


a = (h/827A 9) {Ay?( KAgh*v5/44? Cove) 
X ((vatve) /vatve!)?+ro7h?[ (vg —v5)/2 
+ ((vg—v5)?/4+(KAsh/42°Co)?)! }? 
X ((va+vs)/vatvs!)?} { (vs—vs) (h?L(ve—vs) /2 
+ ((vg—vs)?/4+(Krgh/42?Co)?)! }? 
+ (Kyr3h?/42?Co) (v6/v5))} 
and 
B= (h/8x7A 0) {dy7h?L(ve—vs)/2+((ve—vs)?/4 
+(Ky3h/42?Co)?)* (ve +r.) valve)? 
+)o?(KAsh?/42?Co) (v6/v5) ((va tvs) /vatys')?} 
X {(A?[(ve—vs)/2+((ve—vs)?/4 
+ (Kyr3h/42?Co)*)* }? 
+ (Kr3h?/42?Co)*(v6/v4))(vs—ve) } 7, 
both of which are quantities greater than zero. 
We shall not here consider it necessary to ob- 
tain the selection rules, these having been dis- 
cussed elsewhere, but shall be content with stat- 
ing what these are. The transitions which may 
take place and which are of interest to us in 
this discussion are the following ones: from 
the ground state to the states vs; and »., with 


AJ=+1, 0 and AK = +1; from the ground state 
to the state », with AJ= +1, 0 and AK =0. 
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Ill. THe Errect oF THE INTERACTION 
ON THE SPECTRUM 


We are now in a position to inquire into what 
will be the effect on the spectrum caused by these 
interaction terms. We may accomplish this by 
obtaining expressions for the frequency positions 
of lines in the spectrum using the above energy 
relations and selection rules. For the frequencies 
ys and vg only those lines arising from the 
transitions AJ=0, AK = +1 are of interest, the 
other lines being at present spectroscopically 
unresolved. The frequency positions of these lines 
in the two bands are given by: 


v= (vs t+v6)/2—(h/8x*)(1/Co—1/Ao) 
+[(ve—vs)?/4+(Kdsh/4x*Co)?} 
+(J(J+1)—K*)ah/8r*Ay 
+(Kh/4x*)(1/Co—1/Ao), (5) 


a relation which, except for the presence of a 
small term ah*/8xr*Ao, confirms the one used by 
Ebers and Nielsen in their work on ordinary 
formaldehyde. 

For the frequency », we obtain as frequency 
positions of the permitted lines, according as 
AJ=+1, AK=0 and as AJ=0, AK=0, re- 
spectively, the two expressions: 


v=vg(Jh/8x°Ao)(2+a4+ 8) 
+(J?—K*)(a+B)h/8x°Ao 
and v=4+(J(J+1)—K?*)(a+B)h/8x*Ao. (6) 


Equation (5) has already been used by Ebers 
and Nielsen to explain the structure of the two 
bands »; and », in the spectrum of formaldehyde, 
and as we have pointed out, a structure analagous 
to this is found in the corresponding two bands 
in deutero-formaldehyde. Our endeavor from the 
outset was, however, to account for the unsym- 
metrical broadening toward higher frequencies 
of the Q branch of the parallel band in this region. 
We are therefore particularly interested in the 
second of Eqs. (6) which sets forth the com- 
ponents comprising the Q branch. Since a and 8 
are both positive quantities it becomes clear at 





once that, consistent with measurement, the Q 
branch must spread out toward higher frequen- 
cies. The lines of Eq. (6) have been plotted with 
their relative intensities for comparison with the 
branch as actually observed. This diagram is not 
reproduced here, but is in reasonably good agree- 
ment with experiment both in shape and in 
breadth. 

In determining the values of the moments of 
inertia of the formaldehyde molecules, it has 
always been assumed that the spacings of the 
perpendicular C—H and C—D vibrations near 
2875 cm and 2160 cm are the correct ones to 
use, an assumption which has been justified 
principally on the basis that this vibration is a 
single perpendicular oscillation. While this is 
indeed true, it may be pointed out that the other 
CH and CD vibrations, parallel to the axis of 
symmetry, and therefore normal to the first ones, 
lie near 2780 cm and 2056 cm=, respectively. 
One may then advance the argument that these 
frequencies might interact so as to perturb the 
accompanying rotation levels to such extent that 
also these spacings are unsuitable directly for 
determining the values of the moments of in- 
ertia. The extent of this perturbation can be de- 
termined by referring to Eqs. (3). Since there is 
only one perpendicular frequency in this region, 
any other perpendicular frequencies which the 
molecule may have must lie so far away that the 
difference vg—vs will be large compared to 
(Ki3h/8x°Co). The term under the radical may 
therefore be expanded and E(v.) takes the form: 


E(v6) = (h/2)(vat+-vs+3e) + (J(J +1) 
— K*)h*?/8x*Ao(1+a)+K*h?/8r*Co(1+7), (7) 


where y is equal to \°3h4/8x*Co(vs— ve). ¥ will 
under these conditions be small and we have 
already seen that a is not large. To a fair ap- 
proximation then A» and Cy remain unchanged 
as the molecule makes a transition from the 
ground state to this excited state so that these 
spacings may safely be used in making numerical 
determinations of the moments of inertia of the 
molecule. 
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Several arguments are given which support the view that the temperature distribution of 
the vapor in the high pressure mercury discharge tube is not that given by Adams and Barnes 
(line A, of Fig. 1) but has a shape as indicated by the line B, of Fig. 1. 





N TWO recent articles Barnes and Adams! 
have concluded that the temperature distri- 
bution of the vapor in the high pressure mercury 
arc is as indicated by the line A, in Fig. 1, in 
contradiction with the distribution B, which I 
proposed some years ago.’ A, indicates that the 
temperature of the vapor is substantially constant 
over the whole cross section except in the 
boundary layer which is about 0.8 mm thick (in 
analogy to the stationary film of Langmuir).* The 
curve B, was found? from the condition that 
= 2ardr 
0 T(r) 
must have a given value (pressure and mean 
density were known) and from a plausible as- 
sumption on the dependence of the coefficient of 
heat-conduction on temperature (the absolute 
value was not used in the calculation). As the 
convection current in the vertical mercury tube 
is substantially laminar and vertical, it has only 
little influence on the horizontal temperature 
distribution (we do not consider the case that the 
convection is turbulent).‘ 

The distribution A, agrees with the first con- 
dition as well. By assuming a constant tempera- 
ture over the total cross section one finds from 
the same data which led to the distribution B,, a 
temperature of 555/0.223 = 2500°K,? which is in 
agreement with the ditribution A, of Barnes and 
Adams. The difference between the two curves 


thus arises from the second condition according 
to which I used the equation of heat-conduction 





1B. T. Barnes and E. Q. Adams, Phys. Rev. 53, 545 
(1938); E. Q. Adams and B. T. Barnes, Phys. Rev. 53, 
556 (1938). 

2 W. Elenbaas, Physica 1, 211 (1934). 

#1. Langmuir, Phys. Rev. 34, 401 (1912). 

‘ W. Elenbaas, Physica 3, 484 (1936). C. Kenty, J. App. 
Phys. 9, 53 (1938). 


(energy streaming through a cylinder one cm 
high is 2xrAd7T/dr) whereas Adams and Barnes 
use the theory of Langmuir. Furthermore, the 
distribution of the electron temperature, ac- 
cording to Barnes and Adams, is as substantially 
indicated by the line A,., whereas I supposed the 
electron temperature (curve B,) to be only 
slightly higher than the gas temperature, so that 
calculations can be made on the assumption that 
a temperature equilibrium between the gas 
molecules and the electrons exists in the arc. 

The difference between the two conceptions 
is so considerable that we wish to state the 
following objections against that of Adams and 
Barnes : 

(a) Barnes and Adams find the intensity ratio 
for a number of spectral lines to be the same in 
the axis and 8 mm aside (inner tube diameter 36 
mm) ; from this fact they conclude the tempera- 
ture to be the same in both points. However, this 
does not agree with my own measurements® and 
with more extensive recent measurements of 
Fischer and K6nig,® who find that the intensity 
of a line drops more rapidly if its initial level is 
higher. This indicates that the excitation temper- 
ature is a maximum in the axis. As the points 
through which Barnes and Adams have to draw a 
straight line, whose slope determines the temper- 
ature, scatter much, it is possible that the 
temperature difference between the axis and 8 
mm aside, has escaped their attention. 

(b) Even if the temperature in the arc path 
were uniform, we may not conclude that the 
temperature over the total cross section (except 
the boundary layer) is constant (diameter of 
measured portion 16 mm and PQ about 34 mm). 

(c) Langmuir’s theory of the film of stationary 
gas around a heated wire cannot be applied to the 


* W. Elenbaas, Revue d’optique 15, 343 (1936). 


* E. Fischer und H. Kénig, Physik. Zeits. 39, 313 (1938). 
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TEMPERATURE 


high pressure arc, where the energy is conducted 
to the surrounding solid part (tube wall) instead 
of from the solid body (the wire) to the sur- 
rounding gas, as in the case studied by Langmuir. 
This is caused by the fact, that in the case of the 
wire there can exist a practically uniform temper- 
ature field (namely room temperature or tem- 
perature of the bulb) beyond the stationary film 
(by the upward stream of the cold convection 
gas), but inside the mercury tube there must be a 
temperature gradient when energy flows to the 
walls (unless the gas is in a rapid and irregular 
motion in the part PQ (Fig. 1)); this is, however, 
not the case as stated under (d). The temperature 
distribution A, would therefore only be possible 
if no energy were transformed into heat in the 
region between P and Q, but only on the edge 
PQ. This is very unlikely, as the arc burns in the 
central part of the tube. 

(d) Convection is observable on small particles 
which are moving along with the vapor.‘ They 
indicate a chiefly laminar streaming which goes 
upwards in the central part of the discharge and 
downwards in the surrounding part. The con- 
vection cannot cause a uniform temperature 
across the diameter PQ, as for that purpose only 
radial components, which practically do not 
exist, would be useful. If the temperature distri- 
bution were as indicated by A, the convection 
would be upward in the whole part PQ. Actually 
the particles are moving downward at a much 
greater distance from the wall. This phenomenon 
cannot be caused by electrical forces, as it is the 
same with alternating and direct current. 

(e) For the energy which is conducted to the 
wall, Barnes and Adams calculate from the 
Langmuir theory 32 watts per cm of length. The 
total radiation according to Barnes and Adams is 
therefore as a maximum 40—32=8 watt or 
20 percent of the input. The total radiation in the 
above case, however, is greater. For the total 
radiation S per cm of length of the high pressure 
mercury discharge in a quartz tube we found in a 
wide range of diameters and pressures :’ 


S=0.7(L—10) watt, (1) 


where L is the input per cm of length. Applying 
this to the tube of 20 mm diameter, mentioned by 


*W. Elenbaas, Physica 4, 413 (1937). 
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Barnes and Adams, on which I did the measure- 
ments,? we find S=0.7(40—10) = 21 watt (L=40 
and not 43.8 as stated in reference 2, where the 
energy losses at the electrodes had not been 
taken into account). Therefore, at most only 19 
watts could be conducted to the wall. The 
interpretation given of Eq. (1) leads even to a 
smaller amount of energy which is conducted to 
the wall, namely 10 watts per cm of length (we 
supposed that of the input LZ, 10 watts are 
conducted to the wall and L—10 watts are 
radiated, of this radiation only 70 percent passes 
the quartz wall). This energy of 10 watts is in 
agreement with measurements on the influence of 








TUBE AXIS 




















— RADIUS 


Fic. 1. Gas temperature A, and electron temperature 
A, as functions of the distance to the axis according to 
Barnes and Adams. B, and B, are the corresponding 
temperatures according to the author. 


noble gases on the mercury arc*® and with the 
behavior of the gradient as a function of LZ at a 
given mean vapor density.® If we put this energy 
of conduction in the formulas of reference 2, we 
find for the coefficient of heat conduction \ for 
mercury vapor 2.1X10-’7**. Extrapolating to 
476° (203°C) we find a value of 2110~* for 4, 
whereas a value of 18.5 10~* has been measured 


* W. Elenbaas, Physica 3, 219 (1936). 
* W. Elenbaas, Physica 2, 757 (1935). 
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by Schleiermacher.'® Thus the loss of 10 watts by 
conduction is in agreement with the slope of curve 
B,. As the convection currents in this tube were 
laminar, the energy loss by convection may be 
neglected (only at the first part from the bottom 
there appears an additional loss caused by the 
heating up of the colder convection gas). 

(f) According to Adams and Barnes the energy 
of 32 watts, which is conducted to the wall per 
cm of length, can be dissipated just at a temper- 
ature of 900°K (outer diameter of the tube 22 
mm). With a temperature of 900° they calculate a 
loss of 5.3 watts by conduction and convection 
and a loss of 25.6 watts by radiation. They 
obviously treat the glass wall as a black body. 
As the maximum of the black body radiation of 
900°K lies at 3.24 and the transmission of most 
glasses is good till about this wave-length, the 
emissivity of glass must be considerably smaller 
than one at this temperature. If we assume that 
of the 21 watts which are radiated by a quartz 
tube (formula (1)), 18 watts pass the glass wall 
(the tube wall was of thin, ultraviolet-trans- 
mitting glass), the tube wall has to dissipate 22 
watts per cm of length. Subtracting 5.3 watts for 
conduction and convection, 17 watts remain for 
radiation. At 875° (actually measured tempera- 
ture) an emissivity of 75 percent is needed hereto. 

(g) A maximum value for the difference be- 
tween the electron temperature 7, and the gas 
temperature 7, may be calculated according to 
Druyvesteyn" in the following way: The energy 
lost by the electrons per cm’ and per sec. by 
elastic collisions with the Hg atoms, is according 
to Cravath” (with a few immaterial neglections) 


Out on —) (1 =) (2) 
nia’ { m / M T.) 





nm, and m, are the number of electrons and Hg 


10 A. Schleiermacher, Wied. Ann. 36, 346 (1889). 

" Unpublished. See also R. Mannkopff, Zeits. f. Physik 
86, 161 (1933). 

# A. M. Cravath, Phys. Rev. 36, 248 (1930). 
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atoms per cm’, ¢ the radius of the Hg atoms and 
m and M the mass of the electron and the Hg 
atom. We can find , from the mobility equation: 


3n.efPG 


{= ——____-__, 
80°n,(2emk1.)! 


(3) 


where i is the current densiiy and G the gradient. 

Substituting (3) in (2) and putting na, 
=mo(273/T,), if mo is the number of atoms per 
cm’ at a pressure of one atmos. and at 0°C (the 
pressure in the above-mentioned tube was one 
atmos.), we find as 0<iG 


mr? 7k 3e°G?M 
=) ~=)- @ 
ff T,] —128%(273)%o'k?no2m 


If we take for ro? the value 5X10-'* and for 
G=8 volt/cm (measured value), we find 0.17 for 
the right side of Eq. (4). This gives (7./7,) <1.15. 
As only a fraction of the input will be transferred 
to the atoms by elastic collisions, the real differ- 
ence between electron temperature and gas tem- 
perature will be smaller than 15 percent. The 
only uncertain quantity in (4) is o. The value 
5X10-" for ro? may be wrong to a factor 3." 
Taking xo* three times smaller we find (7./T,) 
<1.83, whereas Barnes and Adams need at least 
a value of 2.4 (6000/2500). If, however, the elec- 
tron temperature is 6000° and the vapor tempera- 
ture 2500°, the total excitation is practically 
caused by the electrons, so that the energy trans- 
ferred to the atoms by elastic collisions, will be 
smaller than 41G (the measured total radiation is 
50 percent of the input, Eq. (1)). A ratio of 2.4 for 
T./T,, therefore corresponds to a value of 
mo? <0.8X10-5, which seems to be a very 
improbable value. 

For all these reasons we believe that it is much 
more likely that the temperature distribution is 
as indicated by the line B, than as indicated by 
the line A,. 





13 W. Elenbaas, Physica 5, 573 (1938). 
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A Study of Beta-Brass in Single Crystal Form 


Wayne WEBB 
Department of Physics, University of Iowa, Iowa City, Iowa 
(Received September 21, 1938) 


The resistances of twenty single crystals of brass in the beta-phase have been measured as a 
function of temperature from room temperature to about 550°C for the purpose of studying the 
order-disorder phenomenon. The results obtained are compared with the theory of Bragg and 
Williams. The resistivities at 25°C as a function of composition place the limits of the beta- 
phase at 51.3 and 55.5 atomic percent copper. A study of the thermal e.m.f. of the crystals 
against copper shows only a slight discontinuity in the second derivative of the e.m.f.-tempera- 
ture curve at the temperature of complete disorder. The e.m.f. curves themselves are quite 
smooth. Young’s modulus, EZ, for ten of the crystals, measured at room temperature, gives the 
elastic constant Sj, (44.8X10~" cm*/dyne) and the combination of constants, 2Si:.+Si, 
(0.028 x 10-" cm*/dyne). A graph of 1/E against the orientation function, (?m*+-m'n*+-n'*?P) 
is linear. From it the ratio of the maximum to minimum E is found to be 8.94. Below the elastic 
limit the stress-strain relation is linear with no permanent set nor hysteresis. 





INTRODUCTION 


ESEARCHES on the properties of metal 

alloys capable of forming superlattices have 
been stimulated, particularly in this country and 
England, by the recent theoretical developments 
of Bragg and Williams,! and Bethe.? Improve- 
ments in details of the theory have been made 
by a number of workers.* On the experimental 
side specific heat and electrical resistance have 
been determined as functions of temperature and 
have been correlated with the state of order 
existing.* An excellent review of both theory and 
experiment has been given recently by Nix and 
Shockley.§ 

Beta-brass is one of the simplest of alloys 
showing the order-disorder phenomenon. It is the 
beta-phase of the copper-zinc system, with a 
composition in the neighborhood of the one-to- 
one ratio of atoms. The crystal structure is 
body-centered cubic. In the disordered state the 
two kinds of atoms have equal probabilities of 
being at any lattice point. In the ordered con- 
dition the zinc atoms occupy the cube centers 


'W. L. Bragg and E. J. Williams, Proc. Roy. Soc. 
A144, 340 (1934). 

*H. A. Bethe, Proc. Roy. Soc. A150, 552 (1935). 

*R. Peierls, Proc. Roy. Soc. AlS4, 207 (1936); N. F. 
Mott, Proc. Phys. Soc. 49, 258 (1937); T. S. Chang, Proc. 
Roy. Soc. London A161, 546 (1937); E. E. Easthope, Proc. 
Camb. Phil. Soc. 33, 502 (1937). 
need Sykes and F. W. Jones, Proc. Roy. Soc. A157, 213 

6). 

‘F.C. Nix and W. Shockley, Rev. Mod. Phys. 10, 1 

(1938). 


and the copper atoms the cube corners or vice 
versa. This fact, though surmised from other 
data, was definitely established by Jones and 
Sykes in an x-ray study in which they were able 
to obtain superlattice lines for the ordered alloy*® 
by a judicious selection of the x-ray wave-length 
used. 

The present research includes a study of cer- 
tain properties of single crystals of this beta-brass 
in the ordered condition (at room temperature), 
and of the order-disorder transformation by 
following the changes in certain properties as 
the temperature is gradually raised to a point 
above the critical temperature. 


PREPARATION OF CRYSTALS AND DETERMINATION 
OF ORIENTATION 


The beta-brass from which the crystals were 
grown was prepared from very pure commercial 
materials. The copper was cathode copper from 
Anaconda Copper Company and was stated by 
the Company to be 99.95 percent copper. The 
zinc was Bunker Hill from Platt Brothers and 
Company and contained not more than 0.01 


percent impurity.’ 


*F. W. Jones and C. Sykes, Proc. Roy. Soc. Al61, 440 
(1937). 

7 An estimate of the impurities present in this brand of 
zinc is given by H. E. Way, Phys. Rev. 50, 1181 (1936). 
Actually the writer's material was all taken from a slab 
which showed distinctly less lead than Way's, but other 
impurities in the same amount. See Way’s footnote 14, 
which refers to the slab used by the writer. 
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The crystals were grown from the molten 
state in a horizontal mold by the method em- 
ployed for some time in this laboratory.* The 
preliminary polycrystalline casting was prepared 
by melting together in a crucible® suitable pro- 
portions of copper and zinc and pouring the 
molten metal into the crystal mold. This was 
done in a furnace in which a nitrogen atmosphere 
was maintained throughout. The crystal was 
grown in a second furnace also in nitrogen. 
The orientation of the crystal was not pre- 
determined by the use of a nuclear crystal, since 
it was found that, by making the end of the 
casting pointed, a single crystal would start. 
There seemed to be no particularly preferred 
orientation, so that thus a suitable variety of 
orientations was assured.” The crystals were of a 
trapezoidal (nearly square) cross section, about 
0.7 cm? in area, and 14 cm long. All the crystals 
were grown with a temperature gradient of 
about 6°C/cm and a rate of growth of approxi- 
mately 6 cm/hr. 

During the cooling following the formation of 
the crystals the temperature was checked at 
600°C for twenty-four hours and then allowed to 
fall slowly to room temperature. The purpose of 
this anneal was to secure homogeneity of com- 
position throughout the crystal. 

After completion of the measurements the 
crystals were analyzed for copper content by the 


5C. A. Cinnamon, Rev. Sci. Inst. 5, 187 (1934). 

® The crucible and molds were made from lava obtained 
from the American Lava Company. The writer is indebted 
to the company for furnishing the unfired material and 
later firing the molds, etc., after they were made. 

1 Tt is doubtful if a seed crystal could be fused to the 
casting with a flame since attempts to fuse two pieces of 
beta-brass together resulted in evaporating so much zinc 
that the joint became alpha-brass. 


TABLE I. Orientation of crystal lattice, measured angles. 


standard method of electrolytic deposition from 
a sulfate solution. In most cases the sample 
analyzed was about 3 g, 1.5 g from each end of 
the crystal. In three cases separate analyses were 
also made for the two ends. None of these showed 
a difference between the two ends greater than 
0.1 percent of copper content. The compositions 
obtained extended throughout the beta-brass 
phase and somewhat beyond at room tempera- 
ture as can be seen from Table II. 

The crystals, as removed from the mold and 
without further treatment, showed in a very 
marked manner the phenomenon of strong light 
reflection in definite directions, precisely as de- 
scribed by Bridgman" for a number of metal 
crystals. These reflections could be improved in 
some cases by heating the crystals from two to 
three hours in a hydrogen atmosphere. Etching 
the surfaces damaged the reflections and made 
them more diffuse. 

These reflections were used to determine the 
orientation of the length of the specimen relative 
to the space lattice by a slight variation of 
Bridgman’s method. The general idea of the 
method is to locate the normal to each reflecting 
plane by causing a fixed beam of incident light 
to be reflected back on itself. When the normals 
to a number of planes have been located, the 
angles between the planes are then easily com- 
puted. To identify the planes it is then necessary 
to see whether each plane can be assigned a set 
of Miller indices so that measured and computed 
angles between every possible combination of 
pairs of planes form a completely consistent 
scheme. Table I shows the degree of consistency 


" P. W. Bridgman, Proc. Am. Acad. Arts and Sci. 60, 305 
(1935). 























REFLECTED 
1 2 3 4 5 6 7 8 Beam No. 
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obtained for eight of the reflections from one 
of the crystals. The numbers across the top and 
down the side refer to arbitrary numbers 
assigned to the reflections. In parentheses are the 
Miller indices which give the type of plane. 
The angles in the body of the table are those 
observed between the normals to the planes in 
the corresponding column and row. It can be 
seen that rarely do the observed angles disagree 
by more than one degree with one of the geo- 
metrically possible angles given in the table. 
More complete tables involving all observed 
reflections, were fully as consistent. Besides the 
two types of planes shown in the table, planes 
of type {112} were frequently observed. 

The above measurements were made by use 
of an especially constructed goniometer, which 
need not be described in detail. It permitted one 
to obtain a particular reflection from the whole 
of the crystal length without disturbing the angle 
setting. Only specimens for which the direction 
of the reflection remained constant the full 
length were considered to be single crystals and 
satisfactory for use in subsequent measurements. 

Naturally, once the positions and types of 
several lattice planes are determined, it is 
merely a matter of geometry to compute the 
orientation. The orientations so obtained are 
given in Table II by stating the two larger of 
the three direction cosines of the long axis of the 
specimen with reference to the three principal 
crystal axes (cube edges). While it has been 
stated that the angles between the lattice planes 
did not involve an error greater than one degree 
the error involved in the orientation angles may 
be as great as three degrees. This difference is 
caused by inability to set the specimen in a 
precise manner on the goniometer. 


RESISTIVITY AND RESISTANCE 


Two sets of resistance measurements were 
made on each crystal; one to obtain an accurate 
value of the resistivity at room temperature, and 
the other to study the temperature dependence 
of resistance. A potentiometer method employing 
a Leeds & Northrup, Type K2, potentiometer 
was used to make both sets of measurements. 

For the resistivity measurements needle points 
were used as potential contacts which made 


TABLE II, Resistance characteristics of the crystals. Limits 
of the beta-phase are indicated by the arrows. T’ and T”’ 
represent the temperatures between which the resistance- 
temperature curves were nonlinear. 








Two LaRGEst 





DiREcTION 
CosINES 
Atomic | aX10 e”" —p”" 
PERCENT PER x 108 
CopPpEeR © OHM-CM T’(°C) | T" (PC) i ay 





50.22 | 3.54 4.40 193 | 490 | 0.806 | 0.560 
50.70 | 3.85 4.31 250 | 480 807 | .562 
50.81 | 3.47 4.41 200 | 490 622 | .561 
50.89 | 3.41 4.37 165 | 485 868 | .396 
$1.25 | 3.61 4.46 207 | 490 914 | .294 
51.63 | 3.50 3.62 
$2.13 | 3.43 3.76 250 | 475 : d 
52.16 | 3.53 4.03 205 | 487 908 | .357 
53.20 | 3.02 3.76 160 | 475 -788 | 488 


53.34 | 3.32 4.35 160 | 485 876 439 
53.65 | 2.96 3.46 175 | 490 914 334 
53.77 | 3.16 3.04 185 | 490 801 435 
53.81 2.42 4.30 150 | 480 813 583 
54.25 | 2.87 2.68 212 | 480 717 633 
54.26 | 3.07 3.09 210 | 480 716 716 
54.60 | 3.12 2.57 240 | 465 .739 584 
54.89 | 3.08 3.34 150 | 475 .900 453 
55.50 | 2.96 2.99 200 | 485 779 624 
55.85 | 3.00 2.67 200 | 490 933 291 


























possible a precise measurement of the length. 
The average cross-sectional area between the 
potential contacts was determined by finding the 
volume of this portion of the crystal. This was 
done by sawing out this section, at the com- 
pletion of the work, and getting its mass and 
density. 

To measure the resistance as a function of 
temperature the crystal was placed in a furnace 
in a holder equipped with three thermocouples, 
one near the middle and one near each end of 
the crystal. The potential and current leads were 
spot-welded to the crystal. Readings were taken 
as the temperature of the furnace was raised at 
a rate which varied from a degree per minute for 
the lower temperature range to about one-fourth 
degree per minute at the higher temperatures." 
Readings were also taken with the temperature 
decreasing gradually to 400°C, and a final reading 
at room temperature. There was a slight tem- 
perature gradient along the crystal of such 


12 These runs were made at this slow rate to eliminate any 
lag between the temperature indicated by the thermocouples 
and the actual temperature of the crystal. The agreement 
for rising and falling temperatures shows that this was 
accomplished. It shows further no lag between order and 
temperature ; such a lag was not expected, however, since 
it is known that the relaxation time for this alloy is very 
short. 
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Fic. 1. Resistivity as a function of composition. Curve A 
at 25°C; curve C at 500°C. Curve B, Haughton and 
Griffiths’ curve; they state their values are for room 
temperature. The ordinate scale for C is on the right. 


amount that the temperatures at the potential 
contacts usually differed by about three degrees, 
never, however, exceeding ten degrees. Thermal 
e.m.f.’s caused by this difference were corrected 
for by taking potential readings with the current 
through the crystal both direct and reversed; 
also, with no current. 

In curve A, Fig. 1, the resistivity at 25°C is 
shown as a function of composition. Curve B, 
given for comparison, represents the data of 
Haughton and Griffiths" for polycrystalline ma- 
terial at “room temperature.” The two curves 
are seen to agree in a general way. Some of the 
difference is undoubtedly caused by difference in 
purity of materials and to the mechanical treat- 
ment given the samples of Haughton and 
Griffiths. 

The boundaries of the beta-phase occur at the 
end points of the middle straight section of the 
resistivity-composition curve. The writer's data 
(curve A) indicate 51.3 and 55.5" atomic per- 
cent copper whereas Haughton and Griffiths 
estimate 50.7 and 54.2 atomic percent copper 
from their data. Values given by metallurgists, 
obtained usually by thermal and micrographic 
analyses, vary considerably among themselves. 
The exact limits of the beta-phase at room 
temperature are extremely difficult to determine, 
since, while an alloy may be examined at room 


sti S. Haughton and W. T, Griffiths, J. Inst. Met. 34, 
245 (1925). 

™“ This boundary is not closely established by the data. 
The break in the curve might well come anywhere between 
54.5 and 55.5 percent Cu with an equally good fit of the 
curve to the observed points. 
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temperature, the structure observed may be 
typical of a higher temperature if the cooling to 
room temperature has been too fast, or perhaps 
at any finite rate. Some typical values for the 
beta-phase limits in atomic percent copper 
are : Gayler :'® 55.0—51.0; Hansen :** 55.6—51.0; 
Saldau and Schmidt :'’ 54.4—49.8; International 
Critical Tables: 54.7—51.2. Sykes and Wilkin- 
son'* obtained phase boundaries by measurement 
of the energy difference of the ordered and dis- 
ordered state as a function of composition. 
Their limits are: 55.3 and 51.1 atomic percent 
copper, in very good agreement with the writer’s, 

It is interesting to note that practically all 
determinations of the beta-phase of brass just 
fail to include the ideal, equiatomic ratio. It is 
only for such a ratio that perfect order, or a 
perfect superlattice can exist. Other compositions 
in the beta-phase may have a maximum order, 
but never perfect order since the number of 
copper atoms is always too large. Sykes and 
Wilkinson in order to find the transformation 
energy for the ideal CuZn, extrapolated their 
data to the one-to-one ratio. The writer's 
resistivity data may be treated in the same way 
by extending the straight middle portion of 
curve A, Fig. 1, to 50 atomic percent copper. 
The value of the resistivity so found is 4.02 K 10-* 
ohm-cm. This is very near the mean value of the 
resistivities of the constituent metals, which is 
3.93X10-*, derived from Way’s'® value for 
Bunker Hill zinc, and the International Critical 
Tables value for pure copper. This shows that 
the completely ordered structure behaves in re- 
sistivity like a pure metal, with value inter- 
mediate between those of the constituents, and 
does not show the enhanced resistivity charac- 
teristics of an ordinary, or random, solid solution. 
The temperature coefficient, as will be seen 
later, is also about that of a pure metal. This 
simpie relation seems to be approximately satis- 


%M.L. V. Gayler, J. Inst. Met. 34, 235 (1925). 

%M. Hansen, Der Aufbau der Zweistofflegierungen 
(Julius Springer, Berlin, 1936). 

17 P. Saldau and I. Schmidt, J. Inst. Met. 34, 258 (1925). 
The specimens of Saldau and Schmidt were annealed 84 
days at 440°C, 

18C. Sykes and H. Wilkinson, J. Inst. Met. 61, 223 
(1937). These authors give a good discussion of the 
uncertainties of these boundaries. 

1* Reference 10. An average value for random orientation 
is given by ps0 and 4po. The value so obtained is reduced to 
25°C, aad Gea the value in the text. 
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fied in other similar cases® but in this case is 
really surprising in view of the very different 
lattice structures of the alloy and either con- 
stituent (alloy, b.c.c.; copper, f.c.c.; and zinc, 
hex.). 

Curve C of Fig. 1 gives the resistivity at 
500°C. These data are taken from the resistivity 
at room temperature combined with the re- 
sistance temperature results, to be described 
later. It is seen that the phase boundaries are 
no longer defined by this curve. All the crystals 
were within the beta-phase at this temperature. 
This curve is similar to those for continuous 
solid solutions, such as shown by Mott and 
Jones.” 

The very regular variation of the resistivity 
with composition shown by curves A and C, 
Fig. 1, indicates that there is certainly no very 
large dependence of resistivity on orientation. 
In view of the fact that the lattice structure is 
cubic it is almost certain that no such variation 
should be expected. 

The curves of Fig. 2 show typical resistivity- 
temperature relations. Both curves are for 
crystals inside the §-phase, one, A, near the 
boundary on the copper rich side, and one, B, 
near the other boundary and very close to the 
ideal equiatomic ratio. These curves have the 
following features: a straight portion at lower 
temperatures followed by a portion of rapidly 
increasing slope, a short “‘tail’’ of decreasing 
slope, and finally another linear portion extend- 
ing to higher temperatures. The curves for all 
the other crystals were similar to these two. 
Without presenting all these curves, the differ- 
ences between crystals may be indicated as 
follows. If the straight line portions are extended 
(curve B), the point 7’ at which the curve de- 
parts from the lower tangent marks the onset of 
detectable disorder, while 7’’ marks the com- 
pletion of the principal disordering process. 
Furthermore, the difference, p’’—p”, is the 
contribution of disorder to the resistivity under 
the assumption that the temperature dependence 
of resistivity, aside from disorder, may be 


* See F. C. Nix and W. Shockley for CuAu and CuPd 


systems. 
“= N.F. Mott and F. W. Jones, The Theory of Properties 
of Metals and Alloys (Clarendon Press, Oxford, 1936), 
ig. 103, page 299 and Fig. 104, page 300. 


linearly extrapolated. Values for such quantities 
taken from corresponding curves for all the 
crystals are given in Table II. This table inclvdes 
also the composition and temperature coefficient 
of resistivity of the ordered state taken from the 
resistance-temperature curves. 

It is seen from this table that the total con- 
tribution of the disorder to the resistivity is 
roughly greater the closer the crystals are to 
the ideal equiatomic ratio. There are, however, 
easily noted exceptions to this general behavior. 
These arise largely because of the irregularity of 
occurrence of the temperature 7’, which is seen 
to range from 150° to 250°C. When it comes at a 
high temperature, less resistance change is 
attributed to disorder, than when it comes low. 
However, for want of a better criterion it seems 
fair to consider the “disordering”’ as starting in a 
noticeable fashion at 7’ as chosen. The general 
trend to lower values of p’’—p” for increasing 
copper content seems reasonable because the 
greater the excess of copper (above equiatomic 
proportions) the less the completeness of the 
order which must be destroyed at the high 
temperature. The alloys beyond the beta-phase, 
at the top and bottom of the table, behave quite 
similarly to alloys just within. These alloys in 
the two-phase field are quite interesting in that 
the disordering process is accompanied by change 
of the alpha- or gamma-brass present to the 
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Fic. 2. Resistivity-temperature. Curve A for crystal of 
55.50 atomic percent copper. Curve B for one of 51.25 
atomic percent copper. The ordinate scale for curve A is on 
the right. 
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Fic. 3. Order versus temperature. Curve A, experimental 
curve for crystal 51.25 atomic percent copper. Curve B, 
Bragg-Williams theoretical curve for equiatomic ratio. 


beta-phase as the resistance-temperature meas- 
urements are made. 

The completion of disorder is far more definite 
than the onset, as shown by the smaller range 
of values of 7’’, 475°C to 500°C. The existence 
of the “‘tail’’ at the upper portion of the curves 
is quite definite and was shown by all the 
crystals. Points on it could be repeated on 
different runs, on the same run for heating and 
cooling, as well as by holding a constant tem- 
perature at a particular point. Only a small 
part of this “‘tail’’ is caused by the slight tem- 
perature gradient along the crystal since the 
width due to this cause should not exceed the 
difference in temperature between the points of 
contact, which was in most instances certainly 
less than 5°C with an observed “‘tail’’ of the order 
of 15°C wide. Thus in no case was there any 
indication of a discontinuity at the critical 
temperature in the slope of the resistivity- 
temperature curve such as is indicated for CusAu 
by Sykes and Jones‘ and by Kurnakow and 
Ageew.™ As far as the writer can tell from the 
figures published by Haughton and Griffiths® for 
a single polycrystalline specimen of brass con- 
taining 47 percent copper by weight, no pro- 
nounced “tail” is present, but on the other hand, 
there is no indication of the very sudden upturn 
of the resistance-temperature curve and the 
consequently very marked discontinuity just 
mentioned. 

The temperature coefficient of resistivity of 
the ordered state is seen likewise to decrease in 


#2 N.S. Kurnakow and N. W. Ageew, J. Inst. Met. 46, 
481 (1931). 
*3 Reference 13, Figs. 2 and 3, pp. 248-249. 
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Fic. 4. Curve A, ordinate scale on left, shows a typical 
thermal e.m.f.-temperature relation. Curve B, ordinate 
scale on right, thermoelectric power-temperature. 


general as the composition varies from the 
neighborhood of the equiatomic ratio toward the 
copper-rich side of the phase. This is to be 
expected because the behavior should depart 
from that of an ordered alloy and approach that 
of a disordered solid solution. From the curves 
of Fig. 2 it is seen that the slopes of the upper 
linear portions are greater than those of the 
lower. This tendency was general, the difference 
in the slopes being greater the greater the 
copper content of the crystals. There seems to 
be no obvious explanation for this unless perhaps 
it be attributed to a gradual disappearance of 
local order after the main disordering process is 
complete. 

A very rough estimate of the dependence of 
resistivity on order has been made by Bragg and 
Williams. They assumed essentially that the 
amount of added resistance was proportional to 
the disorder present. This may be illustrated 
from curve B of Fig. 2 in which the order at any 
temperature is given by the ratio of the ordinate 
segments, ab/ac. Curve A, Fig. 3, is a plot of 
the order obtained in this way, for the crystal 
of 51.25 atomic percent copper. For abscissae a 
reduced temperature, 27/7., is used, in which T 
is the absolute temperature and 7, is the critical 
temperature at which long range order dis- 
appears. As the critical temperature for this 
composition the writer has adopted the carefully 
determined value of Sykes and Wilkinson," 
741°K (468°C). For comparison the Bragg and 
Williams! curve for long range order and for 
equiatomic proportions is shown, curve B, Fig. 3. 
The general agreement is fairly good. The 
assumption that has been made throughout that 
order is complete at room temperature naturally 
causes the experimental curve to lie above the 
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Bragg and Williams curve in the high order 
region. The order shown beyond the critical 
temperature is largely due to the “tail” of the 
resistance-temperature curves, a possible ex- 
planation for which has already been suggested. 
Similar curves for the other crystals are some- 
what as follows: (1) crystals near the equiatomic 
ratio agree very well with curve A; (2) crystals 
of high copper content differ from curve A, and 
lie about ten percent or less above or below it in 
regions of order greater than 0.5, and distinctly 
to the right in the region of rapidly decreasing 
order, below 0.5; (3) crystals with copper less 
than the equiatomic proportion lie very close to 
curve A for high order and either close or 
slightly to the left in the region of lower order. 

The theory of T. Muto*® leads to a connection 
between order and resistivity given by: 


p=pitp2St+psS*. 


A comparison of the above experimental results 
with this expression is not possible because of 
the inability to evaluate the quantities pi, pe 
and ps which are themselves temperature de- 
pendent. The Bragg and Williams approxima- 
tion is equivalent to the first two terms of 
Muto’s expression. 


THERMAL e.m.f. AND THERMOELECTRIC POWER 


The thermal e.m.f. of the crystals against 
copper was measured from room temperature up 


TABLE III. Elastic characteristics of the crystals. 














| Ex10™" 1/E X10" ATOMIC 

tL, m, n) (DYNES/CM®*) (cM®/DYNES) PERCENT Cu 

0.143 3.62 27.6 53.65 
183 4.02 24.9 53.34 
243 6.07 16.5 55.50 
.009 2.29 43.8 52.13 
.193 5.39 18.7 50.89 
.286 8.25 12.1 54.60 
231 | 6.00 16.6 50.22 
.323 14.3 6.97 50.81 
.162 4.18 24.1 52.16 
a iy 21.6 51.25 

| 











*% This deviation could be lessened if the Bragg and 
Williams procedure had been exactly followed. Thus, 
instead of using the experimental resistivity near room 
temperature as typical of the completely ordered state, an 
ideal resistivity would have been used, namely, the average 
resistivity of copper and zinc. 

T. Muto, Sci. Pap. Inst. Phys. Chem. Res. 30, 99 
(1936). 
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Fic. 5. Curves A and B, ordinate scale left, show load- 
strain relationship. Curves C and PQ, ordinate scale right, 
show the effect of plastic flow. 


to 550°C. The curves obtained for e.m.f.-tem- 
perature difference were quite smooth, with no 
noticeable effect in the region of disordering. 
A typical result is given in Fig. 4, curve A. 
The e.m.f. values for all the crystals ranged from 
0.4 to 0.7 millivolt for a hot junction tempera- 
ture of 500°C and cold junction at 18°C. The 
variation from crystal to crystal did not correlate 
with the composition nor orientation of the 
crystals. The mean thermoelectric power in a 15° 
temperature interval was directly measured for 
six of the crystals. When plotted against tem- 
perature it shows two straight line portions 
intersecting at 460°+15°C with the higher tem- 
perature part having the smaller slope. Curve B, 
Fig. 4 is typical. Although these thermoelectric 
investigations were not very extensive they 
show conclusively that the degree of order 
affects the thermoelectric properties very slightly. 
The only effect is a slight change in the second 
derivative of the thermal e.m.f.-temperature 
relation, which occurs apparently very close to 
the temperature of completion of disorder. 


ELASTIC BEHAVIOR 


Ten of the crystals with the most uniform and 
regular cross section were selected for a pre- 
liminary study of their elastic behavior. Young's 
modulus was measured for each crystal on the 
bending beam apparatus described by Hanson.** 


26 The apparatus has been somewhat modified by Pro- 
fessor Tyndall. An improved and more sensitive optical 
system for observing the bending is now in use. Further- 














Two typical stress-strain curves are shown in 
Fig. 5, A and B. The curves were all taken with 
increasing load (circles) followed by decreasing 
load (crosses). It will be seen that there is no 
sign of hysteresis nor permanent set. In these 
curves the maximum load reading corresponds 
to a tensile stress on the top surface of the beam 
of about 20X10* dynes/cm*. The values for 
Young’s modulus for the ten crystals are given 
in Table III. 

The limit of proportionality between stress and 
strain and the beginning of plastic flow were 
obtained for one crystal. The crystal was first 
taken up to a load of 2500 g and on unloading it 
returned to zero. (See Fig. 5, curve C.) It was 
then again loaded to 2500 g, and the preceding 
readings were duplicated. An additional 500 g 
caused a rapid drift (plastic flow). The whole 
load was immediately removed,”’ and the crystal 
returned to a new zero as shown by the point 
marked P. A 1500-g load then gave the point Q. 
When this was removed point P was again 
obtained. From these results it seems clear that 
the limit of proportionality and the beginning of 
plastic flow are very close together, if not 
identical. The increase in the modulus (about 20 
percent) shown by the decreased slope of PQ is 
very marked after such a small permanent 
deformation. The tensile stress on the upper 
surface of the beam for the 2500-g load is 
70X10° dynes/cm*. The strain is 1.410~. 
A second crystal (50.05 atomic percent copper) 
was strained until the maximum tensile stress on 
the upper surface of the beam was 80X10° 
dynes/cm? without, however, exceeding the pro- 
portional limit. The “‘yield point’’ determined by 
Elam* for beta-brass is of the order of 10 times 
the above values. This difference is not surprising 
in view of the fact that Elam’s values are for 
permanent extensions sufficient to cause slip 
bands to appear. The extreme values of tensile 


more, a cylindrical roller on a flat block is used in place of 
one of the knife edges, since it was found that with two 
knife edgesa frictional drag occurred. This was undoubtedly 
the cause of the small hysteresis observed by Hanson. The 
fact that no such hysteresis exists for zinc crystals when the 
roller is used is now definitely established. Van Allen also 
obtained no hysteresis when measuring Young’s modulus 
for a zinc crystal by a direct stretch method (Thesis, State 
University of Iowa, 1936). 

*? The extra 500 g was acting at the most, for 10 seconds. 
**C, F, Elam, Proc. Roy. Soc. A153, 273 (1935-36). 
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Fic. 6. Reciprocal of Young’s modulus, 1/E, versus 
f(l, m, n) =(Pm*+m'n*?+n*P). The arrows at the bottom 
show directions in the lattice that correspond to certain 
values of f(/, m, nm). 


stress obtained by the writer would correspond 
to a resolved shear stress of about 35X10* 
dynes/cm, i.e., half the tensile stress provided a 
slip plane*® is near the most favorable orientation 
of 45° to the length of the specimen. 

For a cubic crystal Young’s modulus in a 
given direction may be expressed in terms of 
the elastic parameters, S;;, Siz and Sy and the 
direction cosines, /, m and n, as follows: 


1/E= Su - (2Si1 - 2Si2— Sua) 
X (Pm?+min?+n*P), (1)” 


where (/?m*?+m?n?+-n’/*) is called the orientation 
function, f(l,m,mn). Since 1/E is linear in 
f(l, m, n), a straight line should result when the 
reciprocals of the values obtained for Young's 
modulus are plotted against f(/, m, m). Fig. 6 is 
such a plot for the values given in Table III. 
The points, with the exception of two," lie 
within experimental error on a straight line. 
This is true in spite of the fact that these crystals 
differ widely in composition, some even being 
outside the beta-phase. The value of Si, the 
intercept on the 1/E axis (f=0), is 0.44810 
cm?/dyne. This corresponds toa Young’s modulus 
of 2.2310" dynes/cm? in the [100] direction, 
in which direction the modulus is a minimum. 
The maximum value of the modulus obtained 
from the reciprocal of the ordinate at f(/, m, n) =} 


**Slip planes are not definitely known for beta-brass. 
The glide direction is probably rit) according to Elam 
and the slip plane possibly (110) or (112). 

8 E. Schmid and W. Boas, Siruktur und Eigenschaften der 
Materie X VII, Kristallplastizitat. 

% The discrepancies for these two points are definitely 
greater than can be attributed to experimental error. The 
writer is unable to account for them. 
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js 19.910" dynes/cm*. This is for the [111] 
direction, or space diagonal. The modulus at 
f=0.25 includes among other possible directions, 
the face diagonal [110]. This value is 6.54 10" 
dynes/cm*. The material shows an unusually 
high degree of elastic anisotropy as may be 
seen from the large value, 8.94, of the ratio of 
the maximum to the’ minimum modulus in 
comparison with other metal crystals. For in- 
stance some metals of face-centered cubic 
structure have values of this ratio: Al 1.20; 
Cu 2.85; Au 2.72; Ag 2.64; alpha-brass 3.41. 
For body-centered cubic alpha-iron the ratio is 
2.15 while tungsten (also b.c.¢.) is practically 
isotropic.” Zinc, which is hexagonal, has a ratio 
of 3.6.% That there is no simple connection 
between the degree of anisotropy and the par- 
ticular type of cubic lattice, is perhaps not 
surprising because of the fact that the elastic 
behavior does not depend primarily on the 
geometry of the structure, but on the electro- 
static energy of the valence electrons and ions, 
plus an exchange energy caused by overlapping 
of ions.* 

The slope of the line of Fig. 6 is obviously, 
from Eq. (1), equal to the combination —(2S,, 
—2S\2:—S,). From this slope and the value of 
Si, the value of the quantity (2Si2.+5,,4) may 
be computed. It is 0.2810~-" cm*/dyne. The 


® Values of cubic metals from tables in E. Schmid and 
W. Boas: Kristallplastizitat, pp. 200-202. 

% This is the ratio of Ego/Eo=S,/Sz. E. P. T. Tyndall, 
Phys. Rev. 47, 398 (1935). The max./min. ratio is slightly 
larger. 

* Reference 20, Chapter IV. 


two constants 5S); and Sq cannot be found 
separately from measurements of Young's mod- 
ulus alone. 

Although this determination of the elastic 
parameters is not complete and cannot be 
regarded as of high precision chiefly because of 
the lack of sufficiently uniform dimensions of 
the crystals, it is so far as the writer knows, the 
first attempt to study the elastic behavior of 
single crystals of this alloy. In spite of the 
disabilities mentioned, the type of stress-strain 
relation is clearly shown. The small (or no) 
dependence of the elastic parameters on the 
composition is very remarkable particularly 
when it is remembered that crystals distinctly 
outside the beta-phase fall in line with the others. 
The very large anisotropy of Young’s modulus 
with direction of the lattice is a particularly 
interesting feature. It seems very likely that this 
is largely because of the high degree of order 
and that it might become considerably less above 
the critical temperature. 
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Viscosity in an Expanding Bubble 
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The rate of doing work per unit volume in a homogeneous medium is shown to be dw/di 
=V(v : ¥) where v is the velocity vector and © the stress dyadic. For an incompressible fluid 
this becomes dw/dt =2u¢ : @+(p/2)dv*/dt, where uw is the coefficient of viscosity, @ the sym- 
metric rate of strain dyadic and p the mass density. From this relation, it is shown that the 
pressure difference between the inside and outside of a steadily expanding bubble composed of 
an incompressible viscous fluid, is given by p: — p2=2y(r.-' +172") +4uK(r,-*—1r2*) where ¥ is 
the surface tension, r; and r, the inner and outer radii of the bubble, respectively, and K is the 
expansion constant r*v for any point in.the material of the bubble. 





INTRODUCTION 


BULB method has been used by one of the 
authors! for the determination of the sur- 

face tension of molten glass. Under static condi- 
tions the pressure difference between the inside 
and outside of the bubble is balanced by the 
surface tension of the liquid. When, however, a 
bubble is expanding at a uniform rate, a portion 
of the pressure difference is employed in main- 
taining the rate of expansion against the viscous 
forces then introduced. After one has determined 
the surface tension for the static condition, one 
may then determine the coefficient of viscosity. 
Experiments are under way for the determination 
of the viscosity of molten glass by this method. 
The theory has been developed for the general 
case of an expanding bubble with the assumption 
that the material is mobile to the extent that 
forces of restitution caused by what is generally 
termed elasticity may be neglected. The method 
of approach is from the standpoint of dyadics 
as employed by Page.? When we have occasion 
to refer to this text explicitly in the course of 
the paper, we will employ the abbreviation (7). 


RATE OF DoINnG WORK IN A 
Viscous FLuip 


In the notation of Page (7 P, Chap. VI) where 
W is the stress dyadic and v the velocity of the 
medium, the force on a surface element ds of the 
volume + will be (ds)- &. Since W is symmetric, 
this may be written as W&-ds and hence the rate 

1W. B. Pietenpol, Physics 7, 26 (1936). 


*L. Page, Introduction to Theoretical Physics (Van 
Nostrand, 1935), second edition. 
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at which the surface forces do work in the volume 
ris J/v-W-ds. The rate per unit volume is then 
given by 


dw/dt=\im (1 ) [ve weds= 0-9 (1) 
r=0 


or, in terms of the Gibbs double dot product 
dw/dt=Vv : ©. (2) 


As the body force in the medium is V- & per 
unit volume (7‘P, 52-1) the time rate of increase 
of kinetic energy per unit volume is 


(d/dt)(jpv*)=vV : W, (3) 
where p is the density. Hence 
dw/dt=(d/dt)(4pv?)+W : Vv. (4) 


From (7'P 76-12), 


Ww : Vv= —(p+9uV-v)V-v+2u®@ : Vv 
= — (p+ FuV-v)V-v+2u®@ :@, (5) 


since the double dot product of a symmetric with 
a skew-symmetric dyadic is zero, and ® is*the 
symmetric part of Vv. Hence (4) becomes 


dw/dt = (d/dt)(4pv*) 
—(pt+iuV-v)V-v+2lu®@ :@. (6) 


For an incompressible fluid V-v=0, and (6) 
reduces to 


dw/dt+4p(dv* dt)+2u® :® (7) 
while, if the motion is such that the rate of 
change of kinetic energy is zero, the first term 


drops out and 
dw/di=lyp®: ®. (8) 
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If, in addition, the motion of the fluid is irrota- 
tional, ®=Vv and (8) may be written 


dw/dt=2y(Vv) : (Vv). (9) 


RATE OF DoInG WoRK IN AN 
EXPANDING BUBBLE 


Consider the material of a bubble of thickness 
fy—1,, Where r; is the inner radius and rz the 
outer radius of the bubble, divided into two con- 
centric bubbles by a sphere of radius r. As the 
bubble expands or contracts let r vary in such a 
manner as to maintain the original division of 
material. Then for an incompressible fluid 


(d/dt)(r®—1,°) =3r°dr/dt—3r;dr,/di=0. (10) 


If the only motion of the bubble is that of radial 
expansion, v=dr/dt for the velocity of the fluid 
at a distance r from the center of the bubble and 


ryv=r7,;0,=K, (11) 


where K is a spatial constant. Since the vector 
velocity v at any point is directed along the 
radius vector rf, 

V=rw, (12) 


and 
Vv= (Vv)rot+vVro, ( 1 3) 


where the zero subscript is used to denote the 
unit vector. From Eq. (11) 


Vr=KV(1 r?)= —2Kry r= —2rw, ’, (14) 


and from (7P, 79-4) 


Vro= (1/1) (OOo + dogo), (15) 
hence 
Vv=(v r)( — 25 oo t+ OpOo+ dodo) (16) 
and 
Vv : Vv=60"/r*. (17) 


For a bubble expanding at a slow steady rate 
so that the rate of change of kinetic energy is 


negligible, Eqs. (9) and (17) give 
dw/dt=12yv"/r? = 124K*/r* (18) 


as the rate per unit volume for doing work 
against the viscous stresses in the fluid. Inte- 
grated throughout the material of the bubble, 
this gives 

dW /dt=16rpK*(r,* —r2"*), (19) 


for the total rate of doing work in the bubble. 


If p; is the gas pressure on the inside of the 
bubble, and #,’ the opposing pressure exerted by 
the surface tension on the inside of the bubble, 
the net pressure on the inside of the bubble is 
given by 

pi’ = pi-— fi’ = pi— 2y//n1, (20) 
where y is the surface tension. The work done 
per second by this pressure is 


dW, /di=4rp,"r,*0,;=40 Kp, —82Ky/ri. (21) 


Similarly, the work done per second by the net 
pressure on the outside of the bubble is given by 


dW, dt= —4nK p.—8xKy/re. (22) 


But since the work done by the net pressures 
must equal the work done against the internal 
stresses, Eqs. (19), (21) and (22) may be com- 
bined to give 


Pi— P2=2y(ri ttre) +4uK (ri —re*). (23) 


If the thickness d of the bubble is small compared 
to the average radius, Eq. (23) may be replaced 
by the approximate relation 


pi— po=4y/r+l2pod/r’, (24) 


where r is the average radius, and v the average 
radial velocity of the bubble at any given instant 
and ~:— p2 is the pressure difference required to 
balance the viscous forces in a steadily expanding 
bubble of finite thickness. This formula is par- 
ticularly applicable to bubbles of extremely 
viscous materials such as molten glass. 




















FEBRUARY 1, 1939 


PHYSICAL REVIEW 











VOLUME 55§ 


Actinoelectric Effects in Tartaric Acid Crystals* 


James J. Brapy AND WiLt1AM H. Moore 
Oregon State College, Corvallis, Oregon 
(Received October 26, 1938) 


Tartaric acid crystals illuminated with light from a carbon arc exhibit a flow of current when 
connected to a galvanometer without the aid of an impressed battery. This actinoelectric effect 
is greatly dependent upon the orientation of the crystal with respect to the illuminating beam 
of light. It is shown that a section of the crystal can be chosen which shows a series of equi- 
potential lines when its surface is examined with two probes while the crystal is illuminated. 
The measured currents are about 10~" ampere, when the crystal is illuminated with the light 
from a carbon arc. The effect is a linear function of light intensity for the lower light intensities, 
but increases more rapidly for higher intensities. The effect is practically independent of 
temperature in the range from 20°C to —40°C. The maximum spectral response is at approxi- 


mately 10,500A. 





NUMBER of experimenters' have shown 

that certain crystals, when illuminated, 
produce an electrical current which can be 
measured with an electrometer or galvanometer 
connected in series with the crystal without the 
aid of an impressed battery. According to 
Coblentz? the term actinoelectricity was used by 
Hankel to designate the above phenomenon 
about 95 years ago. Present theories' are not 
entirely satisfactory, and it appears desirable to 
investigate the effect in more detail. 

Most of the crystals used in this investigation 
were chosen from a supply of commercially pre- 
pared tartaric acid crystals.* The results were 
checked, however, with crystals grown in our 
laboratory. The ‘dark resistance’”’ of a crystal 
three mm thick with a cross-sectional area of 
approximately 50 sq. mm is 210° ohms. 

The experimental arrangement for carrying out 
the investigation is shown in Fig. 1. No results 
are reported when an external source of e.m.f. 
is used. A carbon arc and a 150-watt gas-filled 
tungsten filament lamp are used as light sources. 

An FP-54 tube is used in a Brown and Du- 
Bridge‘ circuit for amplifying the actinoelectric 





* Published with the approval of the Monographs Pub- 
lication Committee, on State College. Research Paper 
No. 13, Department of Physics, School of Science. 

! For bibliography and a review of theories see B. Lange, 
Photoelements (Reinhold Publishing Corporation, 1938). 

*W. W. Coblentz, Scientific Papers of Nat. Bur. of 
Standards, No. 486, April (1924). 

* Maximum limits of impurities: SO;—0.01 percent; 
Ca—0,001 percent ; Fe—0.001 percent ; other H. M.—0.000 


percent. 
*H. Brown and L. A. DuBridge, Rev. Sci. Inst. 4, 532 
(1933), 
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currents. The current sensitivity of the amplifier 
is 10-“ ampere per mm deflection of the gal- 
vanometer. 

The crystal is held between two brass disks 
with spring tension on one of the disks. The 
crystal holder is designed so that the light can 
be sent either perpendicular or parallel to the 
direction of the current through the crystal. 

Contradictory evidence has been reported by 
other experimenters as to the origin of the effect. 
Coblentz,? working with molybdenite shows that 
the seat of the effect is at certain points along 
the crystal surfaces. Geiger’ finds, in the case of 
argentite, no observable effect when the contacts 
are shielded from light and that the currents are 
measurable only when at least one of the elec- 
trodes is illuminated. Schneider*® finds, in the 
case of argentite, that the seat of conversion of 
light into electronic energy occurs both at the 
contacts as well as at certain spots on the 
crystal. 





TTT FP-S+ 
Crystal/] } 














= 





—— 


Fic. 1. Experimental set-up. 


5 P. H. Geiger, Phys. Rev. 22, 461 (1923). 


®*W. A. Schneider, Phys. Rev. 31, 82 (1928). 
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The experiments with tartaric acid crystals 
show none of the irregular variations in current 
as the light beam is moved over the surface of 
the crystal as reported for instance by Coblentz.? 
There were no “sensitive spots’ found in the 
case of a good clear single crystal although the 
response in many cases varied in a regular manner 
as the surface was explored with a light beam of 
approximately one mm diameter. 

Figure 2 represents the way in which the 
current varies with the time of illumination when 
the light beam is at right angles to the flow of 
the current. The maximum deflection is reached 
within the period of the galvanometer and then 
drops down to a final value. This final current 
remains constant even though the crystal is 
illuminated for one-half hour or more. As soon 
as the light is cut off, the current immediately 
reverses and then drops to zero in the manner 
indicated by the curve. The decrease in current 
which takes place a few seconds after the light 
is turned on is obviously caused by the building 
up of a space charge or polarization within the 
crystal which limits the flow of electrons. When 


Light off 





Deflections 








Galvanometer 
> 
qT 








-6 
Time in seconds 


Fic. 2. Variation of actinoelectric current with time 
of illumination. 


the light is turned off, the space charge is 
neutralized by the rush of electrons in the 
opposite direction. The residual current, which 
remains constant after the light has been shining 
on the crystal for some time, indicates that 
electrons enter and leave the crystal through the 
electrodes. 

Crystals which appear to have flaws and are 
not very transparent give results which are 
notably different from those from good clear 


crystals. The current, in these crystals, starts 
out in one direction just after the light is turned 
on, quickly reverses, and then slowly reaches a 
maximum (approximately after one minute) in 
the opposite direction. 

Figure 3 represents a single crystal of tartaric 
acid. When the light beam is parallel to the 
direction of the current (longitudinal effect), the 
following results are recorded. With the c axis 
vertical, the electrodes against the orthopinacoid 
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Fic. 3. A monoclinic crystal of tartaric acid with the 
crystal faces labeled —— to crystallographic notation. 
a q—clinodome; c—basalpinacoid; r— 
orthodome. 


faces, and the clinodome face to the right looking 
at the illuminated orthopinacoid face, the elec- 
trons in the crystal flow in the same direction as 
the light energy. In this case, a residual current 
remains after long periods of illumination. When 
the crystal is turned through 180° about the c 
axis so that the opposite orthopinacoid face is 
turned toward the light, the electron flow is still 
in the direction of light propagation. When the 
crystal is turned from this position through 180° 
about the 6 axis, the electron flow is still in the 
direction of light propagation, but no residual 
current remains after 30 sec. or more of illumina- 
tion. If the crystal is now turned through 180° 
about the ¢ axis, the electrons flow in a direction 
opposite to the direction of light propagation. 
In this case, no residual current remains. These 
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Fic. 4. A section of the crystal cut parallel to the 
clinodome face. 


results are surprising because the crystal, to all 
outward appearances, is perfectly symmetrical 
about the 6 axis. 

With the two parallel electrodes against the 
orthopinacoid faces, a light beam falling on the 
upper clinodome face (transverse effect) produces 
a deflection of the galvanometer which is just 
equal to but of opposite sign from the response 
of the lower clinodome face. The electrons in the 
crystal travel in opposite directions as the light 
is shifted from one clinodome face to the other. 

A very interesting study is made by cutting a 
section of the crystal parallel to the clinodome 
face and exploring its surface with two probes 
while it is illuminated from above the surface 
with a carbon arc lamp. A sheet of Cellophane is 
marked off in two-mm squares and holes made at 
the corners of the squares to receive the probes. 
The Cellophane is then cemented to the crystal, 
and the probes are used to measure the current 
between two points on the crystal surface as the 
probes are inserted into the holes in the Cello- 
phane. No current is observed when the two 
probes are placed anywhere on the line AB 
(Fig. 4). The same is true for all lines parallel 
to AB. These lines are very nearly parallel to 
the side CD. Appreciable galvanometer deflec- 
tions were noted when one probe was placed at 
a point on the line AB and the other probe was 
placed at a point on one of the other equi- 
potential lines. 


LiGHt INTENSITY INVESTIGATION 


The way in which the current varies as a 
function of light intensity is shown in Fig. 5. 
The initial value of the current (Fig. 2) is used. 
The experimental arrangement shown in Fig. 6 
is used to obtain this curve. The light intensity 
is varied by changing the number of glass plates 
interposed between the carbon arc and the 


AND 





W. H. MOORE 
crystal. The light reflected from the glass plate A 
falls on a light shielded photronic cell which is 
connected to a galvanometer. The light intensity 
falling on the photronic cell is in a region where 
the response of the cell is linear. The light 
intensity falling on the crystal is then directly 
proportional to the deflections of the photronic 
cell galvanometer. The curve is similar to that 
obtained by Schneider* for the actinoelectric 
effects in argentite. 
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Fic. 5. The variation of the actinoelectric current with 
light intensity. 


From the data given by Barnes and Forsythe’ 
for the spectral distribution in energy from a 
100-watt gas-filled tungsten lamp, the total in- 
tensity falling on a sq. cm at one m distance is 
determined. In this way the thermocouple is 
calibrated in microwatts/cm* of radiant energy 
falling on the junction. After the thermocouple 
is calibrated in microwatts/cm*, it is used to 
determine the intensity of light falling on the 
crystal. This value is compared with the response 
from the photronic cell. The output from the 
photronic cell is then used to measure the light 
intensities in microwatts/cm*. 


EFFECT OF TEMPERATURE 


In order to study the effect of temperature on 
the current response, a crystal is mounted in an 
evacuated Pyrex tube to prevent moisture from 
collecting on the crystal as its temperature is 
lowered. The glass tube is covered with tinfoil 
except for the portion near the electrodes and is 


7 B. T. Barnes and W. E. Forsythe, J. Opt. Soc. Am. 26, 
313 (1936). 
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then surrounded with solid carbon dioxide. The 
results show no effects of temperature greater 
than experimental error in the range from —40 
to +30°C. 

This is in agreement with the findings of 
Coblentz® for molybdenite. 


SPECTRAL RESPONSE 


Figure 7 shows the way in which the actino- 
electric current varies with the wave-length of 
the light. A carbon arc is used as a light source 
and is dispersed by means of a double-prism 
monochromator. The arc has a maximum in- 
tensity at approximately 8500A. The energy 
distribution of the light from the monochromator 
is determined by means of a_ thermocouple 
designed for this purpose. The actinoelectric 
currents are then corrected for the energy distri- 


2 Light source 
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Pile of glass plates 
y Cardboard tube for light shield 
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Fic. 6. The experimental set-up for measuring the variation 
of current with light intensity. 
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Fic. 7. Spectral response curve. 


bution, and the currents are plotted on an equal 
energy scale. Because of the small galvanometer 
deflections the monochromator slits must be 
widened so the wave-length band has a width of 
300A at 5500A. 

The maximum response can be seen to be at 
approximately 10,500A. The light is in a direc- 
tion perpendicular to the electron flow through 
the crystal in this case. The actinoelectric re- 
sponse is linear with light intensity for the 
intensities emerging from the monochromator. 

In conclusion, the authors wish to express their 
thanks to Dr. W. D. Wilkinson of the Geology 
Department of Oregon State College for orienting 
the crystals. The senior author wishes to express 
his appreciation to the Physics Department of 
the University of California for the loan of the 
monochromator and to the National Academy of 
Sciences for the grant-in-aid which made possible 
the purchase of some of the apparatus used in 
this research. 
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Shape of the Domains in Ferromagnetics 


E. H. Kennarp 
Cornell University, Ithaca, New York 


(Received November 28, 1938) 


T is generally agreed that a demagnetized 

ferromagnetic crystal consists of domains each 
magnetized uniformly to saturation, the appa- 
rent demagnetization being the result of varia- 
tion in these directions from one domain to 
another. The domains are estimated to contain 
10" to 10" atoms each, but there is no agreement 
as to their shape. The common assumption that 
they are roughly cubical can scarcely be justified. 
It seems worth while to restate the objection to 
this assumption in a more concrete form; then the 
condition which must be satisfied by the equi- 
librium positions of domain boundaries in an 
ideal cubical crystal will be thrown into a new 
form. 

In a demagnetized iron crystal one-sixth of 
the domains are commonly assumed to be mag- 
netized to saturation in each of the directions 
parallel to the three axes. Now if such a domain 
were more or less round or cubical in form, it 
would necessarily be surrounded by a local field 
of the order of Jo, the intensity of magnetization 
within a domain; this becomes evident upon re- 
calling the formulas for a sphere. This local field 
might happen to be offset at some points by 
fields due to other sources, but the cancellation 
cannot be universal; for over the part of the 
surface where one domain abuts upon another 
in which the direction of magnetization differs 
by a right angle or more, there is a sheet of 
poles of density roughly equal to J», and hence 
a jump in the normal component of H roughly 
equal to 4xJ». In an iron crystal, with J>= 1800, 
the resulting local fields would be of the order 
of 2000 oersteds. Yet it is observed that an 
applied field of less than 10 oersteds suffices to 
magnetize the crystal as a whole to saturation. 
The conclusion seems inescapable that cubical 
domains could not be stable; they would tend to 
alter each other’s magnetization so as to become 
lined up into rows forming long slender domains 
of larger size. Thus we reach by elementary 
reasoning the old conclusion that the domains in 
iron, and presumably in all ferromagnetic crystals, 
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must have the form either of slender filaments or of 
thin plates, magnetized longitudinally.' The plate 
form was suggested for cobalt by Landau and 
Lifschitz.? 

The domain boundaries are assumed to be so 
located as to minimize the total energy due to 
the Weisz-Heisenberg or quantum field and to 
magnetic interaction. Considering only cubic 
crystals, let us assume all fields to be weak 
enough so that each domain is magnetized sen- 
sibly along one axis. Then the quantum energy 
varies only because there is an excess of it in 
the domain boundaries, proportional to their 
area. Presumably there will be a transition layer, 
through which the direction of the local mag- 
netization varies rapidly from that of one axis 
to that of another; additional excess energy will 
be present here due to the occurrence of spins 
lying oblique to the directions of easiest mag- 
netization in the crystal. These excess energies 
together are equivalent to a surface tension T 
tending to contract the boundaries. 

The effect of the magnetic energy can some- 
times be seen directly from the usual expression, 


U= — ft-Hedr+ fe 8x)dr, 


where I is the magnetization and H, the field 
due to it (both macroscopic, but on the scale of 
the domains), and Hy is the field due to all other 
sources. In the absence of Hy the magnetic forces 
will cause the domains to shift so as to obliterate 
H;. But in many situations a criterion involving 
only conditions at a single point on the boundary 
is more useful. 

To obtain this, suppose that unit area of a 
boundary, with principal curvatures C’ and C”, 
shifts a distance 6s normal to itself. Then by 
surface-tension theory the boundary energy is 


1 Cf. R. H. de Waard, Phil. Mag. 4, 641 (1927); F. Bloch, 
Zeits. {. Physik 74, 295 (1932). 

* L. Landau and E. Lifschitz, Physik Zeits. Sowjetunion 
8, 153 (1935). 











anaoewe  elUcrcrcrlCiteelC rr lC CO 





or of 
dlate 
and 


SO 
e to 
d to 
ubic 
yeak 
sen- 
rgy 
t in 
heir 
yer, 
lag- 
axis 
will 
ins 


ries 


\z 


ne- 
on, 


old 


ler 


ite 











increased by (C’+(C")Tés. The increase in mag- 
netic energy is 


6U= — f Hy: star— 36 [T-Bar = — f H-atar, 


where H=H,+H,, and the last expression re- 
sults from the following transformations, in 
which we introduce the potential due to I, so 
that H,=—VoO, VQ= —div. H,=42 div. I, and 
then carry out integrations by parts as in the 
divergence theorem, the integrated terms vanish- 
ing always at infinity: 


rata: - ~ ft-viadr= f a0 div. Idr 


= f savendr, tr= fo V*5Qdr /4r 


- {2 div. ildr= f H-dldr. 


Now at any point in the transition layer 6/, 
=—TI,’is, 51,,= —I,,'6s; subscripts denote, re- 
spectively, the perpendicular component in the 
direction of 5s and the vector parallel component, 
and primes indicate the space derivative. Hence 


6U= is f (Ht ‘ln’ +Hil,")dz 


integrated through the transition layer. But 
H,,'=0, and H/,’= —4xI,’ by Gauss’ law; thus, 
if subscripts 2 and 1 indicate values on the two 
sides of the boundary, 


feet z=Hi- tu'ds= hu (nan): 


jf ina- ~ f atta: 4dr 


= — (F742? — H,,;°) 8a = 4( Hist Hay) (La2—T). 


Hence we can write, since H,, =(H,,2+H,,:)/2, 
6U=3(H, +H) - (I,—1,)és. 


The condition for equilibrium is now that the 
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total change in energy must vanish, which re- 
quires that 


}(H,+H,) -1, = }(H, +H) -bh+(C'+C")T. 


If the quantity on the left in this equation 
exceeds that on the right, the boundary shifts 
away from medium No. 1, and vice versa. (C’ 
and C” are positive when they mean convexity 
toward domain No. 2.) A plane boundary, 
therefore, shifts away from the side on which 
3(H,+H,)-I is the greater, and is in equilibrium 
when 


}(H, +H) -1, = $(Hi+H)) -h. 


Here }(H,+H:) represents the mean field, due 
to all other sources than free poles on the part of 
the boundary under consideration; for equilib- 
rium of a plane boundary this mean field must 
be equally inclined to I in the two adjacent 
domains. 

The magnitude of the surface tension is not 
definitely known. According to Bloch* the tran- 
sition layer between domains represents a com- 
promise between the tendency of the quantum 
forces to decrease the gradient of the local mag- 
netization and the tendency of the crystal forces 
to turn all spins into the direction of a crystal 
axis, and T should be of the order of the work 
required to turn them halfway from one axis to 








another in a layer about 30 crystal spacings 
thick. In iron this would make T=2X10° X30 
xX3X10-*=0.18 erg/cm’, roughly. If the surface 
tension is really of this order of magnitude, the 
shape of the domain boundaries must be deter- 
mined largely by the local magnetic fields. 

The rules obtained above are useful in making 
qualitative studies. It is readily verified that 
cubical domains would tend, because of magnetic 


*F. Bloch, Handbuch der Radiologie, Vol. 6, Part 2 (1934), 
p. 481. 
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interaction, to change as suggested above. Equi- 
librium in an ideal crystal seems, moreover, to 
require that the domains extend to the surface; 
for any domain ending in the interior tends to 
shrivel up, both because of the surface tension 
and because of magnetic action. The field due 
to the domain itself tends to collapse its sides, 
and the field due to the poles at one end tends 
to pull in the other end. In an actual crystal, 
however, there appears to exist a certain re- 
sistance to the motion of domain boundaries, so 
that it might be sufficient for a domain to 
minimize its own magnetic field by tapering very 
gradually to a point, or to a sharp edge.‘ 

At the surface, again, the magnetization will 
adjust itself so as to keep the local fields within 
bounds. A possible arrangement below a surface 
perpendicular to one crystal axis seems to be as 
suggested in Fig. 1. If the surface tension were 
absent, the magnetization, shown by arrows, 
would be carried round through triangular prisms 
lying parallel to the surface, with their oblique 
faces at 45° to the axis (as suggested for cobalt 
by Landau and Lifschitz’), and there would be 
no free poles at all. In reality, however, the 
boundaries must be pulled up a little by their 
surface tension, with the production of free poles 
as shown. The local fields due to such poles 
might be the cause of the powder patterns that 
are often observed; a magnetic colloid would 
collect in the regions of strongest field at ccc: --. 
If an external field were then applied in the up- 
ward direction, the strongest field might be 
shifted to aa---, and the powder ridges would 


*K. J. Sixtus and L, Tonks, Phys. Rev. 37, 930 (1931). 
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occur there; whereas in a reverse field they would 
be at 6b---. Thus the ridges in the presence of a 
field would be in locations halfway between the 
no-field positions and twice as far apart; this is 
just what is observed.° 

The application of the field would also change 
the domains somewhat, of course. In the figure 
an upward field would broaden those under aa- - - 
and should probably tend to make their tops 
dome-shaped, with an increase in the pole 
strength; whereas the tops of the domains under 
bb--- would become more peaked, with a de- 
crease in the pole strength there, or even a 
reversal of its sign. These changes would intensify 
the phenomena just described, at least so long 
as the field did not become too strong. 

Even if the effect of surface tension is neg- 
ligible, an imperfection in the triangular prisms 
might result from the resistance which is ob- 
served to impede the motion of domain bound- 
aries. It is not so easy, however, to understand 
how the plate-like domains can turn square cor- 
ners as freely as they must do to fit the powder 
patterns. Furthermore, isotropic demagnetiza- 
tion is hard to imagine, for spindles and plates 
lying in all three directions can scarcely fit to- 
gether without gaps. In a demagnetized crystal, 
do the shortest dimensions of all domains perhaps 
lie in the same plane? If so, the observed mag- 
netization curve should differ according as the 
crystal is demagnetized by use of a parallel or 
a perpendicular field. This point might be worth 


testing experimentally. 


5’ W. C. Elmore, Phys. Rev. 51, 982 (1937); 53, 757 
(1938). 
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Communications should not in general exceed 600 words in length. 


On the Interchange Energy of Two Free Electrons 


In this note, a comparison is made of the Coulomb and 
interchange energies for a system consisting of two free 
electrons moving arbitrarily with respect to each other. 

The following equation given by Dirac! represents the 
perturbing energy for a system of similar particles, and 
contains the expressions for Coulomb and interchange 
energies. 

V=Vi-}3 > Val it(e,, ¢)}, (1) 
< 


™ 


where V is the perturbing energy, ¢ is the spin dynamical 
variable, and the subscript rs refers to an interchange of 
the r and s electrons. For any permutation, P, the following 
equation is given: 


Vp=(Pa| V\a)= Sf Sf (Pa\ g’)dq'(q’| V\q"’)dg’’(g”’ |) 


where a refers to the states and g to the coordinates of the 
electrons. 

In order to evaluate the integrals involved in Eq. (1), 
the wave functions are replaced by wave packets repre- 
senting free electrons. If we choose the coordinate system 
so that one electron is at the origin, and the other is at 
the point ro on the z axis; and assume that the uncertainty 
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FRACTION OF CLASSICAL COULOMB ENERGY 
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Fic. 1. Energy relations for two electrons at rest relatively to each 
other. C=Coulomb energy classical Coulomb energy = V1:/ (¢*/re); 
] =Interchange energy ‘classical Coulomb energy = V12/(e/re); C —/ 
=Total perturbing energy/classical Coulomb energy = V /(¢*/re). 


is the same for both electrons, and the same in all direc- 
tions, we find, after integrating : 


Vii = (/re) erf (ro/24e), (2) 
Vis = (@h/wlotmn,) exp [(—re/2e%) —2*] if exp (#)dt, (3) 


V= Viu- Vie, 


where V;, is the Coulomb energy, V1. is the interchange 
energy, fo is the distance between the electrons, %% is the 
relative velocity between the electrons, \ = 24mrvee/h; and 
¢ is the uncertainty. 

Having evaluated V;, and V2, we may now determine 
V for several interesting cases. 

When the positions of both electrons are known exactly, 
i.e., when the uncertainty, ¢, is equal to zero, V equals 
e?/ro. Thus in this case the interchange energy is zero and 
the total perturbing energy is simply the Coulomb energy. 

When the distance, fo, between the electrons is large, 
the interchange energy is negligible compared to the 
Coulomb energy, i.e., Vi2/ Vin approaches zero, as ro ap- 
proaches infinity. Thus the interchange energy diminishes 
as the relative velocity between the electrons increases 
From this it is seen, as would be expected, that the inter- 
change energy is dependent upon the time during which 
the electrons are in the neighborhood of each other. 

Finally, since Vi_ is a maximum when % is zero (re and 
o remaining fixed), it is interesting to calculate V for this 
case. It is found by making A--0 in Eq. (3), that 


V = (e*/re) erf (ro/24e) — (e*24/ert) exp (—1r¢/2e*). 


Figure 1 shows the variation of Vi;, V;; and V, for various 
values of ro/e. 

From curve | it is seen that the interchange energy is of 
considerable importance when the uncertainty of position 
of the electrons is of the order of magnitude of the distance 
between them. However, as the wave packets, representing 
the electrons, become separated, the Coulomb energy 
approaches the classical Coulomb energy, and the inter- 
change energy approaches zero. 

I wish to thank Dr. Boris Podolsky, University of 
Cincinnati, for suggesting this problem, and for many 
helpful discussions. 

Cart A. Lupexe* 

Graduate School, 

University of Cincinnati, 
Cincinnati, Ohio, 
November 30, 1938. 


* Now at John Carroll University, Cleveland, Ohio. 
' P. A.M. Dirac, Quantum Mechanics, second edition, p. 228, Eq. (38). 
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Cosmic Rays and Radioactive Potassium 


It is probably well to put on record the results of an 
experiment completed last spring. As A. K. Brewer has 
pointed out,’ the present abundance of radioactive K* 
and its known life indicate that practically all the Ca and 
A now existing on the earth must have been produced by its 
disintegration. This conclusion seems unavoidable unless 
the supply of K* has been replenished in some way. 
Several nuclear physicists stated to the writer their 
belief that K® is formed from K** by the absorption of a 
neutron, one source of which would be cosmic rays. For 
such an hypothesis the present intensity of cosmic rays 
seems inadequate, but two facts give it a certain plausi- 
bility, namely, that radioactivity seems to be confined to 
the outer crust of the earth and that, although there is 
no evidence of large differences in the radioactivity of 
potassium samples, it is true that most of those tested are 
from sources which have been near the surface of the earth 
during a large part of geologic time. It was suggested to 
the writer that we investigate this idea since the necessary 
equipment for testing small amounts of potassium was 
available here. 

For this purpose Dr. H. J. Fraser of the geology depart- 
ment procured for us samples from an old geological forma- 
tion in Canada in which a continuous mass of granite ex- 
tends from the surface to great depths. The rock from 
which these samples was taken had not been subjected to 
the action of water and we can state with certainty that, 
during geologic time, it has never been closer to the surface 
than it now lies. A piece of granite from the greatest 
depth, 2500 feet, was given to Mr. D. K. Beavon, who 
under the direction of Dr. E. H. Swift of the Chemistry 
Department, prepared a 22-mg sample of KCI from it. 

The counter on which the measurements were made 
was the one used for determining the radioactivity of the 
separated isotopes of rubidium.” * It is a cylindrical 
counter one cm in diameter and three cm long with an 
aluminum wall 0.00125 cm thick. The sample was fused 
as uniformly as possible on the inner surface of a half- 
cylinder of platinum foil of one cm radius and two cm 
length, which, for counting, was swung into position 
concentric with the counter so that all parts of the sample 
were equidistant from its center. The bare foil when in 
place gave the standard background count of one per 
minute, which 20 mg of ordinary KCI should increase to 
about two per minute. It was calculated that with such 
thin samples errors due to absorption with uneven distri- 
bution would be negligible. Redistribution of the samples 
confirmed this conclusion by leaving the count unchanged. 
Similar sized samples of chemically pure KCI and of KCl 
specially prepared from K,SO, to avoid radioactive con- 
tamination, were mounted in the same way. Alternate 
fifteen-minute counts were made on the sample and a 
standard screen, until at least 400 counts on each were 

recorded. Such a set constitutes a “run.” Repeated runs 
were made on all samples, sometimes with a moistening 
with HCl and redistribution between runs. The C.P. and 
the specially prepared specimens gave nearly identical 
counts, but the Canadian sample ran consistently lower. 
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The two final samples were turned over to Mr. J. S. Bill- 
heimer for analysis. He found the comparison sample to 
be pure KCI but the Canadian sample contained about 
20 percent of impurity, probably sodium. When correction 
was made for the impurity, the Canadian sample was 
found to be about ten percent less radioactive than the 
laboratory samples. Since the probable error of the experi- 
ment is about ten percent, we conclude that the amount of 
radioactive potassium produced by cosmic rays is unim- 
portant and probably zero. 


W. R. Smyrue 


California Institute of Technology, 
asadena, California, 
December 23, 1938. 


'A. K. Brewer, Ind. and Eng. Chem. 30, 893 (1938). 
?W. R. Smythe and A. Hemmendinger, Phys. Rev. 51, 178 (1937). 
+A. Hemmendinger and W. R. Smythe, Phys. Rev. 51, 1052 (1937), 





Photoconductivity of Metal Films 


R. S. Bartlett'.* in 1925, while working at this labora- 
tory, showed that a photoconductive effect could be 
observed in thin films of copper, platinum, palladium and 
bismuth sputtered on glass. On exposure to ultraviolet 
light at liquid-air temperatures, decreases in resistance up 
to sixteen parts in one million were recorded. At room 
temperature, and with frequencies of incident light lower 
than the ultraviolet range, no effect was noticeable. 

To establish whether this effect could be detected in 
conductors of bulk metal, a further experiment has recently 
been completed on fine ribbons of copper, platinum and 
palladium. The illumination was with x-rays in order to 
insure that appreciable absorption took place throughout 
the entire thickness of the conductor. 

Ribbons approximately 0.01 mm in thickness, 1 mm 
wide and from 2 to 3 meters in length were rolled from 
fine wires of the materials mentioned, and wound upon 
frames so that the majority of the surface could be exposed 
to an x-ray beam (molybdenum target). At room tempera- 
ture the resistance of the element was measured by a 
potentiometer method. No change in the resistance of 
amount greater than one part in one million could be 
detected at this temperature when the element was 
illuminated. At liquid-air temperature, fluctuations of the 
resistance of the element caused by temperature changes 
and various parasitic e.m.f.’s made the potentiometer 
method of detecting resistance changes unsuitable. A 
shutter was therefore interposed in the x-ray beam and 
rotated to interrupt the beam 1000 times a second. A 
steady direct current was passed through the element and 
the e.m.f. across its terminals fed to an amplifier-oscillo- 
graph circuit. No alternating component of frequency of 
1000 cycles could be detected on the oscillograph with 
any of the three metals used. Although the accuracy of 
the alternating current method was estimated to be only 
one part in 250,000, because the elements were at most 
only a few ohms in resistance and the available e.m.f.'s 
across their terminals limited to a few volts, the method 
was still sensitive enough to have detected changes in 
resistance of the order of magnitude of those observed 
by Bartlett. 
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The conclusion therefore may be made that either 
(1) the photoconductive effect present in thin films of 
these metals sputtered on glass surfaces is a property of 
the bulk metal which can be caused by radiation of the 
ultraviolet range only, or (2) that the effect is not a 
property of bulk metal but is restricted to films alone 
and is related to their structure,—the effect may be an 
external photoelectric current passing between discon- 
tinuous patches of the film. 

The latter alternative would seem to be supported 
experimentally by recent work by Lovell,* and Lovell and 
Appleyard‘ on the resistance of thin films of the alkali 
metals. 


Tuomas C. WiLson 


Sloane Physics Laboratory, 
Yale University, 
New Haven, Connecticut, 
January 6, 1939. 


'R. S. Bartlett, iy Rev. 26, 247 (1925). 

2A. T. Waterman, Phys. Rev. 22, 259 (1923). 

+A. C. B. Lovell, Proc. Roy. » on. A157, 311 (1936). 
a" Ne Lovell and E. Appleyard, . Roy. Soc. A158, 
18 (19 





On the Upper Limit of the 8-Spectrum of ThC” and ThB 


In 1937 Lyman! published a paper in which, from the 
character of the 8-spectrum near the upper limit, he points 
to the existence of a neutrino rest mass different from zero. 
Subsequently a number of investigators, such as Alicha- 
nian, Alichanow and Dzelepow* and also Alichanian and 
Nikitin,* undertook a thorough study of the upper limit of 
the 8-spectra. 

According to the last authors, the difference between 
the experimental value for the upper limits and that 
deduced from the theory of Konopinski and Uhlenbeck* 
determines the magnitude yo of the mass of the neutrino, i.e. 


Wo— (Eo/511+1) =no, 


where W, is the extrapolated value of the limiting energy 
(in mc*) while E» is the experimental value of this energy 
in kv. 

If so, it has to be assumed that the 8-decay of different 
elements is characterized by different values of yo or else 
to associate each elementary process of 8-decay with an odd 
number of neutrinos, as has been suggested by G. Breit.* 

However, in our opinion both alternatives are equally 
unsatisfactory. The comparison of the shape of the experi- 
mental §-spectra with that given by Konopinski and 
Uhlenbeck (KU) has shown that the substitution of the 
value W, (obtained by extrapolating the KU plots) in 
the expression for the distribution of electrons according 
to their energies in the 8-spectra does not lead to a closer 
agreement of the experimental 8-spectrum with the calcu- 
lated one. More satisfactory results can be obtained by 
taking Wo.’ between the limits 

(1+ Eo/511) << Woe’ < We. 

The computations with the experimental results ob- 
tained by Alichanian and Zavelsky* for ThC” and ThB 
are shown in Figs. 1 and 2. 


The solid curves refer to the experimental 8-spectra for 
ThC” and ThB. They are in fairly good agreement with 
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Fic. 2. Konopinski-Uhlenbeck plots and experimental results for ThB. 








those obtained when substituting W,’ in the KU theoretical 
spectrum. The dotted lines show the results obtained by 
substituting W». The existence of the above disagreement 
has been already pointed out by the author’ for ThB, RaE 
and ThC”. 

The following conclusions may be drawn: 

1. The equation given by Konopinski and Uhlenbeck 
must be considered as an empirical equation [/(s,W) 
(Wo’—W)*~N)], requiring a suitable choice of the value 
W,’, distinct from that obtained by extrapolating the 
KU plots. 

2. The assumption that the difference Wy— (E»/511+1) 
determines the rest mass of the neutrino, as expressed by 
certain authors, is not sufficiently justified. 


A. ZAVELSKY 


L. Tenteleva 24, Apt. 9 
Leningrad 95, U.S.S.R., 
December 14, 1938. 


x M. Lyman, 4g Rev. 51, 1 (1937). 
I. Alichanian, A. I. Alichanow and B Deelepow, Phys. Rev. 
33, 766 ( 1938). 
. 1. “Alichanian and S. J. Nikitin, Phys. Rev. 53, 767 (1938). 
Konopinski and G. Uhlenbeck, Phys. Rev. 48, 7 (1935). 
.G. jreit, a Sci. Inst. 9, 3, 63 (1938). 
*A. I. Alichanian and A. Zavelsky, C. R. l"USSR, 17, 467, 469 (1937). 
7A. Zavelsky. Information at the 3rd All-Union Symposium on 
Nuclear Problems, in October, 1938 (Russian). 
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An Electrical Quadrupole Moment of the Deuteron* 


The molecular beam magnetic resonance method! ap- 
plied to HD molecules at liquid nitrogen temperature 
gives the magnetic moment of the proton and the deuteron. 
When applied to H; and D, molecules the method reveals 
the close groups of sharp resonance minima shown in Figs. 1 
and 2. The resonance minima for H,; agree in number 
and location with predictions based on the assumption 
of spin-spin magnetic interaction of the two protons 
(1 «2/7? — 3 (1-7) (u2-7)/r*) and a spin-orbit interaction of 
the protons with the rotation of the molecule (2upHMJ-J). 
All symbols here have their usual significance except A 
which is the spin-orbit interaction constant. The only state 
of the molecule to be considered is the lowest rotational 
state of orthohydrogen: J=1 and total nuclear spin 7=1. 
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The nine energy levels which arise from these interactions 
and from the external magnetic field give six possible 
transitions for the nuclear spin because AMy=-+1. The 
pattern can be accounted for completely on the basis of the 
known value of the proton moment (up=2.78) with the 
known value of the internuclear distance and a value for M 
of 27 gauss. 

In the case of D, the deep central minimum arises from 
the states with J=2 and J=0. The six smaller peaks arise 
from the states with /=1 and J=1. The states with J/=2 
are not abundant enough at these low temperatures for 
observation. Since the internuclear distance in the D, 
molecule is the same as in H,; and the mass is twice as great 
the spin-orbit interaction constant AT is half as great. 
The deuteron magnetic moment (up =0.853) is 0.307 times 
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Fic. 1. Resonance minima for H in H2. 
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Fic. 2. Resonance minima for D in Ds. 
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that of the proton and therefore the magnetic spin-spin 
interaction is proportionately smaller. It was expected that 
the theory of the resonance minima for H: when applied 
to D, should give the locations of the minima from the 
constants given above. The displacements of the minima 
from the center should be much less than those oj Hs. 
However experiment shows the displacements to be six 
times greater than the predicted values. 

This effect can be accounted for by the presence of an 
electrical quadrupole moment in the deuteron. The inter- 
action energy which gives rise to the large displacements 
is that of the nuclear quadrupole moment with the gradient 
of the molecular electric field. This form of interaction 
contributes to the nine energy levels exactly as the spin- 
spin and therefore appears as a larger spin-spin interaction. 

To prove that the large displacements in D, are of 
nuclear origin rather than molecular, similar experiments 
were performed on the proton and the deuteron in the HD 
molecule. The group of resonance minima for H was 
narrow as expected and that for D had large displacements 
as in D,. Furthermore, the experimentally evaluated spin- 
orbit interaction constant for D, is one-half as great as 
that for Hz as predicted. We therefore believe that the 
apparent large spin-spin interaction is not magnetic, nor 
is it of molecular origin and must be a nuclear effect which 
behaves like a quadrupole moment. 

To obtain the magnitude of this quadrupole moment 
one must know the molecular electric field. This value 
can be calculated from the various wave functions which 
have been suggested for the hydrogen molecule. The 
result of such a calculation by Dr. A. Nordsieck with 
Wang wave functions when combined with our data yields 
a quadrupole moment Q=(32*—r*)y of about 2x 10~*’ 
cm*. The chief source of error lies in the inaccuracy of the 
wave functions. 

The sign of the quadrupole moment may also be inferred 
from our measurements in two ways. Present indications 
are that it is positive, that is, the charge configuration is 
prolate along the spin axis. Full details of these experi- 
ments will be published later in this journal. 

We are greatly indebted to Dr. Nordsieck for making 
available to us the results of his calculations and to Dr. 
Brickwedde of the National Bureau of Standards for 
preparation of the sample of pure HD used in the experi- 
ments. The experiments were supported in part by a grant 
from the Research Corporation. 

J. M. B. KeLiooe 
I. I. Rast 

N. F. Ramsey, Jr. 
J. R. ZACHARIAS 


Columbia University. 
Hunter College (J RZ), 
New York, New York, 
January 15, 1939. 


* Publication assisted by the Ernest Kempton Adams Fund for 
Physical Research of Columbia University. 

‘J. M. B. Kellogg, I. I. Rabi, N. F. Ramsey and J. R. Zacharias, 
Bull. Am. Phys. Soc. Vol. 13, No. 7, Abs. 24 and 25. 





Thermal Dilatation of Superconductors 


In recent years there has been much speculation as to 
the underlying cause of the superconducting state of 
metals, and many experiments have been performed to 
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investigate possible abrupt changes in other physical 
properties of materials in passing through the transition 
temperature. In 1925 Sizoo and Onnes' reported that the 
application of pressure lowered the transition temperatures 
for tin and indium, which indicated that metals may 
expand on passing into the superconducting state. This 
possibility was investigated by McLennan, Allen and 
Wilhelm? in 1931. They measured the thermal dilatation 
of lead and Rose’s metal in the neighborhood of their 
respective transition temperatures, but when they plotted 
the thermal dilatation against the temperature they ob- 
tained essentially smooth curves and, hence, concluded 
that there was no discontinuity in the thermal dilatation 
at the transition temperature. As recently as 1934* and 
1935* their work was still being quoted as providing evi- 
dence that no discontinuity exists in this region. As a 
matter of fact, when their data are plotted on a magnified 
scale, as illustrated in Fig. 1, we see that there appears to 
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be a discontinuity which is just within the range of resolu- 
tion of their apparatus. Since the magnitude of the ap- 
parent discontinuity is of the order of one unit in their 
last significant figure, their experiments cannot be ac- 
cepted as definitely demonstrating its existence. The mutual 
consistency of the two curves, however, does seem to indi- 
cate that the discontinuity is real, and makes it appear 
desirable to repeat the experiment with somewhat more 
refined apparatus and to obtain a larger number of points 
in the immediate neighborhood of the transition tempera- 
ture. From their report, the Toronto group seem to have 
pushed the sensitivity of the optical lever almost to its 
limit, however, there is the possibility of developing other 
methods such as one involving the displacement of inter- 
ference fringes for monochromatic light reflected from 
mirrors attached to the ends of a superconducting metal 
rod. 


E. C. WESTERFIELD 


University of Colorado, 
Boulder, Colorado, 
May 12, 1938. 


1G. J. Sizoo and H. Onnes, Communications from the Physical 
Laboratory of Leiden, No. 180b. 

2J. C. McLennan, J. F. Allen and J. O. Wilhelm, Trans. Roy. Soc. 
Canada 25, Sec. 3, 1-12 (1931). 

+E. F. Burton, Superconductivily (Toronto, 1934). 

*H. G. Smith and J. O. Wilhelm, Rev. Mod. Phys. 7, 240 (1935). 
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B-Radioactivities of Rhenium 


On bombarding rhenium with slow neutrons, Fermi and 
his co-workers' found a 8-activity with half-life 20 hours, 
and Kurtschatow and his collaborators* obtained 20-hour 
and 85-hour activities. Pool, Cork and Thornton* found 
that the bombardment of rhenium gives rise to an 18-hour 
period, which they ascribed to Re'™ obtained through 
neutron loss Re!’ (,2n) Re!®. ‘ 

We repeated the caperiments by bombarding metallic 
rhenium with slow neutrons (through paraffin) from Be + D 
and fast neutrons from Li+D (up to about 17 Mev) 
produced in the cyclotron in our laboratory. In both cases 
we obtained activities having half-lives 16 hours and 90 
hours, respectively, which correspond to the two periods 
above-mentioned. The ratio, however, of the initial intensi- 
ties corrected for the infinite activation time, or R= 
7(90 hr.) /J(16 hr.), was 1.2 for slow neutrons from Be+D 
and 17 for fast neutrons from Li+D. 

Now rhenium has only two isotopes Re!** and Re'*’. 
From the above large change in the ratio of the intensities 
for two cases, it seems therefore more probable, contrary to 
Pool, Cork and Thornton, to ascribe the 16-hour period to 
Re" and the 90-hour activity to Re™. 

Energy spectra of these two activities were studied with 
the Wilson cloud chamber method. Both emit electrons 
and the energy upper limits reduced from K-U plots were 
1.2 Mev for Re'™ (90 hr.) and 2.5 Mev for Re'** (16 hr.). 

A more detailed report will be published shortly in the 
Scientific Papers of the Institute of Physical and Chemical 
Research. 

We wish to express our best thanks to Dr. Y. Nishina for 
his valuab'e help and encouragement. 

K. SINMA 
F. YAMASAKI 
Nuclear Research Laboratory, _ 
Institute of Physical and Chemical Research, 
Tokyo, Japan, 
December 29, 1938. 


1 E. Fermi et al., Proc. Roy. Soc. 149, 522 (1935). 
? B. Kurtschatow ef al., Physik. Zeits. Sowjetunion 8, 589 (1935). 
(1957) L. Pool, J. M. Cork and R. L. Thornton, Phys. Rev. 52, 239 





The Effect of Hydroxyl Ion Concentration (or of pH) on 
the Film Potentials of Multilayers 


When an X multilayer of calcium stearate is built upon 
a subphase of metal, the resultant increment of film po- 
tential may be as high as 100 Mv per layer if the aqueous 
subsolution upon which the monolayer is spread initially 
is made alkaline to a pH of 9.4 by the addition of ammonia.' 
With potassium ions instead of ammonium ions the film 
potentials are lower, and of the order of 60 Mv per layer.? 

An experiment carried out by us in March, 1938, gave 
evidence of the effect of pH upon film potentials and sug- 
gested the more definite work described in this letter. In 
the earlier experiment the pH of the aqueous subsolution 
was brought to 9.4 by the use of ammonia, but this is 
volatile. As long as the pH was close to 9.4, the increment 
of potential per layer due to the taking up of positive ions 
by the outside of the multilayer was of the order of 85 Mv, 
but when the pH fell below this, the film potential also 
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dropped, and in some instances to zero. If after such a 
drop the pH was restored to 9.4, the increment of the film 
potential was immediately restored to a large value, as for 
example 300 Mv per layer. (The increase in magnitude 
above 100 Mv per layer is due to an increased thickness of 
the multilayer underneath.) 

In order to secure additional and conclusive evidence it 
seemed that X films should be built at a low pH if possible. 
An interruption in the work gave us the advantage that 
during this interval Goranson and Zisman discovered how 
to produce X films on ebonite at a pH of 6.3, and it seemed 
that their technique could be used for deposition upon a 
metal. 

In our work the speed of dipping was 5 cm per minute, 
with an interval of 30 seconds between each down and up 
trip but only a short interval when the slide was in the air. 

Several trends are evident in Fig. 1. 
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Fic. 1. Variation of film potentials with number of layers. ——-———— 
Y deposition; -— — —— Y changing to X deposition; — — — — — 
X deposition. The arrow marks the optical thickness at which the film 
gave an orange interference color with the electric vector of polarized 
white light perpendicular to the plane of incidence and with the angle 
of incidence about 80°. 


(1) The film potentials of the calcium stearate X multi- 
layers on gold at pH 6.3 do not increase linearly with the 
number of layers, but reach a value of about +0.1 or +0.2 
volt in a few layers, as with the Y type of deposition. 
After this the increment per X layer is almost zero (+). 
For similar X layers on ebonite at this pH, Goranson and 
Zisman report a charge which increases rapidly with the 
number of layers. 

(2) At a pH of 6.3, (10)-*m CaCl, gives the X type of 
deposition ; (10)-*m CaCl, gives the Y type of deposition. 
An intermediate molarity gives a deposition which may 
start with the Y type and change over to the X type. The 
contact angle changes continuously from about 0° to 
about 90° as the type of deposition changes from Y to X. 
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The increase in calcium concentration has no definite 
effect on the film potential at this low pH. 

(3) The experiments described in this letter, when com- 
pared with those at a pH of 9.4, prove that the pH of the 
solution is one of the primary factors which determine the 
film potential. 

The data of Langmuir and Schaefer* indicate that our 
monolayer at (10)~*m Ca** consists of 27 molecules of 
stearic acid to 73 of calcium stearate. We have used still 
higher concentrations of the ion in the solution to increase 
its proportion in the film. It is of interest to determine 
the relation between the concentration of the hydroxyl ion 
in the solution, the Ca** ion in the monolayer, and the 
film potential of the X multilayer. 

Witt1amM D. HARKINS 
RicHARD W. MATTOON 
George Herbert Jones Laboratory, 
University of Chicago, 
Chicago, Illinois, 
January 14, 1939. 


1 William D. Harkins and Richard W. Mattoon, Phys. Rev. 53, 911 
(1938). 

2 Eliot F. Porter and Jeffries Wyman, Jr., J. Am. Chem. Soc. 60, 1083 
(1938). 

3 Irving Langmuir and Vincent J. Schaefer, J. Am. Chem. Soc. 58, 
284 (1936). 





Chemica! Separation of Nuclear Isomers 


Many examples of nuclear isomers are now known, and 
it is of interest to have methods to find possible genetic 
relationships between the isomeric states of the same 
nucleus. For instance, this can be readily accomplished by 
standard methods if the lower level has a longer life than 
the upper. But for the opposite case the problem is much 
more difficult because a transient equilibrium is reached 
soon and the only apparent period is that of the parent. 

An actual separation of two genetically related isomers 
of a nucleus can be performed by the following artifice in 
some cases: Consider the case in which the element having 
the isomers can give compounds suitable for the application 
of the Szilard-Chalmers' method of concentration and 
prepare initially such a compound from a mixture of the 
two isomers. When the isomer in the upper state decays 
to the lower state the nucleus emits a gamma-ray and 
recoils. The energy of recoil may be sufficient to knock the 
decayed atom out of the compound and the daughter 
activity may then be separated as in the Szilard-Chalmers 
method. 

However, if the energy of the transition gamma-ray is 
of the order of magnitude of 100 kv, as one might expect,?~* 
the recoil for elements of medium and high atomic weight 
would not be sufficient to break an ordinary chemical bond. 
Fortunately the gamma-ray is highly internally converted 
and the recoil from the conversion electron is higher than 
that from the corresponding gamma-ray so that recoil 
energies comparable to chemical bond energies are ob- 
tained. 

Even if the recoil energy is insufficient to break the bond, 
the method does not necessarily fail because the energy 
given to the bond may suffice to create an activated condi- 
tion within the molecule such that it may undergo chemical 
reaction upon collision with some other kind of molecule 


THE EDITOR 321 


deliberately introduced for this purpose. The daughter 
activity will then be found in the product of the chemical 
reaction. 

We have used this method to separate the 18-minute 
and 4.4-hour isomers of Br®. The separation was fairly 
complete and showed definitely that the 18-minute period 
is formed from the 4.4-hour period. 

About one liter of ethyl bromide was irradiated for three 
hours with slow neutrons from the Berkeley cyclotron and 
the radioactive bromine was removed as solid sodium 
bromide. Radioactive tertiary butyl bromide was syn- 
thesized by saturating tertiary butyl alcohol with hydro- 
bromic acid generated from the sodium bromide. Five 
grams of the butyl bromide were introduced into 200 cc 
of aqueous methyl alcohol (75 mole-percent methyl 
alcohol) at 0°C where it underwent chemical reaction with 
both the methyl alcohol and the water liberating hydro- 
bromic acid. In this manner the butyl bromide molecules 
which contain a bromine nucleus which has recoiled can 
react to liberate hydrobromic acid. Under the experimental 
conditions chosen this type of activation occurs much 
more frequently than the spontaneous activation due to 
thermal agitation, and hence there is a large enrichment of 
the radioactivity of the lower isomeric state (18 minutes) 
in the reaction products. The chemical reaction was allowed 
to proceed for about three times the radioactive half-life 
of the lower state and then, after first extracting the 
unreacted butyl bromide with benzene, silver bromide was 
precipitated from the aqueous methy! alcohol. 

In Fig. 1, a, is shown the decay curve of such a concen- 
trate effected four hours after the end of the bombardment. 
The first reading was taken 27 minutes after the precipita- 
tion of the silver bromide. The activity evidently has two 
periods, a short one of 18-minute half-life shown in curve 
b, due to the recoil atoms, contaminated with an activity 
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Fic. 1. (@) Decay curve of bromine extracted from 4-hour activity of 


Br*®. (6) Curve a corrected for long-lived tail (18-minute half-life). (c) 
Decay curve of the bromine activity from which the extraction was 
made. 
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decaying in the same manner as the starting substance, 
whose decay is shown in curve c. This small long-lived tail, 
which we followed for 24 hours, arises from the bromide 
liberated in the spontaneous chemical reaction. By repeat- 
ing the experiment, but with an elapse of 26 hours between 
the end of the bombardment and the beginning of the 
chemical reaction, a decay very similar to curve a was 
obtained. These two runs also showed that the upper state 
(parent substance) was decaying with a half-life of ap- 
proximately four hours. By comparing curves 6 and c, 
and by taking into account that curve } started 27 minutes 
after the end of the chemical reaction, it can be seen that 
most of the nuclei in the upper state decay by falling to 
the lower state from which the observed disintegration 
electrons are emitted. 

DeVault and Libby have applied our method of separa- 
tion to other bromine compounds and will report their 
findings shortly. 

We wish to thank Professor E. O. Lawrence for his 
interest and encouragement in this work. Thanks are also 
due to the Research Corporation for continued support. 

E. Secret 
R. S. HALForD 


G. T. SEABORG 


Radiation Laboratory, 
Department of Physics, 
Department of Chemistry, 
University of California, 
Berkeley, California, 
January 13, 1939. 


1 L. Szilard and J. A. Chalmers, Nature 134, 462 (1934). 
2 W. Weizsiicker, Naturwiss. 24, 813 (1936). 
+ B. Pontecorvo, Phys. Rev. 54, 542 (1938). 
* E. Segré and G. T. Seaborg, Phys. Rev. 54, 772 (1938). 





Evidence for Gamma-Radioactivity of 4.5-hour Br” from 
Radiobromate 


In the course of work on the chemistry of the decomposi- 
tion of bromoform caused by the recoil from the neutron- 
capture gamma-radiation, it was noticed that the radio- 
active bromine persisted in coming out of the molecule 
over a long period of time after the irradiation ceased. This 
remained unexplained until Dr. E. Segré (cf. preceding 
Letter to the Editor) proposed a method of separation of 
nuclear isomers and that the 4.5-hour Br*® isomer was 
converted to the 18-minute form by emission of a soft 
highly forbidden gamma-ray. It was found that most of 
the activity in the Jater extracts decayed with the 18- 
minute half-life alf~hough no 18-minute bromine not grown 
from the 4.5-hour isomer could have been present in the 
bromoform extracted at the time. 

Simultaneously with the later experiments on bromo- 
form, work on the chemistry of the neutron-capture gamma- 
ray recoil for bromate ion gave evidence for the effect 
which was considerably more definite because it was 
possible to purify the BrO,;~ solution completely for 
bromine of lower valence by precipitation of AgBr from 
BrO;~ solution 2.5 molal in NH,OH. 

Figure 1A is a typical decay curve for a AgBr precipitate 
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Fic. 1. Growth and decay curves. A, decay of AgBr precipitate; B, 
growth curve for AgBr precipitates; C, growth curve for AgBrOs puri- 
hed for Br-; and D, decay of AgBr activity with 4.5-hour mother 
BrO;~. 


obtained from a mother BrO,;~ solution which had stood 
several hours so the 18-minute activity originally present 
when removed from the neutron source (200 mg RaBr.+Be 
powder) could have been present to less than five percent 
of the activity found. Apparently very little of the 4.5-hour 
activity was present in these precipitates. Figure 1B is a 
plot of the activities of AgBr precipitates obtained from a 
single mother BrO;~ solution which had stood various times 
since its last purification for Br~. The solid line is the theo- 
retical 18-minute growth curve. Figure 1C is a plot of the 
activity of a AgBrO, precipitate (obtained by acidifying 
the ammoniacal solution) which was purified for Br~ at 
zero time. The solid line is the 18-minute growth super- 
imposed on an initial activity of 160 counts per minute. 
This initial activity must have been due to 18-minute Br* 
still remaining in the BrO,~ after the activation by the 
kick from the electron ejected by the converted soft gamma- 
ray or to beta-radiation from the 4.5-hour body itself, or 
to both. However, the curve does show that unless the 
4.5-hour beta-rays are softer than about 500 kv (this was 
determined by the thickness of the sample and the wall of 
the counter) at least 75 percent of the 4.5-hour body must 
decay by the gamma-emission process. Further, in view of 
the fact that experiments with BrO,~ showed the efficiency 
of the much more powerful neutron gamma-ray recoil in 
destroying the ion to be only 80 percent approximately, it 
seems probable that most, if not all, of the initial 25 percent 
activity is due to retention of 18-minute activity in the 
BrO,;~ rather than to 4.5-hour beta-emission. Of course this 
indicates that the 4.5-hour body is practically a pure 
gamma-ray emitter, unless it has beta-radiation softer than 
about 500 kv. 

Figure 1D shows the activities of AgBr precipitates ob- 
tained at various times from a given mother BrO,~ solution. 
The solid line is the theoretical 4.5-hour decay curve. 

Don C. DEVAULT 


W. F. Lisry 
Department of Chemistry, 
University of California, 
Berkeley, California, 
January 14, 1939. 
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